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14. Abstract 

I 

These investigations aimed at a better understanding of the 
complex flov. field in high performance centrifugal 
compressors. Newly developed measuring techniques for 
unsteady static and total pressures (based on semiconductor 
transducers) as well as flow directions (special hot-wire 
probe), and a digital data analysis system for fluctuating 
signals (sampling technique) have been thoroughly tested. 
The loss-affected mixing process of the distorted impeller 
discharge flow has been Investigated in detail, in the 
absolute and relative system, at impeller tip speeds up to 
380 m/s. A theoretical analysis proved good coincidence 
of the test results with the DEAN-SENOO theory, which was 
extended to compressible flows. 
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1. Notation 

The meaning of the most importa*it symbols is given below. /lO 
VHiere deviations occur In individual cases, these will be 
explained in the text. 

a Sonic speed 

a^ Quiescent speed of sound 

A Area 

A/Aq Amplitude ratio 

b Diffuser width 

b’ Diffuser with ratio b/bg 

c Absolute velocity 

c^ Dissipation coefficient 

c^ Skin friction coefficient 

Cp Specific heat at constant pressure 

C Relative absolute velocity c/u2 

d Diameter 

D Diameter 

DS Pressure side 

e Base of the natural logarithm 

• 

E Kinetic energy flow of the absolute velocity 

f Relative total pressure fluctuation factor, Eq. (6.2) 

f Excitation frequency 

f Absolute total pressure fluctuation Intensity, Eq. (6.2) 

fjjj Mass flow correction value, Eq. (5.21) 

fj. Inherent frequency 

Pg Skin friction force in the Jet, Eq. (7.1a) 

P,p Skin friction force in the wake, Eq. (7.1b) 

h Ehthalpy 

h^ Total enthalpy 

ht Relative total enthalpy, Eq. (IV. 1) 

IDA Program for instantaneous data evaluation 

INTIDA Program for Integral, Inst^taneous data evaluation 
Loss relation, Eq. (6.5) 

KSO Cathode ray oscillograph 

1 Coordinate along the middle of the Impeller blade channel 
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Lr 

L2P 
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N 

Mr 

MR 
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N 

P 

Pref 

Pt 

Pt,rel 

T 

q 

r 

r 

R 

R 

R 

R' 

R* 

«K 

Rf 

Rg 

RE 


S 
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Length 

Specific work (on the fluid) 

Form changing work 
Drag work 

Two-focus laser, velocity measuring method 

Mass flow 

Sonic maoh number 

Quiescent mach number 

Torque required to overcome the forces of friction 
Mach number diffusion ratio, Eq. (5.55) 

Speed of rotation 
Number of measurements 
Power 

Static pressure 
Reference pressure 
Total pressure 

Relative total pressure, Eq. (6.6) 

Amount of heat added per unit mass 

Jet-wake pressure fluctuation. Pig. 7.1b 

Radial coordinate 

Interface radius of curvature 

Radius 

van Le mixing radius, Eq. (7.22) 

Gas constant fur air 
Simplified notation for R/R2 
Jet-wake mixing radius 
Channel radius of curvature 
Temperature recovery factor 
Reynolds number 
Unit of scanning 

Length of the outer contour of the Impeller through the 
vertical center section 

Jet-wake interface frictional force, Eq. (7.2) 

Summary Interface stability criterion, Eq. (8.4) 

Interface stability criterion under the effect of 
Coriolis accleratlon, Eq. (8.2) 
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Sj. Interface stability criterion xmder the effect of 

normal acceleration due to streamline curvature, Eq. (8.3) 

SRS Rear of Jet (In direction of rotation) 

SS Suction side 

SST Synchronized sanq>llng technique 

S/T Jet-wake 

SVS Front of Jet (In direction of rotation 
t Blade pitch 

t Time 

tg Gate opening time 

t^ Integration time 

t_ Measurement time 

tj, signal repitlon rate 

t_ Channel division In the probe measuring radius 

Tg Cathode ray transit time 

T Static tenqperature 

T Signal period 

T^ Total tenqperature 

TR Total temperature ratio, Eq. (6.4) 

u Circumferential speed 

u Electric a.c. component 

U Relative circumferential speed 

U Signal voltage 

U Electric d.c. component 

V Volume 

w Relative velocity 

W Referred relative velocity w/ug 

x/s Relative external contour of the Impeller In the 
“ vertical center section 

y Coordinate In the direction of the circumferential 

speed 

y/t Relative blade spaelnf 

y'/t* Relative blade spacing between the blade surfaces 

z Impeller blade number 

z Coordinate along the axis 

z/b Relative diffuser width 
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Absolute fi'j-' angle, relative to the positive clr- 
cumferenti - v. . ’ectlon 

Expansio factoi *’or the synchronized sanpling technique 

Relative flow angle, related to the negative circum- 
ference direction 

Blade angle 

Relative velocity ratio, Fq. (5.63) 

Boundary layer thickness 
Interference parameter, Eq. (5«64) 

Pressure norm correction factoi, page (7^) 

Boundary layer displacement thickness 

Difference, e.g. AT = temperature difference /13 

Relative total pressure change, Eq. (8.1) 

Wake surface component, Eq. (5.61) 

Continuity form factor, Eq. (5-57) 

Energetic form factor, Eqs. (5.58) and (5-59) 

Relative pressure fluctuation, Eq. (^.2) 

Interface coefficient of friction 
Internal efficiency, Eq. (4.3) 

Isentroplc Impeller efficiency, Eq. (5.67) 

Polytropic impeller efficiency, Eq. (5.68) 

Isentroplc Impeller total efficiency, Eq. (5.69) 

Cylinder coordinate 

Temperature norm correction factor, page (74) 

Isentroplc exponent 
Twist parameter 

Reduced output factor, Eq.(6.3) 

Wake - mass flow component, Eq. (5 <62) 

Circular constant 
Pressure ratio 
Total pressure ratio 
Density 

Density ratio p/p^ 

Phase relationship of the s8ciQ>les to the signal period 
Plow factor, Eq. (5.54) 
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^ Phase shift 

u Angular velocity 

Subscripts 

a Indicated value 

a Exterior 

A Initial value 

e Diffuser input region 

e Relative to the energy equation 

E End value 

HD Hot wire 

i Internal 

IDA From the program for Instantaneous data evaluation 

INTIDA From the program for integral. Instantaneous data 
evaluation 

k Relative to the continuity equation 

m Meridional 

max Maximal 

min Minimal 

N Standard atmosphere condition 

N Blade wake 

0 Design parameter 

r Radial 

red Parameter reduced to sta* dard condltic.is 
ref Reference parameter 

rel Relative 

s Isentropic 

B Probe 

8 Blade 

S Jet 

SST Snychronized sampling technique 

t Total condition 

T Wake 

u Circumferential component 

VENT Measured width Venturi tube 

0 Cong>re88or Intake > 
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1 

2 

3 

4 

5 
7 

10 

11 

13 

14 


II 


III 


IV 




Iripeller Intake 
Impeller exhaust 
Diffuser Intake region 
Diffuser, R/Rg ■ 1.687 
Compressor exhaust 

Two- focus laser measuring surface In the Inqpeller, 

x/S- ■ 0*08 
m 

Two-focus laser measuring surface In the Impeller, 

x/s„ - 0.43 
in 

Two-focus laser measuring surface In the Inqpeller, 


x/p 


n 


0.59 


Ti;o-focus laser measuring surface in the Impeller, 
x/s« ■0.87 


m 


Two- focus laser measuring surface In the In^eller, 


x/s 


m 


1.01 


Measuring surface In +ne diffuser Intake, 

R/Rg - 1.017 

Measuring surface in the diffuser intake, 

R/R2 - 1.039 

Measuring, surface in the diffuser Intake, 

R/Rg - 1.089 

Measuring surface In the diffuser intake, 

R/Rg ■ 1.151 

Turoulent or fluctuating compwnent 
Impeller 

Mean value, especially mean surface area, temporal or 
statistical mean value 

Center of mass; value calculated by means of h-s 
equations from representative mean values 

Related to the mixing radius 

Calculated by using ilq. (5*66) 




The equations are all homogeneous with regard to dimensional 
analysis, hence they a.-e valid in any coherent system of units. 
Dimensions are given in the SI system. 
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INVESTIGATION OF THE JET-WAKE PLOW OP A HIGHLY-LOADED 
CENTRIFUGAL COMPRESSOR IMPELLER 

Dietrich Eckardt 

Rheinisch-WestfSlische Technische Hochschule, Aachen 
2. Introduction /17 * 

High performance centrifugal compressors are today used 
primarily in mobile and stationary gas turbines — with units 
of power of about 50-1,000 kW — in turbo-superchargers and 
as process compressors in the chemical Industry. The primary 
aim of aerodynamic developments is, in general, to increase 
the efficiency scale. Thus, for small gas turbines, for 
example, with a 1? improvement in compressor efficiency, one 
counts on about a 3 % reduction in specific fuel consumption. 

Along with this, one tries to obtain high scale pressure 
ratios with a relatively large absorption capacity and a 
range of peT»formance characteristics as large and as stable 
as possible. 

The flow i.i the relatively long, narrow and three- 
dlmensionally twisted channels of a centrifugal compressor 
stage is considerably influenced by friction effects. Potential 
theoretical computer methods are not capable of predicting 
the actual flow with sufficient accuracy. Given the present 
state of the art, it .:eems that further, desired Improvements 
will only be possible if one is successful in obtaining basic, 
more profound knowledge of the actual flow conditions in 
heavily loaded centrifugal compressors. This will require a 
transition from the routine method used up to now — by 
necessity — of empirical corelation and extrapolation of 
existing stage designs. This would open up the possibility 
of modifying existing design procedures with simple but 
relevant model ideas — relevant in terms of the physical 
flow — in such a way that fairly optimal configurations would 
result. The precondition for fulfilling this aim is that it 
must be possible to accurately measure and analyze the com- 
plex centrifugal compressor flow under operating conditions. 

This study attempts to accomplish this in an important partial 
region of the flow. 

In the flow field of heavily loaded centrifugal impellers 
a number of flow events occur which are subject to friction. 

The significance of these individual events and their inter- 
action with each other is still largely unexplained (see section 
6 ,^. 2 ). As a result — mostly in the radial portion of the 
blade channel — three-dimensional flow separation occurs in /18 
the region of the suction side of the blades and the housing 


• Numbers in the margin Indicate pagination in the foreign text. 
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wall. Because of Its almost insignificant, relative kinetic 
energy, this region of separation is called the "dead water" 
region or wake and it increases up to the impeller exhaust. 

Thus, for heavily loaded radial impellers, it leads to the 
characteristic shift of the main flow — the "Jet" — in the 
region of the pressure side of the blades and the rear wall of 
the channel. 

This relative speed distribution of the jet-wake has con- 
siderable influence on the conversion of energy in the 
centrifugal compressor stage: 

- Inside the Impeller the theoretically desired flow 
slowdown is decreased, 

- immediately after the impeller exhaust, equalization of 
the disturbed flow is completed with dissipation and 

- finally, the instantaneous absolute flow resulting 

from the disturbed relative flow influences the pressure 
recovery in the following diffuser. 

This also outlines the central role which studies on jet- 
wake flow behind a highly loaded radial In^^eller have for /19 

detailed flow analysis of the entire stage. With exact 
measurements in this region it is possible to draw important 
conclusions on the Impeller flow. With these measurements — 
in conjunction with measurements of the stage efficiency — 
it is possible to make a separate loss analysis for the 
impeller and diffuser, and in so doing the mixing losses and 
equalization losses for the disturbed flow in the diffuser 
intake region can reasonable be charged to the impeller. 

Finally, exact knowledge of the developemnt of jet-wake 
equalization after the impeller exhaust is Important for the 
design of the subsequent diffusers and for the evaluation of 
this design in connection with standing questions on noise 
generation and on the stability of the range of performance 
characteristics for the compressor. The phenomenon of dis- 
turbed jet-wake impeller flow was recognized early on in the 
construction of radial engines. However, because of the 
well known difficulty in making measurements in a rotating 
system, this phenomenon was only rarely Investigated 
quantitatively. As early as 1902 Smith [2.1]^, under certain 
operating conditions, observed the appearance of a "partial 
filling" in the impeller of a radial pun^). The earliest, 


1. This information comes from an interesting bibliography on 
the topic of jet-wake flow which Moore included in his English 
translation of [2.2]. This was published as Technical Report 
CUED/A-Turbo/TR 73 (1975), University of Cambridge, England. 
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detailed measurements in a centrifugal compressor impeller include 
findings from the Brown-Boveri flow laboratory which were 
published by Garrard [2.2]1 in 1923* These results give 
clear references to Jet-wake distributions. Evidently Garrard 
was the first to try to consider theoretically the influence 
of the disturbed flow on the conversion of energy in the 
impeller using a simple Jet-wake flow model. In his termi- 
nology he referred to this as the "le Jet et la zone neutre” 

[the Jet and the nuetral zone]. After a number of further 
impeller studies at lower circumferential velocities (cf. 
section 6.4), it wasn't until I960, after all, before the > 

first quantitative studies were published on Jet-wake 
mixing behind a centrifugal compressor impeller. Pre- 
sumably, this delay was due to technical difficulties in 
making the necessary, exact measurements of the instantaneous 
absolute flow, both in magnitude and direction. Dean and 
Senoo C7»l]» based on a few quantitative hot-wire measurements, 
developed their two-dimensional theory for calculating the 
Jet-wake mixing under incompressible flow conditions. In 1974 
Senoo and Ishida [7*5] first published a Jet-wake flow analysis 
which was based on nonsteady measurements in several 
radii of the diffuser intake region. Two experimental studies 
[7.1, 7.5] were carried out behind centrifugal impellers with 
backward bent blades in the incompressible flow region at 
rotation speeds of n £ 5000 rmp. No results were given on the 
distribution of the disturbed flow over the diffuser width. 

0 until the time that the present study was begun, all 
quanticatlve investigations of the relative flow in centrifugal 
compressor impellers and of the absolute flow behind the 
impeller exhaust were limited to low rotation speeds and the /20 
region of incompressible flow because of the technical dif- 
ficulties mentioned above. At higher rotation speeds, measure- 
ments were only possible in the absolute system with conventional 
averaging probes . Their readings in the region of the non- 
steady flow strinlng the probes can be affected with con- 
siderable errors. 

In the first place, the aim of the studies presented 
here was to precisely measure and theoretically analyze the 
development of the disturbed, compressible flow behind a 
centrifugal compressor impeller with radial-ending blades 
under operating conditions approximating real conditions, i.e. 
with rotation speeds of up to n * 18,000 rpm (circumferential 
velocity at the impeller exhaust uo ■ 377 m/s). This would 
then provide sure data for improved methods for computing 
and designing high-performance centrifugal compressors. In 
order to achieve this program it was necessary to solve the 
following individual tasks: 


1. See fotenote, p. 2 
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1. Develop and test new, high-resolution (with respect to 
time) measuring and evaluation methods which are suitable 
for an exact study of a high level of rtonsteaay , aosolute ^ 
and relative flow. 

2. Experimental analysis and interpretation of the physical 
flow characteristics of the disturbed flow lines in the 
diffuser intake region: 

- in several points in the characteristic field of the 
compressor, 

- in several radially staggered measuring surfaces 
behind the imepeller exhaust, 

- with exact resolution in the circumferential direction 
over the flow region of a blade channel and over the 
diffuser width in the axial direction. 

3 . Expansion of the Dean and Senoo theory [7*1] i*or the 
region of compressible flows and testing the theoretical 
predictions on the basis of the two-dimensional jet-wake 
flow model by comparison with experimental results on the 
actual formation of the disturbed flow. 

3. Development of Measurement and Evaluation Methods for Studying /21 
Unsteady Flow Processes in Turbomachines 

We have already seen in the Introduction that the development 
of suitable measuring methods is. an important prerequisite 
for the analysis of the complicated flow conditions which exist 
at the exhaust of a highly loaded centrifugal compressor impeller. 

The experimental conditions in this region are distinguished by a 
number of technical difficulties involving measurements — dif- 
ficulties which, to a similar extent, are seldom met with in the 
experimental techniques of flow mechanics. Prom the outset, this 
considerably limits the number of fundamentally suitable measuring 
methods and evaluation methods. 

After extensive laboratory studies which were used create 
or further develop individual methods, a total measurir ; system 
was set up. With this system it is possible to solve the 
present problem — and also a number of other problems — in the 
study of nonsteady flow processes in turbomachines. 

3.1 Special Technical Reouirements for PI 'v Studies in 

Highly-Loaded Centrifugal Compressor St 'ges 

Consideration of the flow conditions after the Impeller 
exhaust reveals the highly nonsteady nature of rht flow in 
this region. This presents a considerable difficulty in making 
a precise and detailed determination of the flow condition. 

Essentially, the variations in the flow field over time — 
as observed in the absolute system — can be divided into three 





groups which are explained by means of Fig. 3*1* The figure shows 
a graphical representation of changes in pressure over time, as 
might be plotted, for example, with a wall pressure sensor placed 
in the coir^ressor housing Immediately behind the Impeller exhaust. 



Fig. 3«1 Changes in wall pressure over time to 
explain various forms of instantaneous flow in 
centrifugal compressors. 


The periodic variations observed in the stable operational 
region of the compressor — @ — have a basic frequency cor- 
responding to the product of the blade number and the rotation 
frequency of the impeller. The reason for these variations is /22 
the change in the condition of flow in the circumferential 
direction within an in?>eller-blade channel. This instantaneous 
flow pattern, in the case of a finite number of blades and 
friotlonless flow, represents the necessary condition for any 
form of energy transfer (cf. , for example, C3*l]» P* 209). Real 
flow effects — such as turbulence, boundary layers, secondary 
flow or flow separation and the manufacturing accuracy of the 
inq;>eller can also affect the variation in the condition of 
flow from blade to blade. These periodic variations also apply 
to the flow field behind the Impeller exhaust. Here they are 
balanced out by mixing and the exchange of energy after a cer- 
tain radial distance. 

An additional instantaneous phenomenon in the operating 
behavior of the compressor stage occurs when the mass flow is 
Induced below certain limit values. Rotating stalls or pumps 
® — are terms used to describe the partial or complete 
separation of the orderly compressor flow. The disturbance? 






caused by this in the absolute flow field can appear periodically, 
but their appearance can also be completely Irregular. 

The third type of Instantaneous flow shows up if we consider /23 
the changes in pressure for several successive blade passages. 

For example, the enlarged portion of the graph for the changes in 
pressure labelled © , which repeat themselves periodically with 
the blade frequency, have a generally similar curve which is 
superimposed by irregularly distributed Interferences. These 
deviations from a statistically average pressure graph are 
especially clear if the Individual signal characteristics are 
superimposed on the screen of an oscillograph by triggering with 
the blade frequency (3) . 

Irregularly occurring fluctuation magnitudes — referred to 
below as flow fluctuation or rough turbulence — are due to the 
overlapping of the general flow turbulence with nonsteady 
effects of the relative flow which, for example can sometimes 
be triggered by unstable regions of separation in the impeller- 
blade channels. 

For the studies made within the scope of this work in the 
stable range of characteristics, the statistical average Instan- 
tane^s flow pattern was determined in each case, as suggested 
by @ . In addition, a parameter referred to as a "fluctuation 
factor f" was determined. This describes the intensity of the 
rough turbulent fluctuations in the region of one blade spacing. 

A quantitative impression of the measuring techniques 
required for analyzing these nonsteady flow forms is given 
by the figures listed in the following table. Table 3»1» These 
refer to measurements made mm behind the exhaust of the 
centrifugal impeller studied with a diameter of 400 mm. 


A few fundamental requirements can be deduced from this 
concise presentation of the extreme test conditions at the 
exhaust of a highly-loaded centrifugal compressor impeller. 

The sensors and measuring methods must satisfy these requirements 
for detailed flow studied in this region. 


High Time flesolatlon 


For the example given in Table 3.1, the blade and /24 

thus the frequency of the periodically recurring flow pattern 
turns out to be 6 kHz. Suitable measuring systems must be 
capable of exactly plotting the change in flow parameters 
in the region of one blade spacing. In so doing, as shown in 
Fig. 3»le for the transition from the blade pressure side to the 
blade suction side, signal rise times can occur which exceed the 
blade frequency by a factor of 5-10. The inherent frequency of 
suitable measuring systems must then be at least 3 times greater 
than the test frequency, as, for example, shown in Fig. 3*3 for 
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TABLE 3.1 RANGE OF VARIATION FOR THE PERIODIC CONDITIONS OF FLOW 
BEHIND THE IMPELLER EXHAUST UNDER NORMAL CONDITIONS IN THE 
COMPRESSOR INTAKE (CP. FIG. 3.1e) 


Example ; 


Rotation speed 
Mass flow 
Blade number 
Signal frequency 


18 000 U/min 
7.^»0 kg/8 

20 

6 kHz 


Total pressure 

Static pressure ^*7 

Total temperature 591 

Static tenperature 338 

Static density 

Absolute velocity 337 

Absolute flow direction 9 

Total pressure fluctuation with 
respect to the local average 10 


4,3 

1.9 
448 
359 

1.9 
445 
31 

25 


bar 

bar 

K 

X 

kg/m^ 

m/a 

o 

% 


the example of the miniature pressure sesor. This Is to prevent 
Impexnnlsslbly high errors In the signs, anplltude and rlgnal 
phase relationships. For measurement!. In '^e centrifugal 
compressor design point considered above, t.ils yield necessary 
Inherent system frequencies on the order of 100-150 kHz. 

3.1.2 gbparat Ion of the Factors Affecting the Signal * 

Table 3*1 shows that all of the t^osslble parameters 
for describing the disturbed flow field at the Impeller exhaust 
vary considerably In the region of 1 blade spacing. Of the 
basically suitable measuring methods — suitable because of their 
high time resolution — most of them, because of their physical re 
latlonshlDB — address themselves simultaneously to several of the 
variable flow parameters. 


For example, the semiconductor strain gauge bridges of the 
pressure sensor can be put out of tune not only mechanically by 
compressive stress, but they also react to changes In ambient 



ten^erature. Likewise, the heat transfer equations for hot-wire 
anemometry Involve a large number of rapidly varying flow condition 
parameters in the exhaust region of a centrifugal Impeller. There- 
fore, a measuring method for the case under consideration can 
basically only be used if either: 

- the test signal depends only on the variable under study, 
or 

- the influence of other variables can be eliminated by 
special measures or 

- the signal variables can be separated by making supplementary 
measurements and by means of precisely known physical 
relationships. 

3.1.3 Small Probe Dimensions 

The requirement for high time resolution already par- 
tially contains the necessity of using miniaturized probes in 
order to reduce the slowness of the system. Extremely small 
probe dimensions are also necessary, however, so as not to 
falsify all too much the diffuser flow field which tends to 
separate. On the other hand, only probes with high mechanical 
resistance and corresponding dimensions possess the required 
service life for extensive measurements in the highly disturbed, 
roughly turbulent flow In the diffuser intake region. /26 

Therefore, studise on larger compressor stages represent a 
certain compromise on this point. In these larger stages the 
larger probe dimensions required for mechanical reasons stand in 
a favorable relationship to the stage dimensions. 

3>2 Review of Possible Measuring Methods 

With respect to the special requirements of the measuring 
technique discussed above, only a few measuring methods 
are suitable for detailed, quantitative flow studies behind a 
highly loaded centrifugal impeller. 

The following three measuring systems first of all possess 
the high time resolution mentioned above as a basic prerequisite; 

- miniature pressure transducer with high Inherent frequency 

- hot-wire and hot-film anemometer and 

- optical velocity measuring methods. 

The use of these measuring systems in various procedures also 
permits the building of sufficiently small probes -- the optical 
methods operate completely free of contact — so that, in prin- 
ciple, it is possible to record changes over time of the following 
flow parameters: 

- wall pressures, 

« total pressures, 

- velocities, magnitude and direction, and 

- flow fluctuations. 



In the last analysis, the most critical criterion for the 
usability of a measuring technique is the test as to whether 
the Initial test signal Is affected by only one condition-of- 
flow variable or whether possibilities exist for separating the 
signal variables. In the case of the pressure transducer and /27 
hot-wire measuring technique, extensive laboratory tests were 
required to answer this question. These will be discussed In 
the appendix. 

The optical velocity measuring methods will not be included 
In the further discussion, since these methods were not yet 
available at the time of the experimental Investigations behind 
the centrifugal Impeller, which is the object of this study. 

In the meantime, using a "dual- focus laser method" developed 
by Schodl [3.2] especially for measurements in turbomachines, 
it has been possible to record velocity distributions Inside the 
blade channels of this impeller (cf. section 6.4.1). Besides 
the possibility of non-contact measuring, the essential ad- 
vantage of this method consists In the fact that. In principle, 
the measured variable is not falsified by other flow variables. 

3.3 Method of Instantaneous Pressure Measurement 

The measurement of rapidly varying pressure changes observed 
in the absolute system within the sphere of influence of a 
turbomachine inpeller requires measuring systems whose Inherent 
frequencies are four times greater than the blade frequencies. 
Miniature pressure transducers on a plezoquartz or semiconductor 
base fulfill this condition. 

In the first case, by means of the underlying Pezo effect — 
release of an electric charge during deformation of a quartz 
crystal due to compressive forces — it is only possible to 
measure the dynamic pressure components. Pressure transducers 
which operate on this principle have Inherent frequencies of up 
to 400 kHz. To be sure, there are limits to how far they can 
be miniaturized because of increasing difficulties in the 
processing of small electric charges. The currently available 
plezoquartz pressure transducers have a pressure-sensitive front 
surface of about 5 nun in diameter. Thus they are unsuitable 
for insertion into Pitot tubes. In addition. In measuring in- 
stantaneous changes In wall pressure, the large membrane diameter 
is a drawback for the desired measurements which, as much as 
possible, should be point measurements. 

By contrast, with the advances made In the miniaturization /28 
cf electric components In recent years, it has been possible to 
develop strain gauge bridge pressure transducers based on semi- 
conductors which present themselves for use In turbomachines. 

This is because of their small dimensions — in general, they 
have a diameter between 1.2 and 3 mm. 
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After detailed laboratory studies had produced Information 
on the operational behavior of these transducers, measuring 
systems were developed from them to make Instantaneous wall 
pressure measurements and total pressure measurements. These 
systems take Into account the special requirements for use In 
highly-loaded centrifugal compressors . 

3.3«1 Operational Behavior of Miniature Semiconductor 
f>ressure Transducers 


The extremely small dlmesnlons of these pressure 
transducers are achieved by means of a special semiconductor 
photomasklng technique by which a complete, 4-arm Wheatstone 
bridge Is vapor deposited on the pressure-sensitive silicone 
membrane. The ”blocklng” p-n transfer (between positively and 
negatively doped semiconductor material) resulting In this way 
obviates separate Insulation of the Individual resistor elements. 

Fig. 3*2 shows. In principle, that when the membrane Is put 
under compressive stress, the parallel arrangement of the 
resistor elements leads to variable deformation (elongation or 
compression) of the Individual resistors and thus puts the 
bridge out of tune. 

The following points are especially Important In the use 
of these pressure transducers for flow studies on turbocom- 
pressors: 

- the frequency resolution and the acceleration sensitivity 
of the transducer, 

- the Installation conditions and 

- the Influence of thermal effects on the pressure trans- 
ducer output signal. 

3. 3 « 1.1 Frequency Resolution 

Miniature pressure transducers are a nearly /29 

undamped oscillating system, comparable to a simple spring-mass 
coupling, which is excited to forced oscillations by external 
pressure forces. In this process It Is aimed at to keep the 
mass of the membrane and the bending of the membrane under 
stress as small as possible. 

In the case of the pressure transducers (made by Kullte, 
basic type CQL-080-25, outside diameter 2 mm) for these centri- 
fugal conqpressor studies, the actual pressure-sensitive part of 
the membrane with a diameter of 1.2 mm Is about 0.02 mm thick 
and at the maximum allowable pressure load It Is deflected about 
0.007 mm. At the same time, the density of the silicon membrane 
Is less by a factor of 4 than that of the usual steel membranes 
(cf. C3.3, 3.^13). 

This combination of a small membrane mass and a high spring 
constant determines the very small mechanical acceleration sen- 



Pressure transducer 



Si -substrate 

p-doped (GA) 

,pn -transition 
e- doped (Sb) 




Pig. 3*2 Design and principle of operation of the 
pressure-sensitive membrane of a miniature semi- 
conductor pressure transducer. 


sitlvity and the high mechanical Inherent frequency of these 
pressure transducers — f » 150 kHz — as can be seen from the 
principle frequency curve shown in Pig. 3* 3* Prom this curve 
of the amplitude ratio A/A^, i.e. the phaae shift over the 
excitation frequency f, it can be deduced that with this 
measuring system even rise times corresponding to an upper test 
frequency of fmax ^ 0.33 fp measured without impermissibly 

high falsifications. 

In this connection it seems Important to point out a few 
results of dynamic pressure calibrations. Weyer C3*5] checked 
the frequency behavior of these semiconductor transducers with 
a specially developed pressure generator which produced trape- 
zoidal pressure curves up to a frequency of about 4 kHz and a 
doubled pressure amplitude of 0.8 bar. Kurtz and Kicks C3>3] 
reported on similar experiments in which small sinusoidal 
pressure variations up* to 20 kHz were produced with a loud 
speaker. In both studies, deviations in the dynamic pressure 
transducer sensitivity of less than tl% with respect to the 
static value were found. 

3. 3.1.2 Influence of the Installation Conditions 

With respect to the requirements already mentioned 
in section 3.1.3 for small probe dimensions and measurements as 
close to "point" measurements as possible, it was tested to what 
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Pig* 3*3 Principle frequency 
characteristics of a miniature 
semiconductor pressure trans- 
ducer « f ■ excitation fre- 
quency) A/Ap ■ anq;>lltude ratio, 
^ ■ phase shift. 


extent It Is possible t.o set 
the transducer separate "rom the 
actual measuring site. 


The results of simple 
estimates for the "organ pipe 
arrangements" shown In Fig. 

(upper right)) l.e. the 
Helmholtz rebonator, show, how- 
ever, that Inherent system 
frequencies are to be expected 
with this arrangement on an 
order which is considerably 
less than those of miniature 
transducers. Conversely, It 
follows from this that the 
dynamic characteristics of these 
transducers can be taken advan- 
tage of for de 3rmlnlng rapid 
changes in pressure only when /32 
the pulsating pressure acts 
directly on the membrane. In 
accordance with this. Fig. 3*5 
shows favorable Installation 
arrangements for the pressure 
transducers : 

- for wall pressure 
measurements the pressure trans- 
ducer membrane should be in as 
close contact as possible with 
the housing wall; 

- for total pressure 
meCLSurements It is advisable to 
set the membreme cbout 0.5 - 
1 d (d '' probe head diameter) 
behind the Input opening for the 
probe in order to achieve 
sufficiently high direction 
sensitivity for the probe. This 


point becomes especially inqportant In measuring the nonsteady 
total pressure directly behind the centrifugal Impeller exhaust 
because of the large fluctuations In direction which occur 
there (see Table 3*1)* 

3. 3. 1.3 Thermal Effects 

In the layout point of the centrifugal compressor 
stage in question mean total temperatures \xp to about 430 K 
occur behind the Impeller exhaust (see Table 3*1)* Because of 
the necessary type of construction for the probes described 
In the preceedlng section, the semiconductor elements In the 
pressure transducer membrane are also exposed to temperatures 
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Pig. 3.M Inherent system frequencies fj. of various 
measuring arrangements with a clearance volume in 
front of the pressure transducer membrane for air 
or helium (He) in the standard state, L » tube length, 
A ■ tube surface area, V ■ volume 


Pressure transducer 



Pig* 3.5 Principle Installation arrangements of 
Kulite pressure transducers for wall pressure 
measurements and total pressure iMasurements . 


on this order. Therefore it is particularly Important to check 
the effects of temperature on the pressiu:*e transducer reading. 

With the various temperature effects on the output signal /33 
of semiconductor pressure transducers, one distinguished between 
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Pig. 3*6 Original calibration 
curves for a miniature semi- 
conductor pressure transducer 
(total pressure probe) at vari- 
ous ambient ten^>eratures. Feed 
voltage: 5 volts, electrical 
anqpllfication factor: 40. 

Key: A) 

A) Pressure transducer 

B) Initial voltage 

C) And^lent temperature 

D) Excess pressure 


studied, the tenqperature-related 
than -0.7S/10 K. These findings 


the following according to 
existing U.S. test conditions 
(cf. [3.6, 3.73): 

- changes in pressure trans- 
ducer sensitivity (rise in the 
calibration curve) due to 
temperature effects, 

- temperature-r#:lated zero 
point drift and 

- effects of temperature 
gradients. 

In extensive laboratory 
tests, whose arrangement and 
procedure are described in 
appendix I, the operational 
behavior of the miniature 
pressure transducer was precisely 
investigated under conditions of 
variable ambient temperature. 

The following figures will 
elucidate the most important 
findings of these experiments 
and the conclusions deduced fro:.i 
them with respect to the method 
finally chosen to make non- 
steady pressure measurements. 

Fig. 3.6 shows pressure /3^ 
calibration curves for a Kulltc 
pressure transducer at 4 ambient 
temperatures . These were 
determined by using a U-oube 
manometer filled with mercury. 

With Increasing temperature the 
sensitivity of the pressure 
transducer steadily decreases. 

In this example, at 140® C it 
is smaller by ab ut 3 % than at 
20® C. In general, wits all of 
the Kullte pressure tx'ansducers 
changes in sensitivity were less 
were quite reproducible. 


The temperature-related zero point drift is expressed in 
the shift in pressure calibration curves to higher voltage values 
observed with increasing temperature. In the exan^le shown, the 
zero point drift between 20®C and 140®C is on average about 
+1.8J f.s./XO K. To be sure, large differences occurred from 
experiment to experiment and from one pressure transducer to 
another, so that in apit«i of considerable efforts we did not 
succeed satisfactorily in setting up mathematical relationships 






for the changes in zero point drift with respect to temperature. 


This can certainly be explained in part due to the very 
large sensitivity of response of the semiconductor elements with 
ryspect to temperature gradients. As a characteristic result 
of laboratory tests (also see appendix 1) which serves to clarify 
this question, Fig. 3*7 shows the changes in output signal over 
time for an unloaded pressure transducer. It was exposed to 
sudden jumps in ambient temperature. In so doing, the detuning 
of the semiconductor bridge reached a maximum value after about 
1 minute. At AT = 19 K this is about 14/^ of the nominal pressure 
reading (25 psl, 20®C). 



Fig. 3*7 Changes in zero point drift over time for 
an unloaded Kulite KQL-093-25 pressure transducer 
under the influence of temperature gradients, 

U * initial voltage at normal load (25 psl) and 
To = 298 K. 


Tests with a water-cooled pressure transducer adapter 
(cf. Fig. 3-9)» as was developed for wall pressure measurements, 
on the other hand produced extensive damping of the temperature 
gradient effects. Fig. 3*8 compares the detuning of cooled and 
uncooled pressure transducers 1 minute after a rise in temperature 
At. 


The effectiveness of water cooling in suppressing temperature 
effects was also found in long-term tests and was tested in 
practice in the compressor. While blowing in a Jet of hot air 
at a temperature of 150°C, the temperature on the pressure trans- /35 
ducer membrane rose only about 6 K with respect to the 
temperature of the cooling water, as in nonsteady wall pressure /36 
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Pig. 3*8 The zero point drift of cooled and un- 
cooled semiconductor pressure transducer under 
the influence of temperature gradients, test time 
t = 1 minute, U„ = initial voltage at normal load 
(25 psi) and Tq = 298 K. 

Key: A) Ambient air D) Uncooled pressure 

B) Initial temperature transducer 

C) Rise in temperature E) Water-cooled pressure 

transducer 


measurements in the exhaust region of the centrifugal impeller 
(up to 18,000 rpm) the deviation from the exact value of the 
mean pressure transducer reading (measured against atmosphere) 
was in the region of ±4^. 

3.3.1 .^ Results of the Studies on Operational 
Behavior 


In summary, the following was found after the 
laboratory testing of the miniature semiconductor pressure 
transducers : 

- very frequent pressure variations can be measured with 
sufficient accuracy up to 1/3 of the pressure transducer 
inherent frequency if the installation conditions are 
chosen in such a way that the pressure pulsation acts 
directly on the membrane; 

- the pressure transducers ;'ermit exact measurement of the 
dynamic wall pressure or total pressure component (pul- 
sating pressure changes) if the change in the transducer 
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sensitivity with ambient temperature Is taken Into account 
In the evaluation; 

“ regularities in the zero point drift behavior of the 

transducers cannot be detected with any reliability. With 
water cooling the falsifying temperature effects can be 
considerably damped, however even with this measurements 
of the absolute pressure level with the desired accuracy 
are not possible. 

On the basis of these conclusions a measurement and 
evaluation method was developed with which the nonsteady absolute 
pressure curve is reconstructed by superimposing separate 
measurements of the corrected pressure mean value over time and 
of the dynamic pressure component. This method is discussed 
in detail in section 5.3*1» while the individual measuring 
systems for nonsteady or average wall pressure and total pressure 
measurement over time are discussed in the following two sections. 

3.3.2 Development of Measuring Systems for Determining /37 

Changes in the Dynamic Wall Pressure and Total Pressure 

The basis of the nonsteady pressure measurement systems 
is formed by various types of miniature semiconductor pressure 
transducers made by Kullte. The most important operational data 
for these transducers is summarized in Table 3.2. In addition /38 
to the manufacturers' specifications, our own results obtained 
from the laboratory tests on the pressure transducers used on 
the centrifugal compressor are also listed. 

T"'/e initial voltages of the pressure transducers were 
amplified by means of a differential d.c. ampl.^'fler with a 
linear frequency curve up to about 100 kHz. This was to facilitate 
the further processing of test signals by means of oscillographs 
or the synchronized sampling technique (see section 3*5) • 

The findings from the tests on the operational behavior of 
the semiconductor pressure transducers were taken into consider- 
ation in designing the nonsteady pressure measurement systems. 

As shown in Pig. 3.9, for the wall pressure measurement a 
cylindrical pressure transducer, Kullte XQL-093-25, was installed 
flush with the housing wall in a cooling adapter. Thus the fre- 
quency resolution of the transducers is fully taken advantage of 
and at the same time temperature effects on the output signal of /39 
the transducers are largely suppressed. As a result of the water 
cooling, the transducer, whose nominal operational range extends 
up to 120®C, can be used up to substantially higher temperatures. 

The preliminary tests also showed that with cooling, the 
temperature on the membrane increases only slightly even at a 
maximum temperature load in the centrifugal compressor (18,000 
rpm) . Thus it was possible to regularly use only the 20®C 
calibration curves for evaluating the dynamic wall pressure curves. 
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TABLE 3.2 OPERATIONAL DATA ON KULITE ::QL-093“25 AND XBL-9-093-25 G 

SEMICONDUCTOR PRESSURE TRANSDUCERS 


Pressure differential measurement 
range (against atmospheric pres- 
sure) 

non. 25 psi s 
1.72^» bar 
♦ 20 

Sensitivity 

nom. 36,25 mV/bar 

at 


feed voltage 

5 V - 

Temperature 

XQL: -26 ®C to 120 ®C 

operating range 

XBL: -20 ®C to 175 

Nonlinearity and 

*) 

max. ♦ 1 X f.s. * 

hysteresis error 

Change in sensitivity with 

Kulite; 0 to -iJ/io K 

with Temperature 

DPVLR; - 0,25 to -0.7JS/10K 

Zero point drift with Temp. 

Kulite: ♦3%f .8.+0.36*f .s./lOK 

(Shift test bed, see App. I) 

DPVLR : max. ♦ 2$f.s./10K 

Zero point drift due to 

DPLRV: see Pig. 3-7 

temperature gradients 


Inherent frequency 

approx. 150 kHz 

Acceleration 

Perpendicular : 0.3^5 Pa/g ' 

sensitivity 

horizontal : 0.034 Pa/g 

•) % f.s.: in percent of the nominai test range (full scale) 

••) g; acceleration due to gravity 




A 


Kuhlwautr 
CO. 0L6 I /min 


KUUTE 
XQL. 093-25 









Schutzsicb vor 


Q Oruckoufnthmcr • 
Mtmbron 


Fig* 3*9 Cooling adapter for measuring changes In 
(dynamic wall pressure with a miniature semiconductor 
pressure transducer, measurement In mm. 

Key: A) Cooling water about 0.6 liters per minute 

B) Protecting screen in front of the transducer 
membrane 


This considerably restricted the calibration costs for the 
l6 wall pressure transducers in simultaneous use. Moreover, 
even after longer operation times, no aging effects were detected 
in the calibration characteristics. 

For measuring the dynamic changes in total pressure the 
hook-shaped probe, shown in Fig. 3*10, was designed and the 
head of the probe was completed with a built-in semiconductor 
pressure transducer according to DPVLR specifications by 
Kulite. The transducer was designated as XBL-9-093-25 G. By 
directly Integrating the transducer in the manufactured pro- 
duct, a favorable probe Intake diameter of 2.2 mm resulted in 
relation to the exhaust width of the centrifugal impeller (26 mm). 

Effects of the probe shaft on the flow field to be studied 
are avoided by the probe head length of l4.5 mm. The design of 
the probe Intake part (see Fig. 3.10, upper right) enables exact, 
direction-independent total pressure measurements in an angular 
range of 30® (cf. [3*8], pp. 10, 11). Our own tests with 
steady and nonsteady oncommlng flow — in a wind tunnel up to a 
mach number M ■ 0.5 and with the probes mounted behind the cen- 
trifugal Impeller — confirms this data given in the literature. 
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Pig. 3*10 Pitot tube with built-in miniature semi- 
conductor for measuring changes in dynamic total 
pressure, measurements in mm. 

Key: A) Device for adjusting the probe 

B) Membrane 

C) Protecting screen 


To protect the silicone membrane of the pressure tr^xnsducer 
from dust erosion a perforated screen was installed in the head 
of the probe on the tip. In studies by Kurtz and Kicks [3*3] it 
was found that this protective screen (also shown in Fig. 3*9, 
perforated form covering the entire area) only insignificantly 
affects the frequency characteristics up to 30 kHz of the 
transducer with variations in pressure. 

The pressure transducer built into these probes can be 
operated without cooling up to a temperature of 175®C. In order 
to correct the temperature effect on the sensitivity of the 
transducer pressure calibration curves were plotted between 20°C 
and lilO®C — after artificial thermal aging of the transducers 
(see Appendix I). These calibration curves were used as the 
basis in evaluatl’‘'g the dynamic total pressure curve in accordance 
with the mean temi ex’ature load in each case. Subsequent cali- 
brations of the probes after being in operation for 25 hours at 
temperatures in excess of 100® C showed no measureable changes in 
the probe characteristics. 

3.3.3 Measurement Systems for Determining Exact Pressure 
Mean Values over Time ~~ 


Simple pneumatic measurement systems for determining 
pressure mean values over time — e.g. the following arrangement: 
Pitot tube or static wall pressure hole — test line — slow 



manometer — can produce substantial erroneous readings when 
used in regions with quickly changing pressures and high pressure 
amplitude (e.g. within the sphere of influence of a turbo- 
machine Impeller). In studies on 'transonic axial compressor 
stages variations up to 10 % of the correct pressure mean value 
over time were found (cf. [3*5]). /^1 

This measurement uncertainty is influenced by a large number 
of parameters, the most Important of which are the following: 

- the ratio of pressure an5)litude to the mean pressure value, 

- the shape of the pressure curve over time and 

- the shape of the probe Intake or hole intake. 

The problem of determining mean pressure values over time 
in turbomachine has been extensively discussed by Weyer [3*5]. 

In the scope of these studies 3 different measuring methods were 
developed. Two of these methods are described below which were 
used in the centrifugal compressor for measuring correct mean 
wall pressure or total pressure values over time. 

3. 3. 3.1 Special Pneumatic System for Wall 
Pressure Measurements 


The difference between the indicated and correct 
mesui pressure value over time in the case of the usual pneumatic 
systems can be attributed to inflow and outflow phenomena of the 
gas in the measuring system during a period of variation in 
pressure. 

Starting with the laws of mass conservation and energy 
conservation (cf. [3*53) » it is possible to calculate this 
difference for a special measuring system in which the conditions 
present can be readily grasped in physical terms because of the 
shape of the measuring hole and the choice of the cross-section 
ratio of the measuring hole to the connecting line. However, 
the dynamic pressure curve (without reference level) must of 
course be known in order to make the correction. 

Fig. 3.11 shows the principle construction of this special 
pneumatic s-stem. The measurement hole is in the shape of a 
semicircular nozzle. In this way, nearly constant and uniform 
flow coefficients are obtained for the inflow and outflow phe- 
nomena in the interesting Reynolds number range. The cross- 
section of the measurement nozzle is very small in relation to 
the connecting measurement line, so that pressure and temperature 
fluctuations within the system can be ignored. 

Moreover, the axial extent of the nozzle is very small in /^2 
relation to the wavelength of the pressure fluctuation, so that 
the flow phenomena in the region of the measurement hole may be 
regarded as quasi-steady. 
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Fig. 3 >11 Special pneumatic measuring system for 
determining exact mean wall pressure values over 
time, measurements In mm. 


For recording the mean pressure value determined with this 
system (which Is still corrected mathematically) a common, slow 
pressure measuring Instrument can be used (U-tube manometer, 
pressure transducer with a low Inherent frequency). Because 
of the necessarily large cross-section of the connection line 
connected to the measurement nozzle, this system is suitable 
only for wall pressure measurements, bat not for total pressure 
measurements in the narrow flow channels of the centrifugal 
compressor stage. 

3.3»3.2 Hydraulic System for Total Pressure 
Measurements 


The principle construction of thls^easuring 
system Is shown in Fig. 3.12. A capillary tube ^ 0.35 mm In 
diameter and 50 mm Ijmg Is rigidly attached to the pressure 
transducer adapter ^ . A pressure transducer with a low 
inherent frequency 0 serves as a signal modulator. The total 
pressure varying over time acts on the open end of the capillary 
tube at the head of the probe (|) . The pressure pulsation is 
damped by the viscosity forces acting In the capillary tube so /A 3 
that nearly constant pressure is adjacent to the pressure trans- 
ducer 0 . 

The entire system Is filled with silicon oil whose kinematic 
viscosity at room tenqperatures Is about 1.5 cm2/s. During the 
pressure measurement , a small amount of oil (about 100 mmVhour) 

Is put In at point Q) In order to make sure that the capillary 
tube Is always filled up to the Input point. 

The measurement principle Is based on the Impulse equilibrium 
between the almost constant pressure over time on the pressure 
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Pig. 3*12 Hydraulic measuring system for determining 
exact mean total pressure values over time, measure- 
ments in mm. 

1 » capillary tube, 2 * pressure transducer adapter, 
3 ® oil input, k » pressure transducer with low in- 
herent frequency, 5 ■ probe head with latteral holes 
for oil outflow, measurement in mm. 


transducer membrane and the integral value over time of the 
variation in pressure at the probe intake. To determine- the 
exact value two corrections are generally necessary : the first 
concerns a possible difference in height between the test point 
and the pressure transducer membrane, while the second concerns 
the slight drop in pressure of about 5 mb which is caused by 
the necessary flow of oil in the capillary .tube. 

Because of the small dimensions of the probe head CE)» 
whose sensitivity is nondirectional up to ±30JJ, the hydraulic 
method is especially well suited for determining the correct 
mean total pressure value over time of the strongly fluctuating 
flow at the exhaust of the centrifugal Impeller. The accuracy 
of the total system depends basically on the pressure transducer 
used. With this system, even when used at high temperatures on 
the centrifugal compressor with appropriate water cooling of the 
pressure transducer, reliable results were obtained. 

3.M Development of a Method for Measuring Direction in a 
rerlodlcally Unsteady l!*low * 

3.4.1 Problems in the Use of Hot-wire Anemoiaetry 

For detailed studies of the strongly disturbed flow 
field at the exhaust of highly-loaded centrifugal compresser 
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in?)eller8, precise information is necessary on the distribution 
of the absolute flow direction. Following the successful 
testing of the miniature pressure transducers in measuring 
rapidly changing wall pressures and total pressures over time, 
it was obvious to use the high frequency resolution of hot-wire 
anemometry to measure other fluctuating flow variables, such as 
the magnitude and direction of flow velocities and total 
tenperatures. 

A glance at the mathematical relationships of heat transfer 
on thin, freely stretched hot wires (cf. for example, [3*93» 
p. 75) shows, however, that direct measurement of the flow 
parameters mentioned above is only possible in a few experi- 
mentally Idealized cases without being falsified by other 
Influences. 

The following are the important variables which affect a 
hot-wire test signal: 

- Reynolds number, 

- mach number, 

- Prandtl number, 

- flow direction, 

- degree of turbulence, 

- over heating rate, /^5 

- dust deposit and aging, and 

- Isentroplc exponent of the flow medium. 

Of these variables at least 3 occur in the impeller exhaust 
flow of a highly loaded centrifugal compressor, sometimes with 
large, rapidly changing fluctuations. Since each probe murt 
be calibrated individually in the total range of variation for 
the different variables, it very quickly becomes clear that 
the required cost of calibrations is out of proportion with the 
relatively short life span of hot-wire probes. 

From this list of fundamental problems involved in the use 
of hot-wire anemometry, it follows that for centrifugal compressor 
studies only those hot-wire techniques can be used with any 
prospects for success in which the number of variables affecting 
the test signal is reduced by special measures. A measuring 
technique for measuring nonsteady directions developed on the 
basis of this requirement is discussed below. This uses only the 
high time resolution of hot-wire anemometry, but with which the 
difficult lnteri>retation of signal amplitudes is not applicable. 

A measuring technique of similar design for quickly changing 
total ten?)eratures is discussed in C3«10» 9. 13* This has not 
yet been tested on the centrifugal compressor, however. 

j. 4. 2 Construction and Measuring Principle of the 
Hot-wire Direction Probe 


The design of the probe is elucidated by referring to 
Pig. 3 . 13 . A tungsten wire 0.6 mm long and 5 ym in diameter is 




Pig. 3*13 Head of the hot-wire shielded probe for 
measuring the nonsteady flow direction, measurement 
in mm. 

Key; A) Hot-wire diameter 5 um 

B) Shield 

C) Direction of oncoming flow at the test point 


stretched parallel to the axis of the probe. The shaft of the 
probe with a diameter of 2.5 mm is formed into an umbrella- 
shaped screen in the vicinity of the probe head (135® cylinU 5 r 
segment). The sharpened leading edge also lies parallel to the 
axis of the hot wire and the probe. The probe is connected to 
a Dlsa 55 D 01 constant temperature anemometer by means of a 
30-meter high frequency cable and is operated in a modified 
1:1 bridge circuit. The overheating rate of the wire is 
selected such that the temperature of the wire is close to the 
upper load limit of 300®C, thus ensuring the greatest band width /46 
and sensitivity of the measuring system. 

Fig. 3*14 helps clarify the ^Inclple on which the measuring 
technique is based. At position © it showed schematically the 
characteristic, disturbed relative flow at the blade channel 
exhaust of a highly-loaded centrifugal Impeller. An idea of the 
real velocity curve is given by the broken line whi^h is 
ajaproxlmated by the Jet-wake velocity block profiles. Position 
© shows the corresponding velocity triangle for this Jet-wake 
distribution. One can clearly see that the Irregu"* ■'.r changes in 
relative velocity in the circumferential direction *.n the 
absolute system are expressed in a change in the flow direction a 
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Pig. 3*1^ Principle of non- 
nteady direction measurement In 
the disturbed flow at the Im- 
peller exhaust using a hot-wire 
shielded probe. 

Key: A) Jet 

B) Wake 

C) Hot-wire 

D) Probe shield 

E) Direction of rotation 

F) Oscillogram 
0} Time 


between the extreme values as 
(for the Jet) eind oqi (for the 
wake ) . 


This variation In direction 
of the absolute velocity vector 
over a blade spacing t between 
the pressure side DS and the 
suct^n side SS) Is diagramed 
at (^. As shown here, the 
hot-wire shielded probe is 
mounted In such a way immediately 
behind the impeller output so 
that the axis of the probe 
lies parallel to the rotor axis 
and the alingment line between 
the leading edge of the shield 
and the straight line between 
leading edge of the shield and 
the hot wire forms a ’’probe 
adjustment angle” oy with the 
circumferential direction. 

If the wake region In the /^7 
vicinity of the suction side SS 
passes the probe with a<ay 
during the rotation of the Impel- 
ler, the oncoming flow strikes 
the hot wire directly and a high 
signal voltage is observed on the 
oscillograph during the time t-y, 
position (5) • passing the 
Jet region with a>ay the hot- 
wire Is shielded from the direct 
oncoming flow. It l.^es in the 
wake region of the shield and 
for this period y a small hot- 
wire signal voltage is recorded. 
Precisely In the transition 
region from a high to low signal 
voltage, the absolute flow 
direction has the preselected 
probe adjustment angle of Oy. 


With approprlat triggering the test signals of several 
sequential blade passages — as shown in Fig. 3*1^> can be 
superln?)osed on the oscillograph screen or plotted on an x-y 
recorder as a statistical mean value (based on 1,500 blade 
passages) by means of the synchronized sampling technique. 

In these signal curves the position of o„ relative to the 
limit of the signal period t on the pressure bide Is determined 
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by y, l.e. by means of the tip :learance signal recorded at the 
same time — and the position relative to the position of the 
blade pressure side at the Impeller exhaust Is determined by Y. 

By gradually changing the probe adjustment angle and by 
moving the probe over the diffuser width In several radial test 
positions towards the Impeller output, the total direction field 
of the absolute flow In the region of the diffuser Intake can 
thus be determined. 
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Fig. 3*15. Test signal curve of the hot-wire shielded 
probe. About 50 Individual signals superimposed. 

Above - Internal reference signal. (Clearance transducer 
installed at a point shifted by t/2 against the direction 
of the probe.) y»40®. n»12,000 rpm, R/R2"1.075 
Oscillograms: horlz. mO ys/cm, vort. 1 V/cra. 

3.^.3 Probe Calibration and Frequency Characteristics 

For developing and testing the hot-wire probes a small 
calibration wind tunnel was setup with a nozzle outlet diameter 
of 20 ram. At the test point within the free jet of the nozzle 
the raach numbers can be set Independently of one another between 
M ■ 0.2 and M ■ 0.9, the densities can be Independently set at 
values between p > 0.8 and p ■ 2.0 kg/m3 and the total temperature 
between room temperature and 250®C. 
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For direction calibrati n the hot-wire shielded probe — 
as indicated schematically on the left in Fig. 3*16 — is 
rotated in the parallel stream of the Jet flow and the probe 
output signal is recorded on "~i x-y recorder. In so doing, 
during the transition from direct to shJ elded flow onto the hot- 
wire, this gives the change in d.c. voltage shown from which 
the zero setting for the probe ^ “ 0° is determined. 



Fig. 3.16 Calibration of the hot-wire shielded 
probe in a parallel Jet: change in d.c. signal 
voltage U/U as a function of the angle 0 at 
M - 0 . 7 , - 288 K, P - 1.5 kg/m3. 

Oscillograms: a.c. component of the hot-wire sign 
signal, horizontal 20 ms/cm, vertical 0.1 V/cm. 



The angulrr range in which the signal transition is 
accomplished is ■ ±4® in the exam; le shown. A subsequent 
test showed no measurable in the poslclon of the probe zero /50 

setting and in the width of the transition region in the Mach 
number range studied. 

The osclllcgrams shown in Fig. 3*16 show the a.c. component 
of the hot-wire signal in each case in the angular position 
Indicated. While at an angle of ♦ “ -30® the oncoming flow 
strikes the hot-wire directly, at t ■ -4® initial smaller signal 
variations occur which reach a maximum in the probe zero setting. 

In the angular positions ^ ■ +4® and ♦ ■ +30® the hot wire 
then lies con^letely in the highly turbulent wake region of the 
probe shield. 

I 

? 
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For an exact evaluation of the hot-wire measurements, 
especially for correlating the test signal with a blade spacing 
data is also required on the frequency curve of the measuring 
system. In obtaining this information in the case of the 
direction measuring technique described, in which only time 
components are determined, only the phase shift of the test 
signals is to be considered. The frequency curve of the 
anemometer is dependent on a large number of parameters, 
such as the geometry of the transmitter, the arrangement of 
the bridge and the band width of the anemometer as well as the 
flow condition and the velocity of the medium in which the 
measurements are made. Under these circumstances, direct 
dynamic calibrations are very costly and time consuming. 

Because of the relatively small Influence of the phase shift 
on the accuracy of the measurement, results, even in the 
present case, obtained with Indirect methods are sufficient for 
making an estimate. Simple, indirect test methods, in which 
the transmitter current is modulated sinusoidally or in a 
square-wave by- means of a frequency generator to simulate a 
fluctuation in flow, have been described by Nielsen and Rasmussen 
in [3-11]. 

Fig. 3.17 shows results determined according to this 
procedure for the hot-’-’lre shielded probe. It is conspicuous 
in the graph that the phase shift decreases very sharply with 
increasing flow velocity. If after a Fourier estimate for the 
steep signal transmission of the shielded probe one applies about 
a 4-fold blade frequency, then it can be estimated that in the 
centrifugal compressor operating range in question (10,000- 
18,000 rpm) excitation frequencies of about 10-25 kHz are to /51 

expected at oncoming flow velocities between 100 m/s and 400 m/s. 

Thus in the evaluation, phase shifts of less than 30° are to 
be taken into account. Under the assumption made above, this 
corresponds to about ' 1 % of the signal period. 

3 . 4.4 Discussion of the Measuring Technique for 

Nonsteady Direction Measurement 

The probe design shown in Fig. 3*13 was selected from 
a number of possible solutions as the most suitable for the 
studies on the centrifugal compressor test stand. First of all, 
it was attempted, by means of continuous recording methods, to 
avoid the main drawback of this technique, namely that about 
50 time consuming individual measurements are required to 
determine the direction field in the region of 1 blade spacing. 

Probes with I'reely stretched individual or crossed wire 
arrangements, which use the well-known cos (p direction effect 
of hot wires (cf. [3.9]), were not considered here because of the 
complex heat transfer conditions (cf. section 3.4.1) and the 
small amount of protection with respect to dust damage. Also 
eliminated was a probe principle which at first seemed promising, /52 
in which hot film transmitters were operated in separation from 
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Fig. 3* 17* Result of an indirect frequency cur\e 
study for estimating the signal phase shift using 
measurements made with the hot-wire shielded 
probe . 

Key: A) Phase angle 

B) Main test range on the centrifugal 
compressor 

C) Frequency 


one another on the sides of a 50'’ quartz glass wedge. This 
method was not used since the risk of different changes in 
drift over time of the two transmitters on the glass substrate 
could not be eliminated. 

Finally, another measuring technique, which likewise 
operated on an Intermittent principle and was developed by 
Giese and Schreck [3*12], was tried out in a practical test on 
the centrifugal compressor. In this probe the hot-wire was 
stretched over a distance of 0.1 mm parallel to the flow direction 
in front of a galss cylinder with a diameter of 0.6 mm. The 
directional characteristic of the probe shows up when the hot- 
wire lies in the sphere of influence of the stagnation point 
flow line of the cylinder, i.e. in a region of low flow velocity. 
In the strongly fluctuating flow at the centrifugal Impeller 
exhaust it v'as very difficult, however, to determine a distinct 
recording medium in the signal curve. Moreover, with this probe 
design, hot-v/lre service lives of only about 15 minutes could 
be achieved. 

In the extensive flow studies on the centrifugal compressor 
test stand — in spite of built-in air filters which separated 
out dust particles over 1 uni in diameter — the last-named point 
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cook on increasing importance. With respect to the hot-wire 
life spans, the shielded probe possesses decisive advantages. 


Because of the principle of the measuring technique the 
hot-wire is in the direct line of flow basically only during 
about 50 s? of the test time, and during test pauses the hot- 
wire can simply be protected in the wind shadow of the probe 
shield. Moreover, the use of tungsten wire is more advantageous 
because of its tensile strength which is 5 to 10 times greater 
than that of platinum wire or platinum-rhodium wire (cf. C3-13], 

Table II ).^ 

Finally, shortening the length of the wire to 0.6 mm 
contributed significantly to an increase in the service life of 
the wire, so that test series lasting several days at an /53 

extremely high load are possible with these probes. Another 
side effect of this measure is that the aim of obtaining ' 
measurements which are as close to point measurements as possible 
is closely approximated and that the transition range of the 
test signal is reduced by about +4® to the value shown in 
Fig. 3 . 16 . 

This signal transition range is set by the fact that, in 
spite of using a stereo microscope and micromanipulator in 
making the probes, it is not possible to setup the hot-wire 
exactly parallel at a distance of about O .7 mm from the leading 
edge of the shield. If, however, one assumes an approximately 
constant degree of uncertainty in the making of the probes 
(angle of Inclination of the hot-wire to the edge of the shield), 
then it turns out that the ±4° signal transition range of a 
0.6 mm long hot-wire at a length of 2 mm would correspond to an 
aperture range of about ±13® and thus the accuracy of the 
evaluation would be considerably limited. 

3»5 The Synchronized Sampling Technique — an Electronic 

Evaluation Method for Test Signals from a Periodic , 

Strongly Fluctuating Flow 

The rapidly fluctuating, periodic flow events within the 
sphere of influence of the impeller flow — as are observed, 
for example with fixed semiconductor pressure transducers or hot- 
wire probes with high time resolution — can usually be recorded 
photographically and evaluated only by means of cathode ray 
oscillographs. 

The image impression resulting from multiple overlapping 


1. A certain disadvantage, to be sure, is the fact that tungsten 
cannot be soldered. For making the hot-wire shielded probe it 
was necessayr to use an electric spot welding technique. 
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of the resulting periodic signal curves is strongly disturbed 
by fluctuations over time in the sequential events, caused, for 
example, by fluctuations in the relative flow of the Impeller. 

The second crucial disadvantage of this evaluation method 
consists in the fact that the analysis of test results recorded 
photographically is very painstaking, time consuming and 
inexact. 

Instead of this, especially with respect to the wealth of 
measurement data in 3~dimensional flow studies, it would be 
desireable to have a method which averaged the stochastically 
distributed signal fluctuations and at the same time slowed /5^ 

down the changes in the signal is such a way that it would be 





Pig. 3* 18. Application, in 
principle, of the synchronized 
sampling technique: 
a - periodic pressure curve, 
b - definition of the local 
instantaneous pressure value p’_ 
and the average pressure value p 
in the region of a blade spacing 
T, c - statistically averaged 
pressure curve p (t/T) and fluc- 
tuation factor f*(t/T) after SST 
separation of periodica lly re- 
curring and stochastically oc- 
curing signal components. 

0 . After passing through the i 


possible to convert the infor- 
mation directly into digital 
form and store the measured 
values for separate data pro- 
cessing. By using the rela- 
tively simple electronic 
principle of the synchronized 
sampling technique (SST) it 
becomes possible to fulfill these 
requirements. A detailed dis- 
cussion of this technique and a 
description of the SST instrument 
used is found in Appendix II. 

The accuracy of the signal 
averaging process basically 
depends on the constancy of the 
teat condition, in particular 
the rotational speed of the 
compressor, during the test period 
of about 2 minutes. The 
small deviations of ±0.03? in 
the desired adjusted value of the 
rotation speed creates extremely 
favorable conditions for this on 
the centrifugal compressor test 
stand. 

As an example, in principle, /55 
for possible SST applications. 

Pig. 3*18 shows schematically the 
changes in wall pressure over 
3 successive blade spaclngs - 
- whereby the real pressure curve 
p'(t) is composed of a periodically 
recurring average pressure pul- 
sation pr (t) and stochastically 
distributed pressure flucutations 
lampllng process the average non- 


steady pressure p (t) and the effective value of the fluctuation 
Intensity f'(t) are available as separate signals ©. Strictly 
periodic test signals can thus be separated from stochastically 
distributed fluctuations and at the same time they can be slowed 
down to the extent that recording them on x-y recorders or — 
after appropriate analog-digital conversion by means of a digital 
voltmeter — on perforated tapes or magnetic tapes presents no 
difficulties . 



Fig. 3 . 19 . Slowed down and statistically averaged 
transfer of a rapidly changing wall pressure curve 
(oscillogram: 0.25 ms/blade spacing t) to an x-y 
recorder (2 minutes/blade spacing tssT) using the 
synchronized sampling technique. 


For example. Fig. 3.19 shows the transfer of a wall pressure 
curve shown on the oscillogram used to an x-y recorder. In 
parallel with this, the slowed down, statistically averaged 
pressure curve was stored In about 250 digital values. The 
special importance of this is In the sampling signal analysis 
in regions with highly turbulent, strongly fluctuating flow, e.g. /5o 
In the viake region at the impeller output where quantitatively 
exact data can be obtained only in this way. 

3.6 Basic Procedure for Determining. Processing and 

Evaluating Nonsteady Measured Variable'? 

After discussing various test methods which have been used 
for investigating nonsteady flow events in centrifugal com- 
pressors and after describing the technique for averaging 
periodic, strongly fluctuating test signals by means of the 
synchronized sampling technique, the entire chain of measurement 
will now be briefly outlined extending from the electronic 
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determination of the flow variables via the processing of the 
test signal to the final evaluation. Detailed Information on 
the measurement technique and evaluation technique in the 
analysis of the total flow field behind the Impeller output are 
given below in Section 5 * 

Pig. 3.20 gives an overview of the various possibilities 
for processing the measurement data in the form of a signal 
flow diagram. 
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Pig. 3.20. Test signal flow diagram for the 
determination, processing and evaluation of 
nonsteady flow events occurring at the centri- 
fugal compressor. 


Key: 

A) Analog signal processing 

B) Centrifugal con^ressor with 
transducers for nonsteady 
flow variables 

C) l 4 -channel tape recorder 

D) Additional Electronics 

E) Cathode ray oscillograph 

P) SST signal averaging 

Q) Camera 

H) x-y recorder 


I) Digital data processing 

J) Operation 

K) On line, Off line 

L) Evaluation programs 

M) Central digital computer 

N) Digital voltmeter A-D 
conversion 

O ) Tape puncher 

P) Protocol printer 

Q) Diagram plotter 
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Directly at the test point in the centrifugal compressor /57 
the flow variables are converted via the transducer — total 
pressure and wall pressure sensor and hot-wire shielded probes — 
into electric analog signals. In addition, synchronously running 
trigger Impulses are also measured. The necessary additional 
electronics basically consist of a voltage supply, d.c. amplifiers 
and the hot-wire anemometer. After passing through this first 
stage, by interposing a l4-channel tape recorder (upper limit 
frequency 80 kHz), the possibility exists of storing the 
analog signals and separating them from the actual compressor 
experiment for further processing. 

Depending on the case, 5 different possibilities present 
themselves for recording and evaluating the signal (cf. Fig. 3*20): 

Direct photographic recording of several superim- 
posed signal curves from the screen of a cathode ray 
oscillograph with an instant-picture camera. This method 
is used to get a quick idea of the data. It is unsuitable 
for precise evaluation of a large number of quickly fading 
measurements . 

(g) With the aid of the synchronized sampling technique 
(SST), periodic, strongly fluctuating test signals are 
separated by averaging their fluctuating components and 
slowed down to the extent that the curvj of the average 
test signal and the fluctuation Intensity can be recorded 
directly on x-y recorders. Operation of the SST Instrument 
requires the saw tooth output of the oscillograph triggered 
synchronously with the test signal. 

From these dispalys are obtained results of the non- 
steady direction measurement, the curve of the pressure 
fluctuation factor and information with respect to the 
temporal correlation of all of the test results behind the 
Impeller output. 

0 Occurring in parallel with the SST scanning is the 
analog-digital conversion of the test signals by means of 
a digital voltmeter and the output of the some 250 values 
per period on a perforated tape puncher. In this way, all 
of the alternating pressure curves are stored on a central 
computer for further exact, automatic processing. /58 

0 Input data to be processed by the extensive evaluation 
programs on the central digital computer TR^40 Include, 
besides the pressure curve perforated tapes, the following: 
the results of the nonsteady direction measurement, the 
correct mean pressure values and temperature values over 
time, the operational data for the compressor €ind the 
necessary geometric data. 


The following were calculated: 

- the exact nonsteady absolute pressure distributions, by 
superimposing the alternating pressure curves with the 
correct mean pressure values over time; 

- the complete flow fields (absolute system and relative 
system) in 4 coaxial test areas immediately behind the 
Impeller output and 

- characteristic Integral mean values for describing both 
the total flow field and the individual jet or wake com- 
ponents of the flow. 

! The results of the flow field calculation are stored 

‘ in magnetic tape files and printed out when needed. 

I In an easy-to-read form, the data for describing 

the complex flow field at the output of the centrifugal 
impeller can be given in three-dimensional representations 
on a diagram plotter. In so doing, the accuracy of the 
digital evaluation is coupled with the clarity of the 
analog signal presentation. 

4. Design of the Test System, Operational Measuring Technique 
and Preliminary Experiments 

4.1 Description of the Test Compressor 

The compressor is designed for a rotation speed of 18,000 
rpm and a mass flow of 7*16 kg/s.l Important data in the design 
point is listed in Table 4.1. 

j 

^ In the design of the compressor stage, besides the demand /60 

j for greater unit output — large absorption capacity with a 

I high stage prssure ratio — points of measuring technique were 

also considered which, in keeping with the aim of the project, 
facilitate exact studies of the impeller flow. To be included 
here are, above all, the diffuser, which is unbladed and connected 
nected to the outlet side of the Impeller — this permits 
unfalsified analysis of the impeller outlet flow — and the large 
stage dimensions, which stand in a favorable relationship to 
the size of the probes. 

Pig. 4.1 shows the centrifugal compressor Impeller studied 
with 20 radially ending blades which extend far into the axial 
inlet region. The Impeller was milled from an aluminum maleable 


1. The compressor was designed during the planning of the 
’’centrifugal compressor” research project by the member companies 
of the Porschungsverelnlgung Verb rennungskraft mas chlnen e.V. 
[Research Association for Internal Combustion Engines]. 
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TABLE 4.1 DESIGN DATA FOR THE CENTRIFUGAL COMPRESSOR STAGE IN 
QUESTION (z«20) ACCORDING TO [4.1], RELATED TO THE STANDARD STATE 
IN THE COMPRESSOR INTAKE (p^o bar, T^q n»288.1 K). FOR 

tE measuring points SEE*FIG. 


INDICATION OF THE 


4.2 


Rotation speed d 

Reduced output factor l, 

Total pressure ratio « 

"t,05 

Total efficiency 

*'i,05 

**i,02 

Total temperature ratio n* /in 

■'■t2 *tO 

Mass flow '^,3 

Power required U 

18 000 
0,83 
2.683 
3,014 
0,80 
0,91 
1,406 
7,16 
844 

Velocity vector diagram 

263,9 

At Impeller Intake 

(a - external) ®mla 

112,0 

*la 

287,0 


0,852 


90,0 

•la “ 

23,0 

Velocity vector dlagreun ^ 3^^3 

377,0 

Impeller exhaust ^ 

335,1 

4 

0,894 

V 

312,9 

»/a 

120,0 

Wp m/s 

136,0 


21,0 

•2 " 

61,9 
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Fig. 4.1. Centrifugal compressor impeller, blade 
number z = 20, blade outlet angle 3s2 = 90°, 
material AlCuMg2P45 • 


alloy. Thus, in terms of the hardness of the metal, curcum- 
ferential velocities at the tip of approximately 525 m/s are 
permissible. 

The main dimensions of the flow channel are given in /6l 

Pig. 4,2. The centrifugal Impeller is taper bore mounted, so 
that the flow is sucked undisturbed through an axial intake 
pipe with a diameter of 28l mm. Because of the high stress of 
centrifugal force, all of the blade units of the impeller stand 
radially. The mean camber line (in cylinder section A-A) has 
an ellptical shape. The approximate curve of the flow cross- 
section A/Aintake Impeller blade channel over the average 

flow path lenght 1/lmax shown in the lower right corner of 
Fig. 4.2. Behind the Impeller outlet extends an unbladed cen- 
trifugal diffuser with a constant orifice area over the radius. 

In order to even out the diffuser flow field in the circum- 



Pig. i|.2. Centrifugal compressor flow channel. 

I^n dime^ions and measuring point positions 
(5) ... (5)* measurements in mm, (Intake width 
of the diffuser which is equal in area, bg* 26. 7 aim). 

Key: A) Throttle ring E) Throttle ring 

B) Unbladed diffuser F) Coi*recting chamber 

C) Centrifugal Impeller G) Area curve of an im- 

D) Mean camber line of peller blade channel 

the blade in cylinder 

section A-A 


ferentlal direction a throttle ring was Installed at the diffuser 
outlet with a radii ratio of R/Ro ■1.91. The manner of operation 
of this ring is explained in detail in section 4.4. 

After a deflection in the axial direction, the air then 
flows into annular collecting chamber of equal area. The com- 
pressor outlet pipes branch off from this tangentially. 

Fig. 4.2 also shows the positions of the different test 
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points or reference points v'lthln the congjressor stage; 


compressor Intake (also see Pig. 4.3, no. 7), 

Impeller Inlet edge. 

Impeller outlet edge, 
diffuser Inlet region (1.0 R/R 2 1*2) 

diffuser outlet (R/R 2 * 1 . 68 ? In performance characteris- 
tics measurements) 
conflpressor outlet. 

4.2 Total Layout of the Centrifugal Compressor Test Stand 

Pig. 4 . 3 . shows an overall view of the complete centrifugal 
coiqpressor test stand. 




Pig. 4.3> Overall view of the Centrifugal ompressor test stand. 

1 - D.C. swivel bearing motors (2 x max. 750 kW) 

2 - Torque scale 

3 - Planetary gearing 

4 - Centrifugal compressor 

5 - Air Intake line NVi 800 with dust filter 

6 - Venturi tube NW 500 

7 - Coiqpressor Intake working planes NW 500 and Nw 28l 

8 - Annular throttle slide NW 300 

9 - Pneumatic membrane valve with pump by-pass line 

10 - Hot air line NW 800 with cooler and noise danq^er at the 
outlet 
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The test compressor is dri^^ by two coupled, separately 
excited D.C. electric motors - (^ - which together have a 
maximum output of 1,500 kW. These can be continuously adjusted 
over a range of 0-3*000 rpm. The motors are moimted in pegdulum 
bearings, so that with the air of the calibrated scale - @ - 
the power take-off torque can be determined. The planetary 
gearing - jf|) - in a contact awangement (gear ratio, 1:9) and 
the centri^gal compressor - ® - are mounted on a separate 
base. This suppresses the propagation of motor oscillations and 
reduces the noise of the rapidly turning part of the system. 

In principle, it is possible to run the test stand with an open 
or closed air circuit. With open operation, the air is taken 
in from the atmosphere and again blown out into the atmosphere. 

In this mode of operation the performance limit of the drive 
motors is reached at a compressor rotation speed of about 
22,000 rpm. If one wants to increase the rotation speed up to 

the m;.xltnum value of 27,000 rpm, then the switch must be made /63 

to close operation, with which the compressor Intake pressure 
drops below atmospheric pressure and the air is recooied behind 
the compressor. 

All of the experiments described with' n the scope of this 
study at rotation speeds up to 18,000 rpm were performed with 
an open air circuit : 

- The air flows through the Intake line - - (Fig. ^.3) 

to the conqpressor after dust particles > 1 um in diameter have 
been separated out by a filter (prc'conditlon for hot-wire 
meacarements) . The nlr mass flow is determined v^h the air of /6'l 

a standard built-in short Venturi tube NW 500 - . In the 

straight pirt of the line in front of the Venturi tube inlet 
a series of straightening blades and vurbulence screens was 
built in. These were spaced about 6 meters apart. Their pur- 
pose is to eliminate the effects of the upstream pipe bend. 

In the hot-air line on the pre^urc side an annular motorized 
slide is built in as a throttle - (8). On the pressure side, 
the compressor can be qulcklv relieved by means of a pneumatically 
activated membrane valve - (P - in case the pump limit is 
unintentionally exceeded. The air circuit terminates with an 
air cooler and noise damper - (l^ 

The entire test stand was equ-,^ped with extensive control 
and monitoring systems which are briefly discussed below. A 
detailed description is found in [^>2]. 

The pu^ose of an axial thrust control is to keep the load 
on the axial Journal bearing of the compressor caused by the 
axial thrust of the impeller (max. about 20 kN) under approxi- 
mately 1.5 kN. The tatic pressure increase in the centrifugal 
Impeller is used as u reference quantity to control the pressure 
of the con«>ressed air which acts upon the axial thrust compen- 
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Fig, 4.4. Centrifugal co.npressor test ..tand with 1,500-kW 
drive motors. 


sating piston co-rotating on the compressor shaft. 

For an alternacive possibility of automatic pump limit 
mon5torlng, when a dynamic pressure limit value is exceeded in 
the compressor intake, the by-pass line mentioned above around 
the tnrottle unit (Fig. 4.3, position 8) is opened and the 
compressor is relieved on the pressure side. This assumes, of 
course, that the location of the pump limit is already known, 
i ccordlng to our experience, ooeratlonal instabilities, such a- 
incipient rotating stall, can ..e most accurately determined by 
installing a wall pressure sensor in the impeller Intake region. 

In the case of close c ’ "cuit operation, two additional 
fixed-point control loops are used to keep the Intake pressure 
and Intake temperature constant. 

With the operation monitoring systems we should, in the 
first place, mention the non-contact control of the clearance 
between the individual impeller blades and the compresser 
housing. This Is accomplished by means o.^ a capacitive measuring 
metho^ which was adapted to the special operating conditions 
encountered during use on highly-loaded compressor stages (cf. 
[4.3]). Thus on the centrifugal com.pressor — after appropriate 
calibration — the changes in clearances at the Impeller inlet 
and outlet — which are about 0.5 and 0.7 mm respectively hen 
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the impeller is not running — where monitored with an accuracy 
of better than 0.05 mm. At the same time, the output signal of 
the clearance transducer serves the purpose of exact coordination 
and evaluation of the nonsteady test signals which are picked 
up within the impeller or behind the impeller outlet. 

To monitor the vibration of the centrifugal compressor a 
commercial measuring system was used which operates on the 
induction principle with seismic transducers. These are mounted 
on the compressor housing in an axial, tangential and radial 
direction. With appropriately careful balancing of the rotor, /66 
the vibration velocities in all of the transducer directions 
were less than the very good value of 3 mm/s. 

Other monitoring included the continuous recording of 
compressor bearing temperatures, separate for cylinder roller 
bearings and axial journal bearings, the recording of to and 
from oil temperatures for the compressor and gears, as well as 
continuous, digital measurement of the flow of oil for the 
circulation of oil for the gears and separate for the supply 
of oil for the ball bearings and journal bearings. All together, 
some 30 different variables are monitored automatically, plotted 
on a recorder and when preselected limit values are exceeded 
this activates both audial and visual alarms. 

^.3 Conventional Operational Measuring Technique and Data 

Determination 


In order to evaluate the overall operational behavior of 
the centrifugal compressor stage, characteristic performance 
studies and a number of detailed measurements with conventional 
time-averaging measuring techniques were carried out prior to 
the nonsteady flow measurements in the diffuser inlet region. 
These measuring techniques, which are necessary for measuring 
performance characteristics and, in part, for setting the 
operating levels, determine the rotation speed, mass flow, static 
pressure, total pressure and total temperature. These techniques 
and the corresponding parts of the data acquisition system are 
discussed in the appendix. 

A trigger disc — attached on the gear side of the com- 
pressor shaft — with 60 radial notches uniformly distributed 
on the circumference of the disc is used to measure the speed 
of rotation. The disturbances in the magnetic field caused by 
these notches are converted into electric impulses by magnet- 
resistor transducers. These impulses are summed by a counter 
with a 1-second integration time and shown in digital form as 
a direct speed of rotation reading in rpm. 

As already mentioned above, the air mass flow is determined 
with the aid of a standard short Venturi tube installed in the 
compressor Intake line. The DIN 1952 specification for flow 
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measurement were observed in the layout for a mass flow range of 
2>10 kg/s. These specifications were also observed in the in- 
stallation of the Venturi tube and during the evaluation of the 
measured values. 

The following, time-averaging measuring techniques were used /6? 
outside the direct sphere of influence of the nonsteady Impeller 
flow. 

The usual wall pressure holes with a sharp-eged inlet 
(diameter 0.8 mm) were used for static pressure measurements. 

The axis of these holes is perpendicular to the housing contour 
of the flow channel. 




Pig. 4.5. Conventional, time -fiveraglng probes. 
Left: total pressure probes; right: NTC total 
tengjerature probe, measurement in mm. 


The probe shown on the left in Fig. 4.5 is used for the 
total pressure measurements. The probe shown on the right was 
developed to measure total temperature. The temperature sensor 
is a NTC thermistor. The reading given by both of these types 
of probes is Independent of direction with subsonic flow and 
in an angular range of t30^ , The use of glass-encapsulated 
NTC temperature sensors in place of thermoelements in the cen- 
trifugal compressor studies is based on a number of advantages: 

- these transducers have long service lives in the strongly 
fluctuating flow. They are insensitive to dust erosion 
and dust deposits; 






the high Internal resistance (1,000 Q) m^es possible 
higher signal voltages by a factor of 10^ during 
operation In a bridge circuit. With less expense 
(without an amplifier) these slgpial voltages csui be 
accurately transfered and further processed. 


On the other hand, a possible disadvantage may appear that, 
with NTC thermistors, the calibration curve le more tenperature- 
dependent in comparison with thermoelements. However, since 
with precise thermoelement measurements a calibration curve for 
the sesnsors must also be teiken, this point does not have much 
bearing. 



In developing the temperature 
probe special attention was 
placed on the fact that the 
ten?>erature recovery factors Rf 
should be easily reproducible 
in operation, l.e. effects due 
to the positioning of the probes 
(e.g. due to increase carrying- 
off of heat in the vicinity of 
the wall) should be largely 
eliminated. The following 
equation applies: 


Rf * 
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with T » static temperature, 

■ true total temperature and 

Tt.a * Indicated total temper- 

_ ^ ® 


ature. 


Because of the necessary 
sensitivity of response, the 
Pig. 4.6. Curve of the tempera- probe has ventilation holes . 


ture recovery factor Rf plotted through which the test gas flow 


over the mach number M of the out axially along the shaft of 

oncoming flow for NTC total the sensor. At the same time, 

temperature probes (see Pig. the front portion of the NTC 
4.3) • sensor (about 6 mrnn long) is 

free-standing. Both of these 
measures reduce the carrying-off 
of heat between the thermistor bead and the shaft of the probe. 
This considerably Increases the Rf value and at the same time 
ensures that the temperature reading remains unaffected by the 
Installation location of the probe. Por this type of probe, 
fairly optimal conditions lie at an area ratio (ventilation 
opening/probe inlet opening) of 0.4. 


Pig. 4.6 shows the corresponding curve of the recovery 






factor Rf as a function of the mach niunber M of the oncoming 
flow. The ^If-tone curve path corresponds to the variation 
width of the Rf curve which depends on the accuracy of manu- 
facture and the effect of contamination after logner periods 
of operation. 

For fast determination of the mean conqpressor Inlet 
temperature at operational level settings, a platinum resistance 
thermometer with an active sensor length of 0.35 m Is also 
mounted In the Intake line NW 500 (see Pig. 4.3» position 7) 
obliquely to the direction of flow. The probe possesses a 
cylindrical dust shield. An Rf correction of the Indicated 
value Is unecessary, however, because of the low mach number of 
the oncoming flow, (M < 0.15). 

With the aid of a data acquis t Ion system, the flow 
variables measured with the conventional measuring systems are 
converted Into electric potential signals. They are scanned 
In a measuring cycle of about 1 minute and stored on perforated 
tapes. The system possesses 20 NTC temperature terminals, 10 
pressure measurement points with ~ In the adapted test range — 
Individual pressure transducers (e.g. for the small differential 
pressures encountered during mass flow determinations) and 2 
pressure measuring point switches ^canl valves) each with 48 
wall pressure and total pressure terminals. The pressure trans- 
ducers measure differential pressures with resryct to atmo- 
spheric pressure. The exact atmospheric pressure Is determined 
by means of a Hg bulb barometer. 

In addition, by means of the data acqulstlon system. It Is 
also possible to record the exact position of up to 20 probes 
simultaneously by scanning the path potentiometer and angle 
potentiometer on the probe adjustment devices. 

4.4 Preliminary Experiments (Performance Characteristics 

Measurements) 


Before the start of the actual test program designed to 
analyse the flow field behind the Impeller outlet, the general 
test conditions were studied for the use of the nonsteady 
measuring technlqires and at the same time the operational 
characteristics cf the entire centrifugal In^ressor stage were 
obtained In the form of a performance graph. This was done with 
the aid of various conventional measuring techniques. 

The measured pressure and ten^erature values revealed that 
the nonsteady measuring techniques described In Section 3 would 
be loaded within tb<>ir operating limits up to 18,000 rpm. 

In the blade! diffuser, the circumferential distributions 
of total pressures, total temperatures and wall pressures showed 
large dysymmetrles . It was found that, starting from the com- 
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Him, , 






pressor outlet (see Fig. ^1.2, position 5) and working back 
through the collection chanber and the unbladed diffuser, 
disturbances of the axially symmetric flow field could be 
detected clear into the compressor inlet. Thus, for all that, 
nonsteady measurements in a certain radii ratio at any desired 
position on the circumference would no longer have been rep- 
resentative for this radii ratio and anyhow the already large 
expenditure for measuring the 3-dlmensional developemnt of the 
flow after leaving the impeller v.'ould have assumed impractical 
dime»' slons . 


To even out the flow a throttle ring was installed a short 
distance in front of the diffuser outlet (see Pig. 4.2). This 
was designed in such a way that in its narrowest cross-section 
a critical state appears and thus downstream dlstrubances cannot 
react on the bladeless diffuser. The effects of these measures 
are explained with the aid of the following figures: 

First of all, using as an example the wall pressure 
distribution in the bladeless diffuser at an operating level 
of 16,000 rpm. Pig. 4.7 shows the dlstrubed flow field before 

Installation of the throttle 

ring. The drop in pressure in 
the region of the outlet pipe 
can clearly be seen (270° 

(270° < 4> < 360°). 


With the variables shown 
in Pig. 4.7, a relative pressure 
variation gp can be defined as 
follows to characterize the 
intensity of the disturbance: 



^^max ” (^« 2 ) 


Fig. 4.8 comprares the 

,, ■ , measured Sp distributions with 

* ** ** ^ and without the built-in 

^ throttle ring within the range 

Pig. 4 . 7 . Distribution of the of performance characteristics 
static pressure ratio over the of the compressor. While without 
circumference of the unbladed the throttle ring relative 
centrifugal diffuser, pressure variations up to a 

n “ 16,000 rpm, maximum of 15 % were observed, the 

m ■ 5.7 kg/s. distribution on the circumference 

is evened out under the influence 
of the throttle ring to such an 
extent (e < 2 . 5 %) that satisfactory conditions exist for making 
the nonsteady measurements. 
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Pig. 4.8. Distribution of the relative variation 
In wall pressure Cp in the compressor performance 
graph, relative to^R/R 2 = 1.1. 


To be sure, as shown In Fig. 4.8, the throttle ring reduces 
the width of the performance graph available for measurements 
by 10^ (l8,000 rpm) up to 40$ (10,000 rpm) of the original mass 
flow difference between the stopper limit and the pump limit . 

While the system stop limit is shifted to a new "throttle ring 
curve" (cf. also Fig. 4.9), the pump limit also shifts for 
the throttle ring configuration by about 25 - 35 % of the original /72 
width of the performance graph towards smaller mass flows. ^ 

For the practical situation, it follows from this that by 


1. Extensive studies on the development of operational Insta- 
bilities (cf. [■^.10], p. 71 f*) have shown that the disturbances 
In the pressure distribution along the circumference trigger 
premature pumping even with mass flows which still lie in the 
stable range of operation with an axially symmetrical flow field. 
As a further difference, it was found that, without the throttle 
ring, continuously increasing, low-frequency pumping developed 
(fgump ^100 Hz), whereas with the throttle ring a rotating stall 
suddenly developed (rotational frequency about 60 % of the blade 
frequency). In both cases the operation? '' instability in the 
region of the diffuser inlet are initiated before the total 
flow upstream and downstream is measured from this. 
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homogenizing the flow field in the circumferential direction — 
in the present case, for example, by means of an adjustable 
ring throttle at the diffuser outlet — the operational range 
of the compressor is substantially Increased. 

For determining the per- 

- ' formance graph for the compressor. 

Pig. 4.9, the following variables 
were measured: 

- corrected rpm 



Pig. 4 . 9 . Performance graph of 
the centrifugal compressor sta 
stage studied, with built-in 
throttle ring. 

Test points at which addi- 
tional nonsteady detail 
measurements were made in the 
region of the diffuser inlet. 
Key: A) Rotating stall limit 
B) Throttle ring curve 


n//5J [rpm]. 


with 


corrected mass flow 
m 4 q ticg/sj , 

total pressure ratio 

®t,04 * Pt4^PtO » 

internal total efficiency 

(4.3) 


/U 


” ^tO^* 

Each of these variables was 
determined at 4 operation levels 
for a total of 5 rotation speeds 
between xi//Sq » 10,000 rpm and 
18,000 rpm. All of the variables 
were related to the standard 
state in the compressor Inlet: 


- total pressure 

»to.N * 

- total 'temperature 

Prom this we get the correction 
factors for 


/?4 


and 


- pressure: 

- temperature 
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In addition to the wall pressure distribution along the 
entire flow channel, we also measured total pressure profiles 
and total temperature profiles in each of the 9 equidistant 
points of the intake line radius in the connpressor inlet (test 
point (oy, cf. Fig. 4.2 and Fig. 4.3, position 7) or over the 
diffuse width in the radii ratio R/R 2 ® I.687 (Pig* 4.2, test 
point (p ) . 

In the profile measurements in the diffuser, first of all 
the distributions of total pressure were measured at 6 uniformly 
distributed positions on the circumference. In so doing, along 
with the probes shown in Fig. 4.5, three 3-hole probes were 
also used. Besides the total pressure, these also show the 
flow direction. This direction measurement is necessary for 
exact averaging. At the same time, it is used for allnglng the 
largely angle-insensitive total pressure and total temperature 
probes with which measurements were made with them in the same 
positions in a separate experiment. By this means, it was 
ensured that the adverse effect on the flow field due to 
obstruction by the probes remains within limits, and such 
measurements of total pressure, direction and total temperature 
which were taken at the same test point were expressed in 
relation to one another. This requires that the experimental 
conditions with respect to the corrected rpm and corrected 
mass flow can be adjusted so as to be exactly reproducible. A 
control program for a table top computer was worked out with 
which the speed of rotation and the mass flow are continuously 
adjusted to the altered Intake conditions. 

In evaluating the performance graphs, the suggestions of 
Traupel (see [3*1]» P* 174 1 . and [4.4]) were taken into 
consideration for the correct averaging. 

By Integrating the flow field, measured in detail, in the /75 
compressor Intake pipe, the mass flow measurements were checked 
in a few points with the aid of a Venturi tube and very good 
agreement was found (deviation < l.^%) . 

To check the measuring and averaging method for total tem- 
peratures — which is especially Important for the evaluation 
of the nonsteady measurements and for exact determination of the 
Internal efficiency — the total increase in enthalpy (l.e. total 
temperature Increase), corresponding to the input of energy into 
the flow medium, was also determined Indirectly. We measured the 
torque and rpm of the motors, as well as the exact flow of oil 
and Increase in oil temperature in the compressor and gears. 

From this we determined the output of the motors and the oil 
dissipation losses. Mathematically , the exteranal disc friction 
loss was taken into consideration or the basis of data from 
Daily and Nece [4.5]. Ignoring the losses by radiation 
on the compressor and gears, the corresponding Increase in 
the total temperature in the test compressor was determined 



in a performance rating of these variables. The deviations in 
the measured tenperature increases were less than 1J{ in com- 
parison with the values determined from the performance rating — 
less than 1.5? in the case where the speed of rotation was 
18,000 rpm. 

A detailed discussion of the results for this complex set 
of factors is found in [4.2], 


5» Experimental Procedure and Evaluation to Determine the Non- /76 
steady ^low Behind the ingelleF 

5,1 Arrangement of Measuring Points and Test Program 


Fig. 5.1* Sectional view of the 
centrifugal compressor with the 
principle built-in locations of 
the various measuring systems: 

(a) conventional total pressure 
and total temperature probes; 

(b) cooled wall pressure sensor 
and window positions for optimal 
velocity measurement; 

(c) probes for time-averaging 
and nonsteady total pressure 


Fig. 5>1 shows an overall 
view of the principle arrangement 
of test points within the cen- 
trifugal comprerjpor stage. In 
keeping with the aim of the pro- 
ject, the main focus of the 
experimental studies is directed 
on the diffuser inlet region, 
where, for e'^ample, probe - © - 
stands for a large number of 
probes for time-averaging and 
nonsteady total pressure 
measurement and for nonsteady 
direction measurement. In each 
case here, the shaft of the 
probe lies parallel to the axis 
of the rotor. The probes are 
handled from the front wall of 
the diffuser over the width of 
the diffuser with the aid of a 
device to slide them into 
position and they are rotated 
into the desired angular position. 

The raea^rlng systems - ^ - / 77 

and - (^ - for the nonsteady and 
time-average wall pressure 
measurement are also installed in 
the Impeller region, each of them 
perpendicular to the housing 
contour. In a few places in the 
impeller region, glass windows were 
Installed for making optical 
velocity measurements in place of 


measurement and nonsteady direction 
measurement; (d) time-averaging 
wall pressure measurement. 
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the cooled wall pressure adpaters. The measurements of the velo- 
city curve in the Inpeller-blade channels and of the blade load \ 

along the outer contour represent a considerable amount of added ’ 

Information for interpreting the flow field, measured in detail, 
behind the impeller outlet. i 

To explain the exact measuring points in the diffuser inlet 
region, we Introduce the rectangular coordinate system shown in 
Pig. 5 . 2 . 

The positive z direction 
lies parallel to the rotor axis, 
i.e. in the direction of the /78 
angular velocity vector w. 

Perpendicular to this extends the 
r,y measuring plane and reference 
plane in which the flow angle 
a of the absolute velocity vector 
c is measured. It is not 
possible to determine exactly 
the flow angle y in the z, r 
plane with existing measuring 
techniques, including the dual- 
focus laser method for optical 
velocity measurement. Thus, in 
evaluating the velocity field 
at the impeller outlet, the z 
components are not taken into 
consideration. 

The existing relationship 
of the diffuser width b to the 
channel spacing t (between the 
radii ratios R/R2 = I.017 and 
1.151, b/t varies between 0 . 4 l 
and 0.32) picbes to be especially 
favorable here. It prevents 
the development of pronounced 
cross-current components in the 
z direction, so that it appears 
permissible to Ignore these 
components (cf. section 5 * 6 ). 

B'lg. 5.3 shows the exact location of the coaxial measuring 
planes and reference planes in which the flow field behind the 
Impeller outlet was measured and evaluated. 

The measuring surfaces designated I-IV lie in the following 
radii ratios; R/R2 - 1 . 017 , 1 . 039 , I.O89 and I.151. 

For purposes of clarity, in Pig. 5.3 only the measuring 
surfaces closest to (I) and most remote from (IV) the impeller 
are shown. 
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Pig. 5.2. Centrifugal compres- 
sor coordinate system and 
position of the r,y working 
plane and reference plane with 
velocity vector diagram behind 
the impeller outlet. 











Pig* 5*3* Definition of coaxial 
measurement and references at the 


outlet of an li®ej.ler._hlade 
channel. I: r /r s 1,017, o_ 

IV: rJ/R 2 * 1.151, aj 


By definition, the measuring 
surfaces in the radii ratios 
Rs/R 2 ii® through the middle 
of the inlet opeinig of the 
total pressure probes. The 
total pressure probes (see 
Pig. 3*10 and 3*12) were all 
installed with as®25® (measuring 
surface I), or 30° (measuring 
surfaces II-IV) and were 
ha.ndled over the diffuser 
width in this angular setting. 

By means of the gradlations in 
the measuring surface radii 
ratios, it is possible to 
exactly measure the development 
of mixing processes and equali- 
zation processes in the dis- 
turbed flow behind the impel- 
ler outlet. In preliminary 
experiments it was found that 


up to measuring surface IV the periodic character of the flow — 
which is required for the use of the total measuringtechnique 
described — still exists, however, as soon as a radii ration of 
R/R2=1.2 is reached, for example, it was not longer possible to 
detect any clear variation in the variables in the circumferential 
direction because of the extensive equalization of the flow in con- 
junction with a greater aiuount of fluctuation. 

Pig. 5.^ shows an overall view of the various measuring points 
in the diffuser inlet region. In each case, the measuring point — 
probe inlet for total pressure probes, shield edge for hot-wire 
direction probes and pressure transducer membrane for wall hole 
for wall, pressure measurements — is indicated by a circle. 

As described above, the coaxial reference planes for the non- 
i steady flow analyses are determined by the position of the 

! total p^ssure probes Q) . For the nonsteady direction measure- 

I ment -©- the hot-wire shielded probes are installed in the cir- 

! cumferentlal position ♦=270° in a wall stop with finely graduated 

I measuring positions in the radial direction. Hence the flow direction 

f can be detected in nearly identical radii ratios. Because of the 

i tight conditions of installation, it was soemtlmes necessary to 

I make ^^romlses in arranging the wall pressure measuring systems - 

that the valid vlaue in each respective reference /80 
plane is found by Interpolating measured values which are determined 
in the immediately adjacent measuring points. 

The test program for determining the nonsteady flow field behind 
the Impeller outlet altogether include 11 operational levels at ro- 
tation speeds between 10,000 and 18,000 rpm. These are marked separ- 
ately in the compressor performance graph. Pig. M.9. The corres- 
ponding quiescent mach numbers of the absolute velocity at the Im- /8l 
outlet lie between approximately ^-O.^l-l.O. For the 10,000 and 
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1^,000 rpni rotation speed characteristic curves the mass flow is 
varied in the 3 operating points from the stop limit up t<'' the vi- 
cinity of the pump limit. By this, it is possible to study, within 
wide limits. Influences if the Mach number and the blade channel 
load on the flow fields in the diffuser inlet. 



Fig. 5.^. Overall view of the measuring points 
in the diffuser inlet region with the measuring 
surfaces and reference surfaces for the non- 
steady flow analysis in the follov/lng h/R. radii 
ratios: I - 1.017, II - 1.039, III • 1.085, 

IV ■ 1.151. Circumferential spacing: 360®/z ■ 18®. 

Measuring systems: 1 « dynamic total pressure; 

2 « time-averaged total pressure; 3 ■ dynamic 
pressure, front wall; ^ » time-averaged pressure, 
front wall; 5 « dynamic pressure, rear wall; 

6 “ time-averaged pressure, rear wall; 7 ■ non- 
steady flow direction; 8 • ln5)eller blade clearance. 


In order to detennlne the flow immediately after the Impeller 
outlet with as much detail as possible, the probes in measuring 
surface I are moved in 8 steps (step width z/b ■ 0.1) between 
test point z/b ■ C.l on the front wall of the diffuser and z/b ■ 

0.9 on the rear wall of the diffuser. In the remaining test 
surfaces II-IV the flow process is determine.’, in the wide positions 








f 

I 

t 




z/b » 0.1, 0 . 3 , 0 . 5 , 0 . 5 , 0.7 and 0.9. 

! 

The entire test program consists of *1 Individual measurement \ 

series which were performed one after the ither: 

i 

a) performance graph measurement, to determine the stage 
characteristics and the mean total temperatures behind 
the Impeller outlet; 

b) determination of the correct, mean pressure value over 
time in the diffuser inlet region with the aid of special 
pneumatic or hydraulic pressure measuring systems; 

c) measurement of the dynamic wall pressure and total 
pressure curves with miniature semiconductor transducers: 
and 

d) determination of the direction distribution of the 
absolute velocity vector in the region of the outlet flow 
of a blade channel. 

This procedv e presupposes the setting of operating points 
which can be accurately reproduced. In so doing, besides the 
reduced rotation speed and reduced mass flow, which are continu- 
contlnuousl) checked during the experiments and adjusted to the 
desired values, in addition all the other variables are uniformly 
related to the standard state in the compressor inlet (total 
intake pressure pto m » 1.0133 bar, total intake temperature 
TtO,N ■ 288.1 K). * 

In order to obtain Information on the flow field behind the /82 
inpeller outlet with as much general validity as possible, the 
nonsteady measuring systems are triggered with the blade frequency. 

By including the sampling signal processing, one thus practically 
obtains one nonsteady flow curve, statistically averaged over 
all of the Impeller blade channels. This characterizes the 
outlet flow behind any desired blade channel of the centrifugal 
Impeller. 

5.2 Determination of Constant Variables over Time 

Mean values for total temperature, as well as total 
pressure and wall pressure, are necessary for the exact mathe- 
matical determination of the nonsteady flow field behind the 
litqpeller outlet. 

The mass averaged total temperature is determined from 
the measurements to define the con^ressor performance graph in 
the radii ratio R/R 2 ■ 1.687, because — as explained in 
Apr sndix III — total temperature probes which are Installed 
immediately behind the impeller outlet in the region of very 
strongly disturbed flow show considerable deviations from the 

I 

I 

I 
1 
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mass average total temperature. 


The mean total pressure values, correct with respect to time 
and measured with the oil-filled capillary probes, are given by 
the summation of the excess pressure reading of the Indlvlc'ual 
pressure transducers connected to each probe and calibrated with 
the barometric reading. 

To compensate for the remaining, small circumferential 
dysymmetrles ~ in spite of the throttle ring — a total pressure 
Isobar distribution in the range 1.0 < R/R 2 < 1.2 is plotted for 
each measuring width pos?.tion z/b, and from this a total pressuie 
profile over the diffuser width is determined by arithrietical 
circumferential averaging for the measuring surface radii ratli.e 
I-IV. 


In a similar way, the wall pressures measured with the 
pneumatic systems are averaged over the circumference and the 
reference values on the front wall and rear wall in the 
measuring surfacos I-IV are interpolated from the curve of th? 
Increase in wall pressure in the radial direction. 

In each case the correction of the pressure reading I’or 
the pneumatic system with respect to nonsteady effects was 
carried out according to the method described in [3.5, p. 67 fj. 

It turns out, however, that because of the relatively small 
wall pressure fluctuations behind the Impeller outlet, the 
reading for the pneumatic systems deviates less than 0.15? from 
the correct mean wall pressure value with respect to time. 

5.3 Determination of Nonsteady Variables 

For accurate reconstruction of the disturbed flow field 
in the diffuser inlet region the various nonsteady data on wall 
pressures a^^d total pressures and on the distritution of 
pressure fluctuations and flow direction must first of all be 
classified with respect to one another and to the impeller over 
time. Furthermore, restricting the eva.luatlon of the data and 
flow analysis to the characteristic outlet flow of a blade 
channel requires chat the location and shape of the blade wake 
be determined relative to the impeller. This Implies, as 
indicated, in principle, in Fig. 5.3t the determination of the 
wake shift ^ in the individual measuring surfaces I-IV and 
over the diffuser width. 

The mutual cooridnation of the ^est signals with reoect 
to time is assured by uniform triggering with the blade frequency. 

The determination oi the blad> wakes over the diffuser 
width is explained, in principle, by Fig. 5>5« The data used 
here are the seunpllng signal plots of the nonsteady wall 
pressure curves and total pressure curves, and the corx'e spending 
fluctuation factors. 



starting from the known 
trigger position (start of the 
sampling process), first of all 
the position with respect to 
time of the outlet edge of the 
blade is determined by simple 
geometric coordination of 
trigger points and test points. 
For example, in the wall 
pressure curves (z/b = 0 and 
z/b = 1) the quantity a cor- 
responds exactly to tg/2, 
because in this case the blade 
clearance transducer used as 
the trigger installed at the 
impeller outlet is shifted 
exactly 1/2 of a blade spacing 
with respect to the spacing 
network for the measuring points 
(see Fig. 5.^, position 8). 

With the total presure curves, 
it logically holds that 
b = ts/2- g because — as lig 
Fig. t).3 snows — the blade 
passes the probe inlet opening 


earlier by a time component of 


/8H 


Fig. 5.5 Sampling signal 
recordings of the curves for 
wall pressure p, total pressure 
Pt and the corresponding fluctu- 
ation Intensity f’ for deter- 
mining the blade wakes over the 
diffuser width, operating 
point n//^ « 10,000 rpm, 
ih/ 0^/60 ® 3-81 kg/s (different 
ordinate scale for wall pres- 
sures and total pressures). 


In the immediate vicinity 
of the impeller — Pig. 5-5 
shows the ratios in the radii 
relationship R/R 2 = I.OI 7 — 
the blade wakes are visible as 
pronounced total pressure peaks. 
Assuming a dip in relative 
velocity and taking as a basis 
the secondary flow curve in this 
region which is explained later 
on, ‘it follows from the 
velocity vector diagram — 

Pig. 5.2 — that there is an 


increase in absolute velocity c on the order of c u, and thus 
with a Comparatively small variation in static pressure there is 


an increase in total pressure. This interpretation is supported 


by the curve for total pressure fluctuation f which, after tne 


relative calm curve in the region of the pressure side of the 


/85 


1. In addition, deviating from the simplified graph shown here, 
the small twist of the blade over the diffuser width and with 
respect to the trigger position is taken into consideration in 
the evaluation. 
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blade, In each case here, clearly shows an Initial Increase In 
the fluctuation Intensity. 

In the wall pressure curves the sharp drop In pressure 
seen In the Impeller region during the transition from the 
pressure side to the suction side Is also still visible at a 
certain distance behind the Impeller, so that here the blade 
wake is set In the middle of the pressure side. The time 
component c thus determined corresponds to the blade wake shift 
6 n In the measuring surface in question. This is shown in 
Pig. 5.3. 

Up to a radii ratio of 
R/R? » 1.039 the dips in the 
blade v/akes can be detected 
in the form of local total 
pressure peaks. For larger 
radii ratios, in which the wake 
dips are obviously largely 
mixed, it is more difficult to 
exactly determine the wake 
shift. A certain approxi- 
mation is first of all obtained 
by extrapolating the blade wake 
in the form of a. logarithmic 
spiral over the known reference 
points In R/R 2 = 1 . 0 , 1.017 
and .039. From the qualitative 
fluctuation curves in R/R 2 = 

1.017 and 1.039 one sees that 
directly behind the impeller 
outlet the pressure fluctuation 
intensity in the region of the 
suction side clearly exceeds 
the values on the pressure side. 
Starting from the notion that 
the flow is mixed with increasing 
distance from the impeller out- 
Fig. 5.6. Example for deter- let and that by means of 

mining the blade wakes in the turbulent exchange the increased 

diffuser inlet region, same flow fluctuation also extends 

operational data as in F1 5.5, into the marginal zones of the 
b' / b since 03 1 Cs II-IV flow region on the pressure side, 

(cf. Fig. 5.6).* * the blade wake separating line /86 

is placed in the region of 
increasing pressure fluctuation 

intensity. By the simultaneous evaluation of several measurements 
over the diffuser width in a measuring surface, the data on the 
shape and path of the blade wake in the diffuser inlet region take 
on additional certainty. As an example of an evaluation performed 
in this way. Fig. 8.1 shows the three-dimensional path of the 
interface between the outlet flow of two adjacent blade channels. 
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After the position of the blade length in the sampling j 

analog recording of the wall pressure and total pressure curves I 

is determined, these points are also determined in the digitalized 
signal curves. The intermediate curve path, which corresponds to 
a channel spacing of cs, is described by approximately 250 
values of the pressure transducer output voltage. 

In an "IDA” computer program lor nonsteady data evaluation /8? 

(cf. [4.2]) the Integral mean voltage value is determined over 

tg, set equal to the correct mean pressure value over time in 

the corresponding measurement point. Based on this, with the i 

aid of the known pressure transducer calibration factors, the ’ 

exact nonsteady absolute pressure flow is reconstructed. Then 

the pressure curves are standardized by mesms of a SPLINE 

interpolation program to exactly 101 values per channel division 

(Interval in the circumferential direction Ay/t = 0.01) and 

the^.o are made dimensionless with the corresponding mean pressure 

value directly in front of the impeller inlet. Assuming a 

linear static pressure curve over the diffuser width, finally 

the static pressure curve is determined from the wall pressures 

for all z/b positions, in which directly measured total pressure 

curves exist, by interpolating between corresponding y/t 

positions on the front and rear wall of the diffuser. 

In the same positions over the width of the diffuser 
measurements are also made to determine the absolute flow 
direction. The measuring technique and the evaluation of the 
data have already been discussed in detail in Section 3.4.2 
with reference to Figure 3* 15. For each probe position z/b, 
the exact position of 5-8 pre-selected probe setting angles within 
a wake channel spacing is determined. These measurements are 
evaluated in the form of an isocline graph (cf., for example. 

Pig. 6.5). Prom this then is obtained' the direction input 
data for the IDA program by interpolation for each z/b position 
at 11 equidistant reference points (&y/t » 0.1). 

With the aid of the SPLINE interpolation program already 
mentioned, intermediate values in an interval of Ay/t * 0.1 are 
calculated, so that in each of the 9 (or 5) x 100 points of 
the z,y screen (in measuring surfaces I, or II-IV) data is 
available on the local total pressure pt (y,z), the static 
pressure p (y,z) and the absolute flow direction a (y,z) to 
describe the flow field behind the outlet of an Impeller 
blade channel. 

In addition, as another nonsteady variable, the distribution 
of the relative fluctuation factor 



in the region of a blade spacing can be calculated from the 
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sampling analog recordings of the total pressure fluctuation 
Intensity f (cf. Fig. 3.18) and from the average nonsteady 
total pressure curve. 

5.M System of Equations for Calculating further Nonsteady 
Plow Variables 


In the "IDA” evaluation program a number of further flow 
variables was calculated from the distributions of the directly 
n^asured nonsteady flow variables p^ (y,z), p(y,z) and 
a(y,z). The exact determination of these variables for the 
outlet flow region of a blade cheinnel was described in the 
preceding section. By calculating these further variables, 
it is possible to give a detailed description and make a 
detailed study of the disturbed flow field behind the in^)eller 
outlet in the absolute system and in the relative system. 

Additional input data include the following: 
with respect to the standard state in the compressor: 

- total intake pressure p^Q jj = 1.0133 bar, 

- total intake temperature » 288.1 K, 

relative operational data: 


- corrected rotation speed 

- corrected mass flow from Venturi tube measurements 
*VENT* 


lass averaged total temperature at the diffuser outlet 
t4 * 


geometric data, such as 


- impeller outside diameter D 2 

- diffuser width at the impeller outlet b 2 » 

- tne measuring surface radii ratio R/Rg. 

The following variables are also included for calculating the 
thermodynamic state in which — as usual — the test air is 
regarded as an ideal gas: 

- the gas constant R, 


1. The position subscript is based on the data in Pig. 4.2; 
the subscript t means "total state." 



- the specific heat at constant pressure cp, which is used 
as a mean value for the range of teii^jerature increase 
under consideration and 

- the isentropic exponent k, which in each case fulfills 
the following condition: 

«c * (5.1) 

The following calculations were made: the circumferential 
velocity at the impeller outlet: 


'*2 ■ ^2 * * ’ ”red» ^5.2) 


with which all velocities were made dimensionless in the 
presentation of the results; and the circumferential velocity 
at the measuring point in question in the diffuser inlet region 

Uj . Uj • R/Rj. (5_3) 


The following equations are valid in each measuring 
surface behind the impeller outlet for any scanning element 
within the spacing of one flow channel. Prom the distributions 
of the static pressure and of the total pressure the local 
quiescent mach number of the absolute velocity is determined: 

(y,z) {1 - [P 3 (y,z)/p^j(y,z)l (5.^) 

_In the following iteration loop, which is run through 5 
tih'^js in the evaluation program, the quiescent velocity of 
sound is firgt of all determined from the mass average total 
temperature T^n — in the following passes determined from 
the local value Tt3» The connection equations serve to 
determine the absolute velocity c^, its circumferential com- 
ponent c „3 and — at the end of the ltL.rative calculation — 
the total^temperature Tt 3 
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1 

s 


’tn 


s [k • R • Tfcij(ytS) 3 ^^^; » (ic • R • 

1 1 

(5.5) 

1 

CjCy.z) s af.3(y,z) • » 

(5.6) 

Cu3(y.z)= Cj(y,z) • cos «j(y,z) ; 

(5.7) 

T^3(y,z)* jj ♦ U3 • c^3(y,z)/Cp . 

1 

(5.8) 


In this process, Eq. (5* 8 ) corresponds to the well knows 
Euler moment equation for non-twisting flow to the impeller, 
in which the component of the rise in total pressure caused by 
Internal disc friction and turbulent exchange movements (back 
flow) at the impeller inlet and outlet are ignored. This 
problem is discussed in detail in secti >i 5-6 within the con- 
text of error estimation. 

The following equations for calculating the radial com- 
ponent of the absolute velocity C 1 . 3 , the relative velocity 
W3, its circumferential component Wu3 and the relative flow 
angle 83 result form the velocity vector diagram (cf. Fig. 5-2) 


Cy3(y**) * Cj(y,«) • sin a^(y,s) • 

^ 1/2 
WjCy.s) » [cpjCy.s) ♦ wJjCy.s)] , 

ftxCy**) • co8[w^j(y,t)/Wj(y,s)] . 


(5.9) 

( 5 . 10 ) 

(5.11) 

(5.12) 


Finally, the static temperature T 3 , speed of sound a 3 , 
relative mach number IV 3 , density P 3 , static enthalpy h 3 and 
entropy S 3 follow from the equations given below: 

T,(y,s) « T.,(y,s) - c|(y,*)/(2c ) , (5.13) 


62 


(5.14) 


* [k • R • ^ 

H„3(y**)* «3(7*z)/a3(7»s) , 


P3(y»z) * P3(y»s)/CR • T3(y,a)] , (5.16) 


h3(y»a) * Cp • T3(y,*) , (5.17) 

•3(7*2) a Cp • ln[Tj(y,z)/T^O,N^ * ^ * lnrP 3 (y* 2 )/P^Q^H^ » (5.18) 

whereby in the last equation s is arbitrarily set equax to zero 
for the inlet standard state; PtO,N> "^tCjN* 

In part using Indirect methods, thus we have developed 
all of the variables for a detailed analysis of inpeller outlet 
flow in the absolute and relative system. For an integral 
check of the accuracy of the measurements and of tho evaluation 
a mass balance Atda is setup over the Impeller outlet cross- 
section A^; 

11 

®IDA * ^3 I [ f »(5.19) 

y/t*0 z/bsO 

whereby because of the diffuser, which is equal in area, the 
following equation applies; 

A3 * Ag « Dg • w • tog I [ I d(*/b)]d(y/t). (5.20) /2i 

y/t*0 «/b*0 

As the error in this mass flow determination, the relative 
deviation with respect to the Venturi tube measurement is de- 
termined. (cf. Section 5.6); 

* ^®H>A " * (5.21) 


In order to ensure the coherence of the experimental data 
the measured flow angle corresponding to the mass flow difference 
are corrected and the evaluation is once again repeated between 
equation (5.5) and (5.19) with the flow angle 013 icorr* 
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®3*korr^^**^ * arc[ain aj(y,a)/(l ♦ 


( 5 . 22 ) 


With this correction, all of the measurement errors are 
unilaterally charged to the nonsteady direction measurement. 
Since this is in fact the moat uncertain, and since fj„ is 
generally smaller than 0.05, this simplified procedure seems 
allowable here. 

5.5 System of Ec^uations for Calculating the Characteristic 
Mean Value 


The variables already defined above and the distribution 
calculated in the program for instantaneous data evaluation 
(IDA) form the input data for an extensive integration sub- 
program (INTIDA) (cf. [4.2]). In the context of this study, 
the calculation equations for inqportant, characteristic mean 
values to determine the condition of flow at the impeller input 
and output are listed here. 


When using the following geometric data 
Dn 


°la»^ll 


- Intake line diameter directly in 
front of the Impeller input and 

- outside and inside impeller inlet 
diameters respectively 


the following equations are valid at the impeller input for 
-static temperature 




- speed of sound 


- density 


• (« R 


*0 * 


- inlet cross-section 

^0 


* « (Dq - » 


(5.23) 


(5.24) 


(5.25) 


( 5 . 26 ) 


- average absolute velocity (with non-twisting on- 
coming flow) 


/93 


s 




(5.27) 


t 


and especially at the outer section for 
- the circumferential velocity 


“la ’ "red 


- or the relative velocity 


( 5 . 28 ) 


(5.29) 


If, in Eq. ( 5 . 28 ) one substitutes the following diameter, 
which halves the intake cross-section, in place of 


* KDq ♦ 


(5. 30) 


then from this follows u,, thus from Eq. (*^ 29) w, , and for th‘ 
average, relative angle of the oncoming flow (measured against 
the circumferential direction) : 


s arc tg(cQ/Uj) 


(5.31) 


As mean values for characterizing the impeller inlet state, 
the following values are determined: 

- the relative mach number in the outer section 


VI, a 


* "la'»0 


( 5 . 32 ) 


and in the diameter D, 


- the average relative mach number 




- or the incidence of the oncoming flow 


(5.33) 


®1»MEAN CAMBER " *1 


(5.34) 


In the various measurement radii ratios of the diffuser 
inlet region the relevant thermodynamic mean values of state 
for the nonsteady flow are again — as w‘th the performance 
graph evaluation — derived from the mean values of the extensiv. 
variables of state — static enthalpy h and entropy s -- (cf. 
[4.4]). 
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The following equations apply 


hj * ^3 ®r3 ^3 * 

*3 


( 5 . 35 ) 


®3 ®r3 ®3 *^^3^^"^ ’ 

*3 


(5.36) 


T, . h, / Cp . 


( 5 . 37 ) 


^3 • ®30 • »tO,ll ' <" • *3> ' 

^t3* ‘ ( *3 • *r3 • *t3 * ■**3>'"n)A 5 


( 5 . 38 ) 


( 5 . 39 ) 


(5.^0) 


Thus for the mean total tenqperature ratio we finally obtain the 
following: 


"t.30 • *30<^t3 ' 


(5. 'll) 


Furthermore, mean velocity values can be determined from 
the experimental data. The mean values to be inserted into the 
continuity equation carry the subscript k and those to be inserted 
in the energy equation carry the subscript e. Mean values which 
derive from the Impulse flow are ol special Importance for the 
formation of the mean velocity vector diagram (cf. Fig. 5.2). 

They do not carry any special subscript. 


The following equations hold: 

*p3k* ^ I ^3 * ®p3 ’ • ^3> • 

*3 

*r3 • U **3' • «?3 • *^3^'"IDA • 


(5.42) 


( 5 . 43 ) 
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^u3 * ^ ( ^3 * ®r3 * ®u3 * ^*3 ^^“ida ^ » (5.^4) 


\j " < *»3 • «r3 • M • ^*3^ '"IDA 


(5.45) 


and for the mach numbers In the absolute system ana relative 
system the following equations hold using Eqs. (5.4), (5.5). 
(5.14) and (5.15): 


^*c3*^ { ^3 * ®r3 ‘ "t,c3“*^3^'“lDA 

*3 

• < / O 3 • «r 3 • "w 3 • *^^^'"ida 


‘3""1DA 


(5.46) 


( 5 . 47 ) 


i *<3**S5 • 


(5.48) 


*3 * *r3 * ^ • 

^ • M ‘»3 • ®r3 • *3 * *S^'"IDA * 
*3 

*3# • M *^3 • ®r3 • "3 • ‘*^3^'"IDA • 


(5.49) 


(5.50) 


(5.51) 


and for the mean flow direction In the absolute system and 
relative system the following equations follow: 


« arc 8in(e^j/o^) » 


(5.52) 


> arc ain(e^^/w^) * 


(5.53) 


In particular, from the mean velocity values the following 
equations are formed: 


1. Since the Individual measuring surfaces lie on constant radii 
ratios, the averaging over the angular momentum R * Cu Is 
simplified In accordance with Eq. (5.44). 


A. 





- the flow through coefficient for indicating the 
operational point 


*3 * «r31c'“2 


(5.5^) 


- the mach number diffusion ratio for characterising 
the retardation of flow in the impeller 


MR . 


(3.5b) 


- or the mach number diffusion ratio for the Jet 
region in the flow pipe along the outer contour 
of the impeller, which is especially threatened 
by flow separation 

“*3 ’ • ( 5 . 56 ) 

- the most important form factors for dt-.-^cribing the 
flov? distribution at tbe impeller outlet (cf. 

[ 4 .^ 3 ) 


c 



• «««'') 

t„(c) . 

»„(«) . 


* l*c»i/c* 
r3k r3 

* (5.57) 



(5.58) 



(5.59) 
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Finally, in order t'' set off the healthy region cf the Jet 
(subscript S) from the wake region (subscript T) vithin the 
disturbed flow distribution the following rei '.t'’ cnships v/ere 
determined; 


( 5 . 60 , 

and Eqs. (5.ly;, (5.20) and (5. 35 >-(5.5^) — maintaining this 
relationship — were recalculated separately fer the Jet region 
and the wake region. 

To characterise vhe typical Jet/weks< velocity distribution 
at the outlet of the centrifugal compressor I’np. ller end the 
maxing of the disturbance in the diffuser inxet region, the 
following equations w?'*e derived; 

- the relative area corrgjoneoc of the wake from Eq. (5.20) 









- the relative mass flow component of the wake from 

I Eq. (5.19) 

"T * ^A,T^“lDA » (5.62) 

- the relative velocity ratio from Eq. (5.^9) 

I 

I ^ 

▼ * ''3,T^*3,S (5.63) 


- and the interference parameter 


« « « (1 - y) • (5.6^t) 


In the efficienyc definitions listed at the end it should /98 
be noted that for the total ten^erature, or the static tem- 

t perature, the folic. -/ing equations were used, for reasons of 

coherence, in place of Eqs. (5.^0) and (5.37) which were »mlfor*»>'.? y 
used in the system of equations setup so far: 

i 

^t3 * * (5.65) 

^ . ( 5 . 66 ) 

Thus as also shown by an estimate in Section 5.6 
the effect of the internal disc friction and the turbulent 
exchange movements is largely taken into conslde^-'s.tlon in 
determining t'.e efficiency. This effect was ignored in the 

1 basic equation (5*8) used to calculate the rise in total 

I temperature. 

I 

The following equations apply for: 

I 

- the isentropic impeller efficiency 

n** « Ah" / Ah" « C’f\ 

• (5.67) 

• *oB«30 • »»0.« ' • >3 / (»; - V • 

? 

t ; 
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( 5 . 68 ) 


- the polytropic impeller efficiency 

P3 

«!"•(/ ▼ <»p> t • 

*• Po 

* (1 " 1 / •)• ln(I^ ' ^ ln(5j / 5 q) * 


- and the Isentroplc total efficiency 


(5.69) 


• ■ ’to,*’ 


5.6 Accuracy df the Selected Measuring Technique and /99 

EiTaluationiechnlque 

In ePtimatlng the error, one basically distinguishes the 
separate errors in the following areas: 

a) channel geometry, arrangement of measuring points and probe 
adj ustment ; 

b) calibration of the measuring system; 

c) experimental procedure ^d 

d) data reduction and te... evaluatl n. 

The adjustment accuracy of the probe adjusting devices ...s 
around ±0.2 mm for the adjustment along the long axis and ±C.2® 
for the angular adjustment. Thus p-' .sible errors In the 
adjustment of tho probes, just like variation in the geometry 
of the flow channels due to manufacturing and mounting 
inaccuracies, are on an order of magnitude which can be Ignored 
in comparison with other measurement uncertainties. 

From the discussion of the measuring program (cf. Section 
5 . 1 ) it becomes clear that the ability to exactly reproduce the 
operating level takes on special Importance with resp.. -^t zo the 
detailed reconstruction of the flow based on varloui individual 
measurements. The uncertainty in determining the rotation speed 
is less than ±0.03lp in determining the mass flow witn a Venturi 
tube installed in the standard manner it is less tnan ‘approximately 
±1!K. In this latter case, only the components due to accidental 
errors have any effect with respf ct to the reproducibility of the 
adjustment. 

In estimating the error* in the flow measurements behind the 
lonelier outlet, the usual statistica.!. ezror analysis for 
multiple measurements at one test point cannot be used directly 
because of the unsteady measurement signal curves and the signal 
processing with the aid of the synchronized sampling technique. 


Therefore, an integral measurement error is estimated for the 
individual measuring emthods. This is based on the results of 
repeated calibrations and laboratory tests, as well as the 
reproducibility of the experimental results. 

With this Integral error Information it is also attempted 
to take into consideration the uncertainties discussed in 
connection with Figs. 5.5 and 5*6 in the exact arrangement of 
the various nonsteady individual measurements uniformly dis- 
tributed over the variables. In so doing, the measurement 
uncerv-alnty increases for larger measurement radii ratios. 

In addition to this effect and the test errors caused by 
the temperature sensitivity of the semiconductor pressure 
transducers, the total pressure reading contains a certain 
error reading which arises from the secondary flow at the 
impeller outlet. The flow components in the axial z direction 
are measured over the total pressure probes, and in the 
evaluation they are added to the components in the r,y 
referenc'=“ plane. For the error estimate it is assumed that the 
C |2 velocity ocmponent (cf. Fig. 5*2) equals at most 20% of the 
absolute velocity c. 

Finally, we must still discuss Eq. (5.8) based on the 
Euler equation of momnets for calculating the distribution 
of total temperature behind the Impeller outlet. 

The general Integral energy equation for the mpeller 
runs as follows: 

= (hj ♦ Z^/2) - (hQ ♦ 5^/2) s + q, (5.70) 

wherein Lp stands for the specific work done on the fluid and 
cj" stands for the heat added. 

This equation is valid only for appropriately formed 
integral mean values, but not along an individual flow line. 

In general, for the individual flow particle, the condition for 
Eq. (5.70): dLp + dLQ = 0, with Lr = the drag work due to 
friction and Lq = the work of deformation, is not fulfilled. 
However, even for the mean value, the validity of this equation 
is not guaranteed under all circumstances, for the character 
of the flow can be so complex that no cross-sectional areas can 
be found within which no work is transferred due to frictional 
stresses (cf. C3«l]» p. 75 f ) . 

This restriction applies to a special extent to a boundary 
surface which is plv.ced around the Impeller circumference at 
the outlet of a highly loaded centrifugal impeller and which is 
mixed with flow fluctuations and back currents. 


/lOO 
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These exchange processes, and also the usual friction of 
the impeller, produce in their totality a torque with respect 
to the axis of rotation turning at an angular velocity of 
w. The energy equation for nontwisting flow striking the 
impeller runs as follows; 

s s 

St3 ~ ho ' ®p ~ * Vtt3 * **t'*'*‘ ^5.71) 

With respect to what has been said above, Eq. (5*3), 
which ic developed for an individual unit area at the impeller 
outlet, violates several of the limiting conditions mentioned, 
in using total temperatures, which have been calculated according 
to Eq. (5*8), static ten 5 )eratures derived from them according 
to Eq. ^5 .iM and their integral mean values (Eq. (5.37) and 
Eq. (5.^0)) care must oe taken that in calculating these values 
the effects of impeller friction or, behind the impeller outlet, 
the effects of friction on the diffuser walls are not ignored. 

The mean order of magnitude of this error is determined 
from the difference between the total temperature Tt 4 measured 
far downstream in the diffuser and the total temperature Tt 3 
calculated from Eq. (5.^0). In so doing, the maximum value 
for the range of error, again for the outer measuring radii 
ratio R/R 2 , equals I. 15 I. 

Tr.ble 5.1 lists the estimated, mean relative measurement 
errors for the individual variables. In addition, for a few 
of the flow parameters derived from these, a mean range of error 
was calculated using the Gaussian Law of Error Propagation 
(cf., for example, [5.1]). 

Consideration of the errors shows that inaccuracies in the 
total temperature determination due to the use of the Euler 
equation of moments, Eq. (5.8) have only a relative small effect 
on the velocity evaluation. At the same time, this equation 
certainly produces the best possible approximation of the real /102 
total tem^ erature distribution, as long as no possibility exists 
for direct measurement. Therefore, it seems justified to use 
Eq. ( 5 . 8 ' to determine the exact jet/wake velocity distribution 
at the impeller outlet. 

In comparison with the pressure measurement, the relatively 
large uncertainty in the direction measurement was taken into 
account in estimating the relative range of error for this 
measuring method. As a result, we obtain the large upper error 
limits for the Cj.^ wj velocity components. However, it 
remains to be considered that, in measuring direction, random 
errors prevail due to the form of the individual measurement 
and evaluation (in contrast to the continuous recording of 
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TABLE 5.1 LIST OP RELATIVE ERRORS FOR A PEW NONSTEADY FLOW 
VARIABLES AT THE IMPELLER OUTLET 
EXAMPLE: n = 1^1,000 rpm, A » 5.31 kg/s, R/R 2 = I.OI 7 . 



pressure measurements). In forming integral mean values of the /lO 
total flow, or of the Jet-wake flow regions, these errors largely 
cancel each other out. 

This also explains the good agreement of the mass balances 
from nonsteady measurements behind the impeller outlet and from 
Venturi tube measurements. For the relative variations in mass 
flows fjn determined according to tq. (5.21) we obtain values of 
2 % Immediately behind the impeller outlet. In the extreme case, 
for R/R 2 * 1 . 151 , these rose to about -7Jf. In this connection, 
the essential uncertainty consists in the consideration of the 
boundary layer portion on the dlfusser walls. 

A qualitative Impression of the accuracy of the selected 
measuring method o.nd evaluation method is also given by a 
comparison with velocity measurements performed in the meantime 
with a dual-focus laser velocimeter (L2F method) after the 
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experimental portion of these^J^vestlgatlons. For the same 
operating level. Pig. 6.13b (l^ , shows the velocity distribution 
Ciu/u 2 for an impeller blade cKannel ae» it was recorded with 
the L2P method directly at the outlet edge of the impeller. 

Pig. 6.8b shows the comparable velocity distribution Cj.o/u 2 in 
the radii ratio R/R2 = I.017. This was determined using the 
nonsteady probe measuring technique. 

If one considers the small differences in the measuring 
position in the formation of the jet-wake distribution and the 
considerable difficulties involved in the measuring techniques 

— above all in the region of the strongly fluctuating wake 
flow — then the agreement of these results, which were obtained 
with basically different measuring systems, seems very 
satisfactory. 

6. The Plow at the Outlet of the Centrifugal Compressor Impeller 

The measuring and evaluation techniques discussed in the 
previous chapter setup for the first time the possibility of 
analyzing more exactly the complicated flow field at the outlet 
of highly-loaded centrifugal compressor impellers. In the 
following discussion of the results of these studies, we will 
first of all restrict ourselves to measuring plane and reference 
plane I in the radii ratio R/R 2 “ 1.017, immediately behind the 
impeller outlet. In so doing, the results of direct measurements 

— nonsteady wall pressure distributions, total pressure dis- 
tributions and flow direction distributions — and calculated 
flow distributions, which were determined on the basis of cer- 
tain simplified assumptions, will be discussed separately from 
one another. 

This flow analysis is Interesting in two respects; 

1. On the one hand, the disturbed flow condition being 
studies determines the initial conditions for the 
mixing and equalization processes occurring in the 
bladeless diffuser behind the ingjeller. These will 
be discussed in detail — theoretically and experi- 
mentally — in section 7 and 8 below. 

2. On the other hand, these studies of nonsteady events 
in the transition region between the impeller and 
the diffuser permit imoortant conclusions to be 
drawn as to the flow mechanisms in the impeller and 
— in the context of a consideration of the losses 
involved — to the aerodynamic-thermodynamic 
efficiency of the individual stage components. 



6.1 Results of the Nonsteady Flow Measurements 


6.1.1 Wall Pressure, Total Pressure and Distribution 
of Total Pressure Fluctuation 


As a lualltatlve example of nonsteady pressure 
measurement with miniature semiconductor pressure transduce: , 

Fig. 6.1 shows wall pressure curves along the housing wall li. 
the impeller region an on the front wall and rear wall of the 


IMPELLER 
Housing wall 



DIFFUSER 

Front wall Rear wall 



OS|SS OS|SS 


/105 


Fig. 6.1 Nonsteady wall pressure curve over approximately 
two blade divisions, measured on the housing wall in the 
Impeller region and on the front and rear wall of the dif 
fuser. Operating level: n//0^ = 14,000 rpm, rti /Sq/^O “ 
5.31 kg/s; Oscillograms: horizontal, 0.1 ms/RE (RE = 
scanning unit); vertical, 0.09-0.14 bar/RE. 
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bladeless diffuser behind the ingjeller. The measuring points 
in the impeller region are specified by the relative outer 
contour of the in 5 >eller in the meridian section x/Sjn (cf. 

Pig. 4.2, impeller inlet: x/s^ « 0), whereas in the diffuser, 
the radii ratio R/R 2 is given. All of the oscillograms were 
recorded with the same sensitivity and in so doing the 
the exposure time was so chosen that 30 signal uns were 
photographically superimposed. A blade clearance transducer 
at the impeller outlet was used to trigger and at the same time 
order the test signal exactly with respect to time (cf. section 
5.3> (see Pigs. 5*4 and 6.2). The signal peaks for this 
transducer are visible on the lower edge of the image for the 
individual oscillograms. 
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Por Interpreting . these experimental results it seems 
appropriate to introduce the nonsteady energy equation Tor the 
absolute flow which reads as follows for compressible, 
frictionless flow: 


Dh^/dt « 1/p 9p/?t , 

i.e. the total change in total enthalpy h^ with respect to time 
is, under the conditions mentioned, directly proportional to the 
change with respect to time of the static pressure at a deter- 
mined point in space. 

In considering the results shown in Pig. 6.1, the following 
phenomena are conspicuous: 

Not only is the absolute flow periodically nonsteady with 
the blade frequency — as is to be expected with the transfer 
of energy to the flow medium c i the basis of the nonsteady 
energy equation — but also the relative velocity is subjected 
to irregular fluctuations which correspond to the width of the 
superimposed "pressure band." 

The following possible reasons for this phenomenon may be 
mentioned: 

- pulsations in the oncoming flow to the compressor; 

- inspiti of the throttle ring, residual circumferential 
asymmetries in the diffuser pressure field and 

- disturbances which can be superimposed in the flow, 
for example by fluctuations of separation regions in 
the impeller. 

A striking fact here is the increase in this signal 
fluctuation intensity behind the impeller outlet, above all in 
the region of the front wall of the diffuser where — as will 
clearly be shown by the results discussed below -- the disturbed 
flow discharged from the Impeller is especially pronounced. 
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This phenomenon in the static pressure curve illustrated the /107 ? 

mixing losses coupled with the equalization processes of the 

disturbed Jet-wake flow. ! 

The periodic wall pressure measurements observed in the 
in 5 >eller extend with increasing amplitude into the diffuser 
inlet region and are almost totally evened out only at R/R 2 s 1.2. 

For this region behind the Impeller outlet it can be deduced 
from Eq. (6.1) that the zone of energy fluctuation extends as 
far as the region of nonsteady flow. Con®)arable to this, a 
reversible exchange of energy between the two partial flows 
is introduced into the ':wo-dlmenslonal theory of Jet-wake 
mixing processes described in section 7. This flow mechanism 
is Jointly responsible for the rapid — in comparison with 
steady conditions — mixing of the disturbed flow behind 
centrifugal impellers. 

The nonsteady wall pressure curves shown in Pig. 6.1 were 
evaluated numerically in Pig, 6.2 for the flow region of an 
impeller blade channel. The isobar distribution of the static 
pressure ratio p/po shown in Pig. 6.2 is determined by inter- 
polating between the Individual measuring positions x/sm on the 
housing wall in the impeller region and on the front wall of 
the diffuser which connects to it. In so doing, the static 
intake pressure pg Immediately in front of the Impeller inlet 
was measured. In this graph, y corresponds to the circumferential 
radian measurement on the housing wall, x is the coordinate 
along the outer contour of the Impeller in the meridian section, 
pressure side and suction side of the blades are, as usual, 
labelled DS and SS respectively, and also the positions of the 
blade clearance trigger and of measuring planes and reference 
planes I-IV in the diffuser inlet are given (see Fig. 5.3). 

The position of the boundary lines of the flow region for a blade 
channel on the front wall of the diffuser (relative to the 
impeller) were determined according to the method described 
in section 5.3* 

Prom the isobar curve it can be deduced that the axial inlet 
region of the impeller up to x/s^j - 0,5 is, depending on the 
layout, relatively less loaded, and e region of greatest 
addition of energy to the flow medium is concentrated in the 
final third of the blade channel. This fact, described by 
Eq. (6.1), is readily illustrated in Fig. 6.1 by the magnitude 
of the pressure gradient dp/dt at x/s^ ® 0.735 and 0.868. /108 

immediately behind the Impeller outlet the sharp pressure 
gradient in the radial direction between R/R 2 " 1.0-1.05 is 
especially noticeable. This expresses another effect of the 
Jet-wake mixing processes in the diffuser inlet region which 
will be discussed in detail in Setion 8. 

Prom studies of the impeller flow field with the laser 
velocimeter it is known that the flow in the preceding operating 
level at x/s « 0.55 separates from the housing wal? and a 
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pronounced wake area Is developed In the region of the housing 
wall and the suction side of the blades up to the impeller out- 
let (cf. Figs. 6.13a and 6.13b). This wake area is indicated in 
the wall pressure curves of Pigs. 6.1 and 6.2 only by the 
tendency towards a smaller circumferential pressure gradient on 
suction aide, however a clear separation of separated and non- 
separated flow — — based slnply on the nonsteady wall pressure 
measurements — is not possible. Other experimental wall 
pressure measurements at other operational levels, which pro- 
i^^formstlon on blade channel load and the supply of energy 
to the flow medium, have been presented, for example, in 
reports [4.2] and [6.1]. 
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Typical examples of the nonsteady total pressure distribution 
Immediately behind the Impeller outlet are shown in Pig. 6.3 
such as they turn out after passing through the sampling process 






Pig. 6. 3* Total pressure distribution Pt3(y>z)/Pto 
the outlet of an impeller blade channel, measured in the 
radii ratio R/R2 “ I.017. 

Operating leyels; 

• « 18 006 rpm, A 

b n/»® • 14 000 irni, A 

c n//§J » 10 000 A 

and the three D plotter program. The perspective graph of 
these results is restricted to the region of the outlet flow 
of a blade channel in the radii ratio R/R2 “ I.017. In this 
representation, the latteral separating lines (y/t « 0 and 1) 
which in the actual flow -- corresponding to the dips in the 
blade wake — are twisted (see Pig. 8.1) were straightened out 
in order to simplify the drawing. What was plotted were the 
nonsteady total pressure ciu*ves Pt3(y»z)» related to the total 
Intake pressure in the compressor inlet Ptp» which was measured 
between the measuring points z/b = 0.1 on the front wall and 
z/b = 0.9 near the rear wall at 9 equidistant axial settings. 

A striking feature is, flrct of all, the basically siml'.ar 
curve — down to the samllest details — for all three operating 
levels shown. This fact reinforces the assumption that the 
flow mechanisms occurring in the Impeller are active in the 
entire range of performance characteristics studies, and as 
a result their relative importance with respect to the development 
of the flow field is largely maintained. To be sure, with in- 
creasing rpm the amplitudes in the pressure curve over channel 
division and channel width increase considerably, and at 18,000 
rpm the range of fluctuation in the circumferential direction on 
the front wall of the diffuser (z/b » 0.1) already equals 40Jt 
(1.3 bar!) of the mean total pressure (note the change in the 
ordinate scale in Pig. 6.3 a with respect to Pigs. 6.3 b and c). 

As a result, the total pressure curve measured in the 
absolute system gives hardly any indication of the strongly 
disturbed , let-wake velocity distribution observed in the relative 
system. The corresponding velocity vector diagram — Pig. 7.2 — /111 
explains this phenomenon which is characteristic for impeller:* 
with radially ending blades: 

In spite of considerable fluctuations in relative velocity, 
the absolute velocities and thus — because of the small 
fluctuation in static pressure (cf. Pig. 6.2) — also the total 
pressures in the two areas are on the same order of magnitude. 

To be sure, a considerable Increase ^n entropy occurs in 
the region of the separated wake flow, but at the saraj time the 
transfer of energy in this region of lower mass flow density 
is considerable greater than in the flow-through region. 

Conversely, the greater efficiency of tne compressor process in 
the Jet region is contrasted with a smaller reduced output factor 


•®/<0 * 7,40 kg/8, 
* 5,31 kg/8, 

/ 5 ®/«o . 3,81 kg/8. 



and thus a smaller addition of energy (cf. Fig. 6.6). 

Especially pronounced In these total pressure graphs of 
the area Immediately behind the Impeller outlet are the pressure 
peaks on the pressure side and suction side. They are caused 
by the wake dips of the Impeller blades and they Increase 
steadily over the channel width — - corresponding to the Increslng 
blade thickness. If one assumes that In this region of the 
blade wake flow only a small or, In the case of back flow, 
even a negative relative velocity prevails, then the total 
pressure levols can easily be derived from the velocity vector 
diagram. 

The signal processing of the nonstead^ total pressure 
curves by means of the synchronized sampling technique offers — 
In addj *:lon to the determination of the statistically averaged 
pressure curve — the possibility of measuring the fluctuation 
Intensity of the measurement signal f In an analog fashion on 
a second output of the sampling Instrument (cf. Figs. 3*18 and 
5.f>). 


As the result of such measurements made directly behind the 
In^'eller outlet. Pig. 6.4 shows t-he distributions of total 
pressure fluctuation coefficient f for the flow region of an 
impeller blade channel In 5 operating points of the compressor 
performance graph. In the equation which defines the local 
total pressure fluctuation coefficients: 

«».i) •'V(Pt * 

p^ stands for the Instantaneous value, while Indicates the 
statistically averaged total pressure. 

In conslderalng the Individual distributions, the first 
striking fact is that the fluctuation intensity — to be sure, 
on average — increases with increasing rpm (Pigs. 6.4 e, c, a) 
and — at consteuit rpm — with Increasing mass flow, however 
the basic shape of the dlstrl. uclon remains approximately the 
seune In €.11 cases. Characteristically, a more or less (steep) 
diagonally running surface separates a region of higher 
fluctuation intensity, the core of which Is shifted on the 
front wall of the channel of towards the pressure side, from 
a region of relative calm along the pressure side of the baldes 
and the re' wall of the channel. Here there Is already a 
hint of ^ ..e basic Jet-wake distribution of the disturbed flow. 
This can be seen In the total pressure fluctuation. In which 
the greatest flow fluctuations occur In the region of the wake. 
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6,1.2 Direction DlstrlbutJ.on of the Absolute Plow 

An exact separation and clearer delineation of the 
jet region and wake region In the area of the outlet flow of 
an Impeller blade channel Is possible by determining the 
absolute flow direction with the aid of the measuring technique 
described In Section 3*4 using sulelded hot-wire probes. 

Fig. 6.5 shows a selection of direction distributions of the 
absolute velocity vector measured Immediately behind the Impeller 
outlet for the region of one Inqpeller blade channel. In these 
graphs the flow direction a Is defined according to the graph 
shown In Pig. 5*2. 


I 
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As with the total pressure fluctuation graphs of Fig. 6.4, 
the characteristic basic form of the Isocline curve Is main- 
tained with small variations in the total performance range 
studied. Scoarated by a steep direction gradient, the jet region 
with laxge angles of flow a and the region of the separated wake /115 
flow with small emgles of fow and a correspondingly small mass 
flow component clearly stau i out. The wake region — nearly 
unchanged for all operating points — lies at y/t « 0.6-0, 8, 
with its core on the front wall of the channel. As a function 
of the speed of rotation, one observes in Pigs. 6.5 e, c, a 
an increase in the width (z/b <* 0.4-0. 7), while the variation in 
mass flow in question at 14,000 rpm (Pig. 6.5 b, c, d) reveals 
Itself only slightly by conqparison. 

Taken together, the results of the nonsteady measurements 
made directly after the in 5 >eller outlet already permit a very 
extensive reconstruction of the disturbed flow discharged from 
the Impeller. In this connectic-i, the nonsteady direction 
measurement deserves special attention. This makes possible 
the inqportant transition into the relative system for the 
theoretical treatment of tiK- jet-wake mixing processes. The 
signal processing by means of the synchronized sanpling technique 
also deserves attention. By means of the statistical averaging, 
this technique sets up the necessary requirement for a quanti- 
tative treatment of the strongly fluctuating flow. 

6.2 Examples for Calculated Distributions 

With the aid of the system of equations listed in Section 
5.4; the required flow variables, for a complete evaluation of 
the flow field at the Impeller outlet, are calculated from the 
measured nonsteady values. 

6.2.1 Conversion of Work and Total Temperature 


Using the Euler equation of moments, the angular 
momer am con?>onent of the absolute velocity Cu 3 and the total 
temp6i»*ture at the impeller outlet T^o follows fron the 
iteration loop of Eqs. (5.5)-(5.8). ^ 
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As an example of this evaluation, Fig. 6.6 shows the 
distribution of the reduced output factor u, 

M = Cu3'“3» (6.3) 


for the flow region of an impeller blade channel, with respect 
to measuring surface I in the radii ratio R/R 2 “ I.OI 7 . 


/II 6 



Pig. 6.6. Distribution of the reduced output 
factor M at the outlet of an impeller blade 
channel, R/R 2 * 1.017. 

Operating point: n//oQ = 14,000 rpm, 

* •%/^0 “5.31 kg/s 


It can clearly be seen that the conversion of work in the 
centrifugal impeller — which is expressed in the curve of the 
reduced output factor — is quite considerable affected by the 
Jet-wake distribution of the impeller flow. In the Jet region, 
the addition of Jet energjj is clearly less than the mean momentum 
value for the total flow \x « 0 . 8 CI, while in the wake region, 
with its small mass flow component, it even exceeds in places the 
values expected with an infinite blade number (p > 1). This 
phenomenon can be explained by the secondary flow field at the 
impeller outlet, where a strong Joint angular momnetum component 
occurs in the region of the front wall of the channel (cf. Fig. 6.9 
and 6.6). 

For the same operating point, Pig. 6.7 shows the corresponding 
total temperature distribution at the impeller outlet. In the 
Isotherm graph for the region of a blade channel a total tem- 
perature ratio TR was plotted: 

TR(y,«) * (T^j(y,z) - » (6.4) 
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TR is thus defined as the ratio of the local Increase In total 
tenq>erature to the mass averaged Increase In total temperature, 
whereby measuring point "4" (cf. Pig. 4.2) Is at R/R 2 * 1.687* /117 

outside the region of noncteady flow striking the probes. 



Fig. 6 , 7 , Isotherm graph of the total ten^jerature 
ratlc TR at the outlet of an ln^)eller blade channel, 

R/R2 '= 1 . 017 . 

Operating point: n//0T » 14,000 rpm, 

& “ 5.31 kg/s. 


For the case shown In Pig. 6.7* the mass average equals 
0.987 Instead of the exact value 1.0. In the section where 
errors were considered (Section 5.6) we discussed In detail the 
restricted validity of the Euler equation of moments, Eq. (5.8), 
for determining the local pressure distribution behind the 
Impeller outlet. In view of the still unsolvable difficulties 
Involved In making a direct, exact measurement of this quantity 
under the given flow conditions, the results shown In Fig. 6.7 
take on even greater Importance. 

As this typical total temperature distribution shows — for 
this disturbed flow at the Impeller outlet — the variable 
addition of energy to the jet region and wake region leads to 
total temperature fluctuations In the area of a blade division 
which are on the order of ±15-2058 of the mass averaged tem- 
perature Increase. 

6.2.2 The Velocity Field In the Absolute and Relative 

System 

With the measured distribution of total pressure, wall 
pressure and direction and the calculated total temperature curve, 
all of the variables are fixed which are required to calculate 



the con 5 )lete velocity field with tte aid of Eqs. (5.9)-(5.12) . /II 8 

For three operating points at li},000 rpm In the radii ratio 
R/R 2 * 1*017> Fig. 6.8 and 6.9 show the thus determined velocity 
distribution In the absolute and relative system. In each case 
for the flow region of an Inqpeller blade cheuinel. 

In the plotter graph of Fig. 6.8 the radial component cf 
the absolute velocity Cp 3 (y,z), with respect to the clrcumfentlal 
velocity of the liqpeller Up, was plotted for the entire flow 
chsumel, whereas Fig. 6.9 snows the curve of the related 
relative velocity vector W 3 /U 2 In magnitude and direction In 5 
circumferential positions y/t oVer the channel width z/b (cf. 

Pig. 5.2). 

In the graph of the radial velocity curve CJ. 3 /U 2 the 
steep velocity gradient clearly localizes the region of 
separated wake flow. On the front wall of the diffuser — 
starting from the suction side SS — this extends over approximately 
2/3 of the blade spacing and In so doing teikes up about 60 ^ 
of the channel width b. The main portion of the mass flow 
leaves the Impeller through the Jet region along the pressure 
side of the blade DS and along the rear wall of the channel. 

In considering Figs. 6.8 and 6 . 9 , the following points stand 
out : 


- the large degree of similarity of the disturbed flow at 
the Impeller outlet for various operating points. As 
already mentioned In connection with the total pressure 
results shown In Pig. 6.3» this implies that the relative 
importance of the individual flow mechanisms, which con- 
tribute to the shape and development of the region of 
separated flow In the impeller, is largely maintained in 
the operating region in question; 

- the relatively small difference In velocity, both In 
magnitude and direction, within the Jet region and wake 
region — this suggests a "block-shaped" model simpli- 
fication of the flow field for the theoretical treatment 
of the Jet-wake mixing processes in Section 7» 

- the small — In spite of the high fluctuation intensity — 
clearly detectable mass flow through the wake region, 
which plays an Important role In the theoretical analysis 
of the Intensity of Jet-wake mixing. 

Has a further Important result of the flow studies at the /121 
Impeller outlet, the distribution of the relative total pressure 
will be presented and discussed In detail In connection with the 
laser measurements In Fig. 6.1^ and on page . 101 . 


• Cr3ly.il/u2 Cr3ly.il/u2 
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Fig. 6.8, Caption on p. 90. 







Fig. 6.8. Plotter evaluation of the absolute velocity 
distribution Cj. 3 (y,z)/u 2 behind an impeller blade channel, 

R/R 2 = 1 . 017 . 

Operating points at n//e7 = 1^1,000 rpm, 

A /¥^/6qI (a) 6.07 kg/s; (b) 5 .31 kg/s; (c) ^.53 kg/s. 


Fig. 6.9. Vector graph of the relative velocity dis- 
tribution W 3 /U 2 at the outlet of an iit^jeller blade c 
channel, R/R 2 = I.OI 7 . 

Operating points at n//0^ = 14,000 rpm, 

A /5^/6q; (a) 6.07 kg/s; (b) 5.31 kg/s; (c) 4.35 kg/s. 


6.3 Analysis of Test Findings at the Impeller Outlet with 

Respect to the Jet-Wake Distribution 

Since it turned out that the qualitative study of the flow 
at the impeller outlet showed a large degree of similarity in the 
velocity distributions and pressure distributions for various 
operating points, in what follows, the detailed, quantitative 
analysis of the Jet-wake flow will be limited to 5 operating 
points of the compressor stage between 10,000 and 18,000 rpm. 

For overall evaluation of the flow condition in these performance 
graph points, the most important mean flow values at the 
Impeller inlet and outlet (R/R2 = 1.017) were listed in Table 
6.1. In this table, the equation numbers refer to the corres- 
ponding definition equations in Section 5.5. 

Figs. 6.10 and 6.11 below illustrate two particularly 
remarkable partial aspects of the test data listed in Table 6.1. 
Starting with the Isentroplc total efficiencies rij^ 03 “ ^t s 
(pos. 3 - R/R 2 = 1.017) and ni o4 (pos. 4 - R/R 2 =*17687) * 
which were formed Identically According to Eqs. (4.3) and (5.69), 
a loss ratio between the impeller and diffuser is setup: 

« (1 - * ( 6 . 5 ) 


Fig. 6.10 shows the curve as a function of the discharge 
coefficient f 3 . In this graph, the experimental points listed 
in Table 6.1 are indicated by the all-black symbols. A striking 
feature is the general decrease in values, i.e. the relative 
Increase in diffuser losses with respect to total losses, as the 
discharge coefficient Increases. Above all, as a function of 
rpm, at the transition from 10,000 to 14,000 rpm we first of all 
see a sharp decline, then at n 14,000 rpm the values vary 
within relative narrow limits between « O. 3 -O. 7 . 

In a first approximation, this curve suggests the conclusion /123 



TABLE 6.1 CHARACTERISTIC MEAN FLOW VALUES AT THE IMPELLER INLET /122 
AND IMPELLER OUTLET FOR 5, NONSTEADY OPERATING POINTS MEASURED 


n/ 

rpm 

10 000 

111 000 


18 000 


kg/s 

3,81 

6,07 

5,31 

4,53 

7,40 


Ql. 

Impeller inlet 


^l,a 

(5.32) 

0,464 

0,669 

0,653 

0,639 

0,861 


(5.33) 

0,363 

0,532 

0,512 

0,495 

0,682 

ABj 0 

(5.34) 

♦ 6,9 

♦ 2,9 

♦ 6,6 

♦10,3 

♦ 3,6 



Impeller outlet 

(R/R2 

= 1,017 

) 


(5.54) 

0,383 

0,395 

0,338 

0,284 

0,327 

"t.c3 

(5.W) 

0,552 

0,776 

0,750 

0,729 

0,926 

®3 ® 

(5.52) 

27,6 

29,3 

25,3 

21,5 

22,9 


(5.47) 

0,283 

0,417 

0,356 

0,303 

0,405 


(5.53) 

67,8 

73,2 . 

71,2 

67,7 

76,6 

MR 

(5.55) 

1,283 

1,276 

1,438 

1,634 

1,684 

MRg 

(5.56) 

1,497 

1,408 

1,641 

1,863 

1,844 


(5.57) 

0,126 

0,191 

0,183 

0,179 

0,164 

e^3<c) 

(5.58) 

0,009 

0,016 

0,013 

0,009 

0,009 


(5.59) 

0,036 

0,060 

0,061 

0,058 

0,057 

^30 

(5.38) 

1,226 

1,413 

1,462 

1,488 

1,779 

^t,30 

(5.41) 

1,529 

2,218 

2,219 

2,209 

3,449 

n*« 

(5.67) 

0,904 

0,926 

0,920 

0,881 

0,925 

n;- 

(5.68) 

0,907 

0,931 

0,925 

0,888 

0,932 


(5.69) 

0,944 

0,962 

0,954 

0,936 

0,949 
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Fig. 6.11. Relative mach num- 
ber diffusion ratio for the Jet 
region MRs as a function of the 
discharge coefficient J3, 
reference surface at Impeller 
outlet R/Rp ° 1.017 * 

see Table b.l. 


Fig. 6.10. Curve of the loss 
ratio Kt between the Impeller 
and diffuser as a function of 


the discharge coefficient <])3, 
(•AB * see Table 6.1). 


that — with the exception of 
operating points with an 
extremely small load on the 
stage — approximately 2.3 of 
the stage losses up to the dif- 
fuser outlet first occur after 
the Inqpeller outlet. The detailed loss analysis In Section 8 
will show, to be sure, that a portion of this Is due to mixing 
and equalization processes of the disturbed Impeller flow (and 
thus can reasonably be attributed to Impeller losses). 


Since the losses In unbladed diffusers with steady flow can 
be determined with satisfactory accuracy (cf., for example, [6.2], 
[ 6 . 3 ] and [6.4]), the loss distribution shown In Fig. 6.10 
underlines the necessity of studying flow mechanisms and loss /124 
mechanisms In the diffuser Input region In greater detail. 


Fig. 6.11 shows the curve of the relative mach number 
diffusion ratio determined according to Eq. (5*56) In the Jet 
retlon MRs as a function of the discharge coefficient 
Besides the almost linear relationship of MRs + f(?3) for the 
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individual rotation speeds, the relatively small effect of the 
rotation speed is, on the whole, striking; 

For operations points with a constant the Jet-wake 
configuration obviously establishes Itself in Just such a way 
that approximately the same retardation of flow is produced 
for the "healthy" Jet flow. A limiting curve of the possible 
flow retardation stand? out which 5.s not exceeded for the 
Impeller design in question. 
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Fig. 6.12. Area component of 
the wake, wake mass flow com- 
ponent p^, wake-jet relative 
velocity ratio and interference 
parameter y» plotted over the 
discharge coefficient 
Reference surfaces at tne impel- 
ler outlet R/R2 ■ 1.017 
, see Table 6.1). 


Detailed studies of the 
impeller flow which have been 
carried out in the meantime with 
the optical dual-focus laser 
veloclmeter (also see Section /125 
6.4), have well confirmed this 
data. Just like the data shown 
in Fig. 6.12 to characterize 
the Jet-wake distribution at 
the impeller outlet (R/R2 = 

1.017). In these graphs the 
following distributions were 
plotted, going from top to 
bottom: 

- the area component e of 
the wake, Eq. (5.61); 

- the wake mass flow 
component p^, Eq. (5.62); 

- the relative velocity 
ratio (wake-jet) y, Eq. (5.63) 
and 

- the interference parameter 
6, Eq. (5.64) 

as a function of the discharge 
coefficient Eq. (5.54), for 
three rotation speeds between 
10,000 and 18,000 rpm. Again 
in these graphs, the test point 
treated in detail in this study 
are Indicated by the blacked- 
In symbols (cf. Table 6.1). 

With the exception of the /126 
test points for the rotation 
speed of 10,000 rpm, where the 
velocity distribution at the 
Impeller outlet is clearly less 
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severely disturbed (cf. y and 6 curves), the characteristic mean 
values of the Jet-wake flow lie within these band widths: 

0,35 < e < 0,41, 

0,19 < li^ < 0,26, 

0,46 < T < 0,55, 

0,17 < « < 0,21. 


Thus these Interference variables are sometimes considerably 
smaller than the extreme values found in the literature, and 
one can assume that in the impeller outlet flow studied here — 
which by comparison is only slightly disturbed — correspondingly 
small losses will arise In the course of the mixing and 
equalization processes. For example, the flow analysis of the 
well known Boeing-Avlabs centrifugal compressor stage, at a 
total pressure ratio of “ 10, gave a wake area component of 
e = 0.7, an interference parameter of 6 * 0.5, and 9 points of 
the stage efficiency were taken Into account only for the 
mixing losses behind the Impt-ller outlet (cf. [6.5]). 

6.4 Discussion of Various Flow Mechanisms in the Centrifugal 
Impeller Which Contribute to the Formation of the Jet-Wake 
Distribution 


In the past, the study of conplicated flow processes within 
centrifugal Impellers was restricted to the range of low cir- 
cumferential velocities for the impeller, primarily because of 
the quite considerable difficulties in making measurements in a 
rotating system. Previous investigators, e.g. Fischer and Thoma 
[6.6], Hamrick and co-workers [6.?], Pujle [6.8], Senoo and co- 
workers [6.9] and Fowler [6.10] have managed to give a qualita- 
tively comprehensive description of fundamental flow forms. 

Only in recent times has the development of new, non-contact 
velocity measuring methods made it possible to make a quanti- 
tatively exact analysis of internal impeller flow processes, 
even at circumferential velocities and blade channel loads 
which are on the order of magnitude of modern industrial com- 
pressors. Thus, the impeller flow of the centrifugal compressor 
stage studied within the scope of this work was exactly measured 
and analyzed with the dual-focus laser veloclmeter developed 
by Schodl [3*2, 6.11]. The conqprehenslve investigation covered 
circumferential velocities of up to 400 m/s [6.12, 6.13, 6.14]. 

The discussion of the flow mechanisms in the centrifugal 
compressor impeller fundamental for the Jet-wake distribution 
is based on a few findings from these investigations. 


Oil 
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6.4,1 Optics L .-''low 'e asurements In a Centrifugal 
Compressor ~ :ipeller~ 

The following selection of characteristic Impeller 
velocity distribution refers to an operating point in the 
vicinity of the optimum line of the performance graph; 

n//5J • 000 A * 5*31 kg/s. 

The velocity measuring method using the dual-focus laser 
which we enqployed makes it possible, at any optically accessible 
measuring points, to determine the absolute and relative 
velocity vectors in a measiirlng plane which lies perpendicular 
to the axis of the laser beam. For showing the change in flow 
from the inlet to the outlet of a centrifugal Impeller, velocity 
measurements were made in 5 measuring planes. The numbers and 
locations of these planes are shown in Table 6.2. The laser 
measuring planes are perpendicular to the outer contour of 
the housing. The exact location of the measurement position is 
defined by the relative outer contour of the Impeller in the 
meridian section x/Sj„ (impeller inlet, outside; x/Sm * 0, 
ln«>eller outlet: x/Sjj^ ■ 1). 

TABLE 6.2 LOCATION OP MEASURING PLANES FOR OPTICAL VELOCITY 
MEASUREMENTS IN THE CENTRIFUGAL COMPRESSOR IMPELLER 
8_; LENGTH OP THE OUTER CONTOUR OF THE IMPELLER IN THE MERIDIAN 
” SECTION 



1 Measurement planes | 

No. 

x/s_ 

n 

7 

0,08 

10 

0,43 

11 

0,59 

13 

0,87 

14 

1,01 


In Figs. 6.13 a and b are plotted the distributions of the 
meridian somnonent of the absolute velocity Cjjj in each respective 
measuring plane. This meridian component is made dimensionless 
with the circumferential velocity at the In^peller outlet Ug* 
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wall 


Fig. 6.13a. The velocity distribution c_./u 2 in the 
inpellor blade channel, measured in the 3 optical 
measuring planes 7* 10 and 11 in the axial in^eller 
region (o ■ measuring points). 

Operating level: n//6n ■ 1^*000 rpm, 

A - 5.31 kg/s.° 






Pig. 6.13b. The velocity distribution Cm/u 2 in the 
impeller blade channel, measured in the two optical 
measuring planes 13 and 14 in the centrifugal impeller 


region (o ■ measuring points). 
Operating point: n//d^ *• 14,000 rpm, 
A - 5.31 kg/s. 


In the prospective views, the measuring planes were represented 
in a simplified way as trapezoids, and the dimensions of the 
channel (spacing and height) were shown in the same scale. 

"Front wall" indicates the limit of the measuring plane (z/b » 0) 
on the side of the housing which lies opposite the hub in the 
background of the figure (s/b ■ 1). The shaded areas along the /128 
pressure side DS (y/t « 0) and the suction side SS (y/t • 1) 
correspond in each case to 1/2 a blade thickness. 

The velocity distribution in measuring planes — 
immediately behind the Inqpeller inlet - and (Tj0 — at the end 
of the axial inlet region — show that the flow largely exhibits 
the expected distribution based on potential theory up to the 
beginning of the meridional bend in the channel. With a com- 
parltiveiy small load on the channel, we observe, a regular ve- 
loci; y errve with maximum values along the housing wall in 
the r . . dian section, i.e. in the circumferential direction on 
the suction side of the blades. 

A first disturtjance in the in^eller velocity dlsti*ibutlon 
is seer in measuring plane in the region of the front 
wall of the channel. In this area, the load on the blade channel 
increases sharply at the same time, as is verified by Fig. 6.2. 

This disturbance indicates an Incipient region of flow s^oaration 
which — as shown by the findings in measuring planes and 
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— becomes still considerably larger in the region of the 
front wall of the channel and the suction side of the blades 
up to the impeller outlet. 

Likewise for the outlet measuring plane , Fig. 6.l4 /129 

shows the loss distribution at the Impeller outlet in the form 
of an isobar graph p^ rel'^^tO N* determined from: 

( 6 . 6 ) /JM 

With p^ , = relative total pressure in measuring plane at the 

* Impeller outlet 

and p^f^ « = standard total intake pressure directly in front of 

* the impeller inlet. 


Rear wall 



Fig. 6.1i». Loss distribution at the outlet of the 
impeller blade channels: isobar graph of the rel^ive 
total pressure p+- rel(y* 2 ) in measuring plane , 
relative to the "Itandard" total intake pr^ sure p^o N 
directly in front of the impeller inlet. * 

Operating point: n//5^ = 14,000 rpm, 

A/oJ/«q * 5.31 kg/s. 

The relative total pressure in this calculation is defined 
analogously to the relative total enthalpy so that for loss- 
free flow pt rel “ ^tO N ® constant. A detailed discussion 
of the relative total pressure p^ ^ei with a derivation of 
Eq. (6.6) is found in Appendix ivl The entropy distribution 
used s(y,z) was determined according to Eq. (5.18) under the 
following assumptions: 
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- that the curve of nonsteady pressure values measured on 
the housing wall is constant over the measuring surface; 

- that the energy transferred to the relative flow due to 
friction on the housing wall can be ignored, so that the 
static temperature distribution can be determined from /131 
the energy equation for the relative flow with the 

aid of the optically measured velocities. 

Finally, as a basis for the discussion of flow mechanisms 
subject to friction in the centrifugal compressor imepller, the 
optically measured velocity curves in Pig. 6.15 are compared 
with the results of a calculation of the impeller flow based on 
potential theory. In the upper graph is plottc ; the relative 
velocity profiles w/u2 between the housing and the hub at one- 
half the blade spacing y/t = 0.5» and in the lower graph are 
plotted the corresponding curves between the pressure side and 
suction side at half the channel height z/b = 0.5. 



Fig. 6.15. Comparison of experimentally determined 
relative velocity profiles w/u2 with calculated curves 
of the Impeller flow based on potential theory. 
Operating point: n//0A * 14,000 rpm, 
m/ 5 T/ 0 Q =5.31 kg/s; 

yV“’: relative blade division between blade surfaces. 

Key: A) Measuring plane E) Theory 

B) Housing F) Pressure side 

C) Hub G) Suction side 

D) Experiment 
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The calculations based on potential theory were made with 
a quasi three-dimensional method. The calculation of the 
meridian flow is based on the flow line bending method given 
by Stanltz [6.15]. In this method, blade-congruent flow is 
assumed to 75% of the merldianal flow path in the impeller 
before a given flow-reducing deflection is taken into con- 
sideration continuously up to the Impeller outlet. The 
solution found in this way is used to determine the flow areas /132 
in which the velocities on the blade surfaces are determined 
with the approximation methods suggested by Stannltz and Prlan 
[6.16]. Other calculations which can be compared with these 
experimental data, among other things based on the "energy 
layering theorem,” have been presented by Traupel [6.27]. 

In summary, the following characteristic picture of the 
flow in the centrifugal compressor Impeller studied follows 
from the experimental findings available so far: 

1. The flow in the Impeller blade channel possesses 
form — into the radial portion — which can 
largely be described in terms of potential theory, 
this is shown by the good agreement between the 
loss-free calculation and the experimgntal results 
(Fig. 6.15). Up to measuring plane (l^ , the observed 
deviations can be explained by the boundary layer 
displacement which is not considered theoretically. 

The same is ture for the calculated velocity gra- /133 
dients in the circumferential direction. These are 
in good agreement with the experimental data up to 
the beginning of flow separation at x/Sj^^ * 0.6. 

2. In the region of the meridional bend in the channel 
observes on the housing wall a sharp retardation in 
flow. As a result of this and due to the rapidly 
increasing region of separation in the region of the 
front wall of the channel and the suction side, 

a displacement of the main flow occurs in the 
radial portion of the impeller In the direction of 
the pressure side of the blades and the rear wall 
of the channel. This effect appears in amplified 
form immediately in front of the Impeller inlet, 
so that the disturbed velocity curve measured in 
the diffuser inlet deviates strongly from the curve 
calculated in advance using potential theory eqatlons. 
This clearly Illustrates the large Influcence of 
the region of separation on the energy conversion 
in the Impeller and underlines the necessity for 
developing Improved calculating methods for this 
flow region. 

3. The position of the separation point in the meridian 
section changes only slightly, as shown by studies 



at rotation speeds between 10,000 and 18,000 rpm. 

In general, (cf. Fig. 6.13) the first sign of 
flow separation Is observed In the form of a 
disturbance In the velocity profile on the housing 
wall In measuring plane 11 . As this happens, the 
core of the Interference region (at constant rpm) 

Is shifted from y/t = 0.8 In the vicinity of the 
suction side up to y/t = 0.^ near the pressure 
side as the mass flow becomes Increasingly smaller 
[ 6 . 28 ]. 

4. After the beginning of flow separation, the wake 

region Increases very rapidly, ^etween measuring 
planes (x/sm = 0.8?) and at the Impeller 

outlet, Tts proportion of the available channel 
cross-section Increase from 20 to 35?» while the 
wake mass flow at the Impeller outlet In the test 
point In question Is approximately 18$. The jet- 
wake relative velocity ratio In the Impeller 
generally lies between y = 0.45-0.7 for the operating 
points studied. In none of the test points studied 
was It possible to detect a complete otagantlon of 
the wake flow by means of the velocity measurement 
method using the dual-focus laser. A striking /134 

feature Is the strongly pronounced and relatively 
stable velocity gradient between the main flow 

(jet) region and the wake region. While In the 
core 01 the area of separation a sharp rise In the 
absolute velocity fluctuation up to 1.5 to 3 times 
the values In the jet region are observed, the 
shape and position of the wake region are, by 
comparison, subjected to only small fluctuations 
[6.14], 

5. In the loss distribution pt rel/PtO N» 6.1'-:, 

Is found the same sharp delineatlon*between the 
wake of lower energy and the main flow. This 
suggests a small amount of mixing of the twr regions. 
Whereas In the wake region a pronounced l-^ss con- 
centration occurs In the corner formed by the front 
wall of the channel and the suction sl-^e of the 
blade. In the jet flow — In addition to small total 
pressure losses — Is found an exte. slve Isentroplc 
core. For this core. In the area along the rear 
wall of the channel and along the pressure side 

of the blades, even slight "overlsentroplc” values 
were measured (1.0 < Pt,rel/PtO N l^OI'. In- 
dependent of the more exact laser measurements 
(for Fig. 6.14 the error In local va3ues was esti- 
mated at a maximum of ±2$), this result had already 
been found using the nonsteady prooe measuring 
method [4.2, 6.1]. Comparable results are also 
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discussed by Prlan and Michel in [6.17] (also see 
[6.4, p. 102]). Traupel [6.29] has explained this 
phenomenon as follows In terms of the "energy 
layering mechanism’-^: 

"Energy layering has its origin in the very 
i sharp fluctuations in the flow velocity. If 

I in a channel large pressure gradients prevail 

I at the same time transverse to the flow 

i direction, then it follows from the law of 

; motion [sic] that particles with above average 

i velocity drift primarily into the region of 

high static pressure, and those with below 
average velocity drift into the region of 
low pressure. Prom this it is clear that 
velocity distributions are found at the outlet 
of a centrifugal compressor impeller which 
contradict the predictions of classical 
theories, and that the total pressures of the 
relative flow can in places even exceed the 
I value which is suggested by a theory of 

' frictionless flow." 

i 

r 

The characteristic loss distribution shown here for 
I centrifugal Impellers is essentially not new. 

Comparable results are discussed in the works of 
Hamrick and co-workers [ 6 . 7 ] and Domm and Hergt 
!' [ 6 . 18 ]. With the results shown in Pig. 6.14, 

however, it was possible — possibly for the first 
j time — to clearly detect an isentroplc flow core 

; in the Jet region for an angular velocity range /135 

such as is typical ^or modern, high performance 
centrifugal compressors. Thus, it was also possible 
! to confrim an Important aspect of the simplified 

flow model proposed by Dean [ 6 . 5 , 6.19] for calcu- 
I latlng the complex centrifugal compressor Impeller 

|| flow . 

I 

6.4.2 Qualitative Analysis of the Impeller Plow 

The flow separation revealed in the above test results at 
approximately 60% of the meridional channel length and the sub- 
sequent development of a pronounced wake area up to the Impeller 
outlet adversely affect the efficiency of the centrifugal com- 
pressor stage (see Pig. 8.25). The separation of flow decreases 
the slow down of the flow in the Impeller and thus the (also) 
desired conversion of kinetic energy into static pressure in the 
in^eller layout. The strongly disturbed Jet-wake velocity dis- 
tribution at the Impeller outlet also reduces the stage efficiency 


1. A detailed description is found in [3.1* third ed.J. 
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as a result of mixing and equalization losses, as well as due 
to loss-affected, nonsteady flow processes in the diffuser which 
comes after it. Accordingly, it is of prime Importance, using 
the results available, to study the formation and development 
of areas of flow separation in order to derive from them 
possibilities for Improved stage design. 

Prom currently available experimental results, which were 
obtained both from rotating individual channels (diffusers) and 
from centrifugal impellers, it follows that boundary layer 
separations in rotating flow channels are particularly promoted 
by two flow mechanisms [6.12, 6.13]; 

- secondary flow effects which lead to a correction of 
flow material of lower energy and 

- turbulence damping effects, i.e. a local reduction in 
turbulent shearing stresses under the effect of flow line 
bending and system rotation (Coriolis forces). 

By secondary flows is here meant all flows within the 
blade channels which possess a velocity component perpendicular /136 
to the blade-congruent main flow. Theoretically, we distinguish 
two causes of secondary flows (cf. Traupel [6.20]): 

- a primary vortex system which is Induced solely by the 
deflection of a (frictionless) flow which at the blade 
channel inlet possesses velocity gradients or relative 
total enthalpy gradients, 

- and a secondary vortex system which is limited to the 
(friction-affected) boundary layer region. The casue 
of this type of secondary flow is the existence of 
pressure gradients which are Impressed on the boundary 
layers by the main flow and which, because of the small 
velocities within the boundary layers, cannot be counter- 
balanced by any dynamic forces. The results is then a 
secondary flow from higher to lower pressure. 

For the case being studied her: of an impeller blade channel 
with radial ending blades, what results, therefore, under the 
effect of Coriolis acceleration and standard acceleration (due 
to flow line bending) is a collection of low-energy boundary 
layer materials in the region of the front wall of the channel 
and suction side of the blades (cf. measured secondary flow 
curve in Pig. 6 .I 6 ). The secondary flow effects occurring in 
centrifugal impellers have been theoretically treated in detail 
by Howard in [6.21], 

Besides the developemnt of the secondary flow field, the 
pressure gradients in a cross-sectional area of the relative 
flow in the ln?>eller induced by ideal field forces also affect 
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the turbulence structure of the wall boundary layers. On a 
concave curved wall or on the pressure side of the rotating 
blade channel, with transverse movements in the turbulent wall 
boundary layer, the faster flow particles are forced towards 
the wall under the effect of centrifugal or Coriolis forces, 
and the slower particles are forced away from the wall. Con- 
versely, on a convex curved wall or on the suction side of 
the blade channel the slow particles become stabilized in 
their path on the wall and thus reduce turbulent exchange 
movements with faster flow particles. This damping of the tur- /137 
bulent mixing also reduces the Reynolds shearing stresses and 
thus — with positive pressure gradients in the flow direction 
— the drag effect of the outer flow on the boundary layers 
close to the wall, so that, on the whole, the tendency of the 
flow to separate increases. In the case of flow separation, 
the turbulence damping effect also Involves the free plane of 
shear between the Jet region and the region of separated flow. 

As a result, the turbulent mixing of high-energy and low-energy 
material is suppressed. 

With the centrifugal impeller with radially ending blades 
in question here, what results in the transition region from the 
axial into the radial portion of the impeller is an intensified 
overlapping of the damping effects mentioned (due to rotation, 
or the bending of the channel in the meridian and circumferential 
direction) in the corner of the channel between the housing wall 
and the suction side of the blade channel [6.12, 6.13]. 

Studies by Bradshaw [6.22] on the effect of longitudinal 
channel bending have shown that already at a ratio of 6 /Rk “ 1/80 
(boundary layer thickness / bending radius of the wall) the 
wall friction coefficient Cf, or the distance to the point of 
separation, decreases by 10| in comparison with the case without 
any bending. 

As Johnston and co-workers [^.23] in boti' theoretical and 
experimental studies have shown, an increase in the angular 
velocity of the rotating system likewise causes an increase in 
the turbulence damping effects which, under certain circumstances, 
can lead to a complete relamlnizatlon of the blade boundary 
layer on the suction side. 

In spite of this fundamental finding on the essential 
flow mechanisms in the centrifugal impeller, quantitative 
treatment of this phenomenon is still in its early stages. Even 
for the comparltively simple case of a rotating diffuser with 
radial, straight side walls, the solutions available at this 
time are still unsatisfactory. This is true both for the con- 
sideration of secondary flows (cf. Moore [6.24]) and for bending 
and rotation effects on the struct ux’e of the turbulence (cf. 

Elde and Johnston [6.2^]. 


In view of this situation, the detailed analysis of the /138 

flow measurements made in the centrifugal compressor impeller 
being studied here [6.12, 6,1^] have up to now been limited to 
a qualitative dlscussslon of the observed flow mechanisms. As 
an example. Fig. 6. 16 shows the main secondary flow pattern 
in the centrifugal impeller region. This pattern was derived 
from relative velocity distributions in measuring surface 
(x/Sjjj = 0.87). In this case, for reasons of simplicity — and 



Fig. 6.16. Principle secondary flow pattern in 
the radial portion of the Impeller blade channel 
(with respect to the blade-congruent main flow). 


deviating from th;. usual definition — all velocity components 
transverse to the blade-congruent flow are designated as 
secondary flow. The movement of the core flow in the direction 
towards the pressure side of the blade channel certainly in 
this channel cross-section already contains components which are 
determined by the flow reducing deflection in the region of the 
Impeller outlet (relative channel vortex). This transverse 
movement ends up in the region of the front wall of the channel 
(arrow 1) as a dominating primary vortex. This removes low- 
energy flow material from the channel surfaces passed over and 
feeds it into the wake region. The same is true for a weaker 
secondary vortex in the region of the rear wall of the channel 
and blade suction side (arrow 2). Finally, due to the gap 
between the outer contour of the impeller and the housing, 
additional low-energy material is brought into the wake (arrow 3)* 

By contrast, the effect of friction with the housing wall on 

the measured relative velocity curve remain remarkably small. /139 

Studies on the development of secondary flow intensity [6.12] 
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have shown that in the axial inlet portion of che^entrifugal 
con^jressor impeller (measuring surfaces Q) and 4^ ) practically 
blade-congruent flow predominates (deviation between blade 
angles amlrelative flow direction < ±1®); in measuring 
surface — at the beginning of the radial impeller region — 

appear the first, clearly measurable signs of secondary move- 
ments (A3 < ±5®) which rapidly Increase in intensity up to the 
impeller outlet (A$ * ±25®). 

The principle pattern and the sharp increase in Intensity in the 
secondary flows towards the Impeller outlet are largely in 
keeping with theoretical expectations [6.12]. It is obvious to 
assume that the rapid Increase in the wake region observed at 
the same time is closely connected with this. The following 
model describes this interaction; 

The slowdown in flow in the axial region of the Impeller 
leads to an increase in thickness in the housing boundary layer 
(cf. measuring surface (^ in Pig. 6.15» upper graph). The 
incipient loss-affected region of separation in measuring surface 
4l) also Increases the total enthalpy gradient transverse to 
the direction of flow. As a result of this and as a result of 
the additional displacement effect of the wake, the intensity 
of the secondary flow Increase in the region of the meridional 
bend in the channel. This in turn leads to an increased feeding 
of low-energy flow material into the wake (arrow 1 in Pig. 6.16). 

The turbulence damping effect present in the free plane of shear 
between the Jet region and wake region suppress the turbulent 
mixing of the two regions. This finds its expression in the 
steep, relatively stable velocity gradients at t he edge of the 
wake region (cf. Pig. 6.13, measuring surfaces and (|5> ) 

and in the loss concentration shown in Pig. 6.14. The process 
described explains the rapid Increase in wake regions even in 
such cases where the flow separates relatively late in the 
Impeller flow channel. 

To be sure, the flow development described does not yet 
permit us to say anything about the exact causes of the separation 
of flow. The interpretation of the available findings here /140 

provides only clues on the relative Importance of the flow 
mechanisms in question. The separation begins in measuring surface 
11 in the region of the greatest channel bending. In actuality, 
measurements here show a clear drop in the fluctuation intensity 
of the absolute velocity vector in the region of the front wall 
of the channel and the blade suction side [6.12, 6.14]. This 
indicates the turbulence damping effects described. By contrast, 
in view of the small secondary flow Intensity in this region, the 
collection of low-energy flow material appears still less impor- 
tant. 


Taylor, in [6.26] has discussed In detail the shapes of 
two-dimensional and three-dimensional boundary layer separations. 
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In the first case, the flow separates at a singular point at 
which the wall friction coefficient Cf * 0. What develops down- 
streeiin then are pulsating back currents and high flow losses in 
the region of separation. By contrast, three-dimensional separ- 
ations have has the only criterion [sic] that the flow leaves 
the blade surface; the wake flow proceeds steadily and less loss- 
affected (also see discussion section in [6.1^1]). 

In the present case, the three-dimensional character of the 
total flow, the absence of complete flow stag^tlon at the point 
of separation (Fig. 6.13, measuring surface 
relative stable form of the wake area can be interpreted as 
characteristics of three-dimensional flow separation. 

Perhaps what happens — as suggested In^ig. 6.6 in t^ 
region of the wake — when secondary flows (p and © + ® 
collide, the boundary layers curly up. This forces the main 
flow away from the wall and thus Initiates the separation. It is 
conceivable that a new state of equilibrium is then established 
in the flow, in the process of which the separation point migrates 
upstream into the region where turbulence damping effects and 
weak secondary flows already Increase the tendency towards separa- 
tion of the channel boundary layers. 

7. Theoretical Treatment of the Jet-Wake Flow behind the Impeller /iHl 
Outlet 


The experimental findings discussed in the previous section 
are characteristic of the disturbed flow behind highly-loaded 
centrifugal compressor Impellers. The strongly pronounced Jet- 
wake distributions of the relative velocity at the outlet of the 
individual blade channels rotate with the Impeller and in so 
doing extend radially so far into the bladeless diffuser until 
the velocity differences in the circumferential direction disap- 
pear as a result of mixing and equalization processes. In this 
region the absolute flow in the diffuser is periodically non- 
steady. With the usual assumption of a steady, i.e. axially 
symmetric diffuser flow, it is therefore not possible to ac- 
curately determine the flow mechanisms occurring in the wake 
region of the centrifugal Impeller. 

Starting with a simplified, two-dimensional model of the 
disturbed flow at the Impeller outlet. Dean and Senoo [7.1] 
developed a theoretical solution method with which it is 
possible to determine, at least approximately, the mixing 
processes in the sphere of Influence of the wake region. From 
this can then be drawn Important inferences as to the inter- 
relationship between the impeller and diffuser, which, for 
example, find their expression, in connection with stage design, 
in the selection of the inlet diameter of the bladed diffuser. 
Detailed knowledge of Jet-wake mixing processes is also Important 
for evaluating noise problems and answering questions concerning 
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mechanical strength In connection with diffuser-blade vibration. 

Dean and Senoo set up their theory « which was published in 
i 960 , for two-dimensional, incompressible flow. Subsequently, 

Johnston and Dean [7*2] simplified, in comparison with the 
original Dean-Senoo theory, the calculation of mixing losses. 

In their simplification the mixing process is treated theoretically 
like the Carnot transition. In [6.4, p. 109 f] Traupel 
discusses a comparable method for calculating mixing losses 
behind a centrifugal impeller. The sudden expansion analysis 
given by Johnston and Dean for two-dimensional, incompressible 
flow was expanded by Dean and co-workers [6.5, app. II] to 
two-dimensional flows, while Watson and Ingham [7«3] setup a /142 

solution for compressible flow. 

The measuring techniques developed within the scope of this 
study made it possible for the first time to exactly measure the 
disturbed flow in the diffuser inlet region up to circumferential 
velocities for the impeller of U 2 < 400 m/s. Therefore, for the 
analysis of flow mechanisms in a complex flow the comparison of 
experimental results with theoretical calculations based on the 
theory given by Dean and Senoo presented itself. Of course, this 
required extending the theory to the region of compressible flow. 

In spite of the strongly idealized model, the thus modified theory 
had the advantage that the continuous development of Jet-wake 
mixing can be determined with it. The comparison with measurements 
made at different distances behind the impeller outlet should 
then provide Information on to what extent the theoretically con- 
sidered flow mechanisms affect the actual flow, and from this to 
derive, if necessary, further refinements of the theory.. 

7.1 Expanded Dean-Senoo Theory 
7 . 1.1 Analytical Flow Model 

The mixing and equalization processes behind the 
impeller outlet were considered in the reference system rotatlngg 
with the Impeller. Within the relative system the flow can be 
considered steady and two-dimensional in a plane perpendicular 
to the compressor axis (6/6z ■ 0). In the circumferential 
direction the disturbed real flow is schematically approximated 
by the periodic curve at the impeller outlet shown in Pig. 7.1 a. 

The flow region of a blade spacing is divided into two areas, 
each with a constant relative velocity vector. 

Analogous to the designations "wake” and "Jet” in the 
original work of Dean and Senoo, here, too, the region of 
separated flow is labelled "wake" and the region of undisturbed 
flow is called "Jet." The wake component 62 extends along the 
blade suction side with low relative velocity W 2 y» while the Jet 
region with higher velocity W 2 s fills up the rest of the cross- 
section up to the blade pressure side. /143 
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Pig. 7«la. Analytical flow model at the Impeller 
outlet; b. Schematic diagram of the* development 
of Jet-wake mixing in the bladeless diffuser. 


The two partial flows are not supposed to mix on their way 
through the bladeless diffuser, but rather — as shown schemati- 
cally in Pig. 7.1b with a constant mass flow component in the 
rotating system, they are supposed to flow next to one another. 
As a compatibility condition for the steady relative flow it 
then holds that at the Jet-wake interfaces the static pressure 
and the direction of the relative velocity are equal on both 
sides . 

As Pig. 7.1b Illustrates, the wake conqponent e steadily 
decreases as the radii ratio increases. Strictly speaking, 
therefore, on the same radii the flow direction at the front 
side of the Jet (in the direction of rotation) must be dis- 
tinguished from that on the rear side. However, in this sim- 
plified, Integral flow analysis a uniform flow direction 0 in 
each rad jus is assumed for the total flow. The nonsteady ab- 
solute velocity c, however, is not made unidlre'^tional, as is 
obvious from the velocity vector diagram shown in Pig. 7.2. 

The last-mentioned, limiting condition of the Dean-Senoo 
theory is not included in the expanded model discussed in Section 
7.2 for extlmating the influecne of variable relative flow 
direction in the Jet region and wake region. 

The expanded theory proposed here takes into consideration 
changes in density in the radial direction, but in so doing 
ignores density differences between the two partial flows. 

With respect to the equalization of the disturbed jet-wedee 
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Pig. 7* 2. Velocity vector diagrams In various 
radii ratios to Illustrate, by way of example, 
the exchange of energy between jet regions and 
wake regions. 


velocity distribution, the Dean-Senoo theory takes Into account 
three flow mechanisms: 

1) Friction on the diffuser side walls; 

2) Interface friction and 

3) Reversible exchange of energy between jet and wake regions. 

The diffuser wall friction force per unit area In their 
theory is expressed In the usual manner by a constant friction 
coefficient cf proportional to the square of the absolute 
velocities, whereby this force acts in both regions against /1^5 

the absolute velocity direction in each case: 



8 

2 

p Cg / 2 

(jet) 

(7.1a) 

»T 

8 


(wake) 

(7.1b) 

To take into 

account 

friction at 

the interface. 

equal and 


opposite turbulent shearing stresses tangential to the jet-wake 
boundaries are assumed. 0\er a plane of shear friction co- 
efficient C these stresses are set proportional to the square of 
the relative velocity difference: 

S . C P («s - V*' * (7.2) 


Besides these two loss-affected flow mechanisms, the reversible 
energy exchange between Jet and wake regions only Influences the 
intensity of the mixing of the disturbed flow. 

This type of energy exchange arises under the effect of 
standard stresses between two (nonsteady) partial flows In the 
abscluce system which In the rotating reference system flow 
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steadily next to one another (as everyone knows, in so doing, 
the transfer of work due to normal stresses, in contrast to 
that due to shearing stresses, is reversible^ 

Poa [7.^] has studied a number f Jet engines which 
operate accordl’^g to this principle. In the present case of 
disturbed Jet-wake flow, the velocity vector diagram shown xn 
Pig. 7.2 illustrate, by way of example, this effect: 

If the Jet-wake partial flows shown in the diffuser inlet 
(R/R2 ■ 1.0) were continued with absolute velocities C2S and 
C2T as in the axially symmetrical case — mutually unaffected — 
then in the radii ra^lo R/R2 ■ 1.5 the aosoiute velocities 
and Cm woul occur. In the rotating system tnese would correspond /146 
to the — no longer unidirectional — relative velocity vectors 
wi and wi. However, to fulfill the compatability condition 
between the Jet and wake, both rtgions must interact in such a 
way that the real relative velocity vectors Wc Wip turn cut to 
have a common relative flow direction 0. As Fig. 7.2 clearly 
shows, this only becomes possible with a simultaneous change in 
angular momentum R*ACy in the Jet and wake. In so doing, the 
energy equilibrium for the example shown requires that the energy 
given off by the wake regions equals t:hat picked up by the Jet 
regions. The force which causes these changes in angular mo- 
mentum and which provides for the fact that the real absolute 
velocities eg and ct occur Instead of cX and ci, Is theoretically 
approximated by a sawtooth-shaped curve^ln the'^circumferential 
direction of the static pressure -- Fig. 7.1b — with equal and 
opposite compressive forces on the Jet and wake region. This 
compressive force rotates together with the distrubed velocity 
distribution in the absolute system where it is recorded at a 
determined place in the diffuser as a periodic alternating 
pressure with a fluctuation range q. 

For a theoretical illustration of she reversible energy 
exchange between relative flow lines, the energy equation in- 
troduced above for nonsteady, compressible and frictionless 
absolute flow — Eq. (6.1) — is rewritten for the casf in 
question. With the designations in Figs. 7.1 ard 7.2 and with 
0 as the circumferential coordiante (positive in the direction 
of rotation), it holds for the change in angular momentum along 
a differential path length ds of a relative flow line: 

d<Re^)/dt ■ - l/(wp) • ip/»6 • (1.3) 


In regions with — in the direction of rotation ~ negative 
pressure gradients (e.g. in the Jet region in Fig. 7.1b) the 
angular momentum along the relative flow line increases, while in 
regions with positive pressure gradients the opposite effect 
occurs and thus the loss-free, energy mixing of the distrubed 




flow is encouraged.^ 

As the measurements made at the Impeller outlet confirm, /l47 
the 90®-centrifugal impeller studied here generally possesses 
the higher kinetic energy in the wake. Therefore, in general, 
with the same mean static pressure at the diffuser inlet, the 
wake must perform work on the Jet region in order to mix the 
flow in an axially symmetric manner on its way through the 
diffuser. However, as is confirmed by sample calculations made 
in Section 8.2, cases are theoretically quite possible in which 
the Jet region first performs work on the wake in the diffuser 
inlet region and thus the difference in energy level between the 
two regions is increased even more. With larger radii ratios, 
the wake then transfers more energy back to the Jet than it 
originally received, and in turn the tendency to axially 
symmetric mixing arises. The conditions shown in Pig. 7«2 for 
R/R 2 = 1*5 (with the corresponding pressure curve in Pig. 7 ‘lb) 
are therefore to be understood only as a typical example which 
Illustrates the necessity for the exchange of energy between 
tlie Jet and wake as a result of the required geometric com- 
patabllity within the disturbed relative flow. 

7 . 1.2 Pundamental Plow Equations 

Since according to our hypothesis the Jet and wake 
regions do not mix, the fundamental equation for steady relative 
flow between a position in the diffuser inlet (subscript e) and 
anyplace downstream in the diffuser (without subscripts) can be 
setup separately for both regions. In so doing, all velocities 
are made dimensionless with the circumferential velocity at the 
Impeller inlet U 2 , e.g. Wg * *'^ 3/^2 U = u/u 2 . Symbols are 
given corresponding to the definitions in Pigs. 7*1 and 7*2. 

In addition, Cf stands for the wall friction coefficient, 5 the 
shearing stress coefficient, b the axial diffuser width, z the 
number of blades, p the density. As for the subscripts, S 
stands for Jet, T wake, u circumferential component, and by way 
of simplification R' = R/R2* b' = b/bg and p' = p/pg. 

Continuity Equation ; /148 

a) Jet region 2 R* b (1-e) p Wg sinB s 

2 IT r; bg (l-€g) Pg Wg 3 sinBg , (7-4) 

with A * R' (1-e.) sinB. it follows that 

0 0 0d 0 


1 . Also see the discussion on the corresponding experimental 
results In the diffuser inlet region. Section 6.1.1, and the 
quantitative estimates on reversible energy exchange by Senoo and 
Ishida in [7-5]. 
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1-e * A / (R* Wg sine b' p') , (7.4a) 

e a 1 - A / (R* Wg sinB b' p') . (7.4b) 

b) Wake Region 

2 ir R* b e p SinB = 2 ir R^ b^ p^ sinB^ , (7.5) 

with B a R' e W sinB^ it follows that 

e 6 61 6 

W.J, a B / (R* € sinB b‘ p') (7.5a) 

and with Eq. (7* 4b) we obtain 

Y a / Wg a B / (R'Wg sinB b* p' - A) . (7.6) 


Angular Momentum Equation ; 
a) Jet Region 

2 ir R* b (1-e) p Wg sinB(d(R* U - R* Wg cosB) / dR’J * 

a - 2 s R'^ Rg Cj. (1-e) p Cg (U - Wg cosB) ♦ (7.7a^ 

♦ C z b R* p (Wg - W^)^ ctgB - z b R‘ q / 

On the left side of Eq. (7.7a) stands the change in angular 
momentum flow through the anular volume with the sidth b and 
the radial distance dR’ , while on the right hand side of the 
equation are arranged the moments of the external forces. The 
first expression stands for the moment due to friction on the 
front wall and rear wall of the diffuser, the second stands for 
the moment of plane of shear friction in the 2z jet-wake 
boundary areas and the last expression stands for the torque 
which arises as a result of the pressure difference q between 
the front and rear boundary (in the direction of rotation) of 
the jet region Uf* Fig. 7.1b). By definition, q is positive 
if the static pressure on the front side of the jet is greater 
than on the rear side. 

With Q » q / (p Ug^) finally follows that /149 


(1-c) (Wg ainB(d(R* U - R* Wg cosB) / dR‘] ♦ 
♦ Of R» Rg Cg (U - Wg cosB) / b) « 

* [C s (Wg - W,p)^ CtgB - z Q1 / (2 ir) 


(7.7b) 
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Analogously, the following applies for 
b) Wake: 

- e (W^ sineCdCR* U - R» cosB) / dR* J 
♦ Cj. R* Rg (U - COSB) / b> 

(C 2 (W« - W_)2 ctgB - 2 Q] / (2 ir) 


(7.8) 


Equation of Momentum In the Radial Direction : 
a) Jet Region 


(7.9) 


2 w R* b (1-c) p {Wg sinB[d(Wg sinB) / dR’]- (Wg cosB)^/R') = 

- 2 ir R’ b (l-€) (dp / dR') / - 2 b q ctgB / - 

- 2 ir R* Rg p (1-e) Cg Wg sinB - C 2 b p (Wg - W^)^ ♦ 

♦ 2 ir R’ b p (l-€) (U^ / R' - 2 U Wg cosB / R’) 

On the left side of Eq. (7.9) stands the change In the radial 
momentum flow with respect to the same control volume as In 
the 'ingular momentum equation. The right side of the equation 
takes Into consideration external forces: compressive forces 
(the radial compressive force averaged over time and the com- 
ponent of the difference In pressure q between the front and 
rear side of the Jet)j friction forces on the diffuser walls; 
shearing forces between the Jet and wake regions and. In addition, 
because the calculation Is made In the rotating system, the 
radial centrifugal force and the Coriolis force through the cir- 
cumferential component of the relative velocity. 


With 


P « p / (p Ug*) 


and Q * q / (p Ug^) 


It finally follows that : 


2 w R* (1-c) {Wg 8 inB(d(Wg sinB) / dR’ ] “ w. 

- (0 - Wg 0080)^ / R* ♦ dP/dR* ♦ Cj R, Cg Wg 8inB / b> ♦ 


(7.9a) 


♦ C 8 (W- - W-)‘ 


- 8 Q otgB 


and analogously for the 
b) Wake Region: 

2 ir R* e (W^ 8inB(d(W^ 8lnB) / dR’] - 


( 7 . 10 ) 


- (U - W,p co 8B)‘^ / R’ ♦ dP/dR’ ♦ Cj Rg W^ 8inB / b> - 

- C 2 (Wg - W^)^ « 8 Q CtgB 





If one replaces Q in Eq. (7.9u.) with Q from Eq. (7‘7t>) then it 
follows that: 


dP/dR» 


/ R' - (dWg^/dR') / 2 - 


®f ^2 ^ COS0) / (b sin0) 

5 z [(Wg - W^) / 8in0]^ / [2 ir R* (1-e)] 


( 7 . 11 ) 


and by combining Eq. (7.10) with Eq. (7.8) according to 


dP/dR' 


/ R* - (dW^^/dR*) / 2 - 


- Cj. Rg (W^ - U CO80) / (b 8in0) 

♦ C z [(Wg - W^) / 8in0J^ / (2 ir R* e) 


( 7 . 12 ) 


in so doing, use is made of the following equations for the 
jet region and wake region: 

s (W sin0)^ ♦ (W co80)^ and 
dW^/dR* s 2 W 8in0td(W sin0)/dR’] + 

4- 2 W co80(d(W co80)/dR'] . 

If one eliminates dP/dR* from Eqs. (7.11) and (7.12), then 
this gives us: 

(7.13) 

d(Wg^ - W^^)/dR’ ♦ C 2 (Wg - W^)^/ [ir R‘ e (1-c) sin^01 ♦ 

♦ 2 Rg fCg (Wg - U CO80) - (W^ - U co80)]/(b sin0) * 0. 

Seting the right side equal to Eqs. (7* 7b) and (7.8) leads to: 

(7.1^) 

Wg (1-e) B(R’ Wg co80)/dR»l ♦ W^ e 61(R' W^ co80)/dR') - 

- 2 U tWg (1-c) ♦ el - 

- R’ Rg (Cg (l-e)(U - WgCO80) ♦ c (U - W^co80)l / 

/ (b 8in0) ■ 0 . 
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Since the absolute velocity C — for the Jet region and /151 

wake region — can be expressed by means of the cosine equation 
by W, U and 3 : 


C = (W^ ♦ - 2 U W cosB)^^^ , • (7.15) 


then Eqs. (7.13) and (7.l4) still contain only the four unknowns 
Wo, Wm, e and p. uf these, e and Wm, for example, can be replaced 
by Eos. (7.4b and (7.6). The equations derived above — (7*13) 
and (7.l4) — thus represent a system of two coupled differential 
equations of the first crder for the dependent variables W 3 (R') 
and p(R'). This system of equations Is solved according to the 
Runge-Kutta method (for differential equation systems) discussed 
In detail In [ 7 . 6 ]. 



As starting conditions, the variables WeS, Wg^, Pg, eg 
and R'e are given. Additional Information on the static pressure 
Pe/Pto and the total temperature in the diffuser T^ii are given 
for the calculation for compressible flow. In addition, the 
friction coefficients Cf or Cf(R) and 5 , the Impeller blade 
number z and the effective width^ b(R) are known. 


The Runge-Kutta approximation method consists of two coupled 
calculations running In parallel in accordance with the speci- 
fication given below for differential equations. Progressively, 
In each case, two function values are calculated at four 
equidistant reference points between the known initial values 
(subscript A) and the end values being sought (subscript E) 

(total Interval AR' = R'e“R*a “ 0.01). Prom the function 
values determined in this way two mean values are then formed 
for calculating the new final values . 

With e = Eq. (?.4b), y = Eq. (7.6), 0311 = Eq. (7.15) and 
the following notation: * 


* (1-e) e sin^B (1 -y^)^ + > 
Cg * (l-€) sin^B (1-Y^) ♦ Y^ * 
Cj * Cg. (Wg - R* cosB) - (W^ - R* cosB) , 



1. Section 7 . 1.4 contains data for taking Into account the boundary 
layer displacement thickness. 
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C|, s Cg (l-€) (R‘ - Wg cosB) ♦ c (R* - cosB) , 

Cj * e (1-e) (1-Y^)^ ♦ » 

Cg s e (1-Y^) ♦ Y^ * 

Cj a 1 - c (1-Y^) » 

Cg = 1-e (1 -y) and 

D a (dpVdR.') / p* 

it holds that: 



dWg/dR' a 

- Y^ cosB (Wg cosB Cy - 2 R* Cg) / (R* C^) 

- Wg 2 C (1-y)^ Cg / C2 s (1-e) R* sin^B C^l - 

- c^ Rg (c Cg Cj - Y^ cosB C^) / (Wg b sinB C^) - 

- Wg sin^B Y^ Cy D / Cj ; 


( 7 . 16 ) 


(7.17) 


dB/dR* * 

(Wg sinB cosB Cg - 2 R* sinB Cg Cg) / (Wg R’ C^) 

- C 2 (1-y)^ (1-Y^) cosB / (2 ir R* sinB C^) - 

- Rg te (1-e) (1-Y^) cosB Cg ♦ C|j CgJ / (Wg b C^) 

- Y^ sinB cosB D / C^ 


If one of the limiting conditions e < 0.01 or 6 = 
e (1 -y) < 0.001 is fulfilled, then the disturbed flow can be 
regarded as "mixed" [equalized], and e is set equal to 0 and 
Y is set equal to 1 in Eqs. (7* 16) and (7.17). 

In Eqs. (7.16) and (7.17) D expresses the variation in 
density in the radial direction. The expanded Dean-Senoo 
theory presented here takes into consideration the compressibility 
effects in an iteration cycle after each step of the Runge- 
Kutta method. In so doing, two equations of the computing method 
given by Traupel ([6.^1], Eqs. ^1.5(16) and ^1.5(21)) are Introduced 
for the compressible, axially symmetric diffuser flow. As /153 

written here, these equations are expressed as follows; 
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, (ie-l)/ie 

/ ip/» . K/lK-1) KPe'P*)' ' - 11 . 

p« 


p£ * **A^ (k-1)/(2ic) - 7g^)u2^/Pgl ^ . (7.19) 


During the first Iteration loop, since D = 0 in Eqs. (7.16) 
and (7.17) the change in density is Ignored to begin -with. 

With the aid of Eq. (7.11), the results found in this way per- 
ini*- the computation of the left side of Eq. (7.18). From this 
follows the static pressure at the end of the Interval pg. 

After determination of the moment urn- averaged absolute velocity 
C of the total flow, Eq. (7.19) makes possible a first approxi- 
mation for the density at the end of the step Pg. For the 
following iteration loops D is then set equal to 

D s (doVdRM / p* « 2 [(Pg - P^)/AR’1/(P^ ♦ Pg) (7.20) 

and thus the radial change in density is approximately taken 
into account. 

Finally, the total pressure in each of the step points is 
calculated from: 


^ « P Cl - (ic-l)/2 Ug^/Cic R j ( 7 .; 


whereby R stands for the gas constant of the air which is 
regarded as an ideal gas. 

Computer results regarding the course of the Jet-wake 
mixing proces are presented in section 8.2 and 8.3. 

7.1.4 Determination of the Wall Friction and Interface 
Friction Coeffic'.ents and the Boundary Layer Displace- 
ment Thickness 

The Dean-Senoo theory contains two unknown parameters 
— the wall friction coefficient cf and the Jet-wake interface 
friction coefficient C. All other paramters are given by the 
conditions in the inlet of the unbladed diffuser and by the 
geometry of the flow channels. In view of the complicated flow 
conditions in the diffuser inlet region, these coefficients can- 
not be determined theoretically with sufficient accuracy. There- 
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fore, one is directed to direct measurements or, as in the case 
of C, to adopting values which have been determined for visible 
flows . 

In Fig. 7*3 are plotted experimentally determined curves of 
cj* as a function of the diffuser radii ratio R/Rp. For the 
centrifugal compressor stage studied here, a typical c^ curve 
(DFVLR) was determined according to the method given by Traupel 
in ([ 6 .^ 1 ], p. 11 3f) for the test point at 1^,300 rpm, A = 5*31 
5.31 kg/s. In so doing, the basis of the Cf calculation is 
formed by_”conventional," time-average measurements of the wall 
pressure p, the total pressure p^, the total temperature ft and 
the flow direction a as a function of R/R 2 * 



1.0 1,1 1,2 1,3 1,4 1,5 

R/Rj — 


Fig. 7.3. Experimentally determined curve of the 
wall friction coefficient Cf as a function of the 
radii ratio R/R 2 of an unbladed centrifugal dif- 
fuser. (Re number at the diffuser inlet established 
with diffuser width). 

Key; A) Measurement C) Hot film 
B) rpm 


; For purposes of comparison, cf values were plotted which 

I were measured by Senoo and Ishida [7*8] in a parallel-walled 

I diffuser behind an impeller with backward bent blades. In 

the region of the approximately axially symmetric diffuser flow, /155 
R/R 2 > 1.15, the Cf values shown at the given Reynolds numbers 
largely agree with the values for technically rough tubing. 

I 

t 

It should be noted that the Cf values of Senoo and Ishida 
were measured directly by means of a Preston tube, or nonsteadlly 
by means of a hot-film probe which was Inserted transversely into 
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the diffuser wall. While for R/Rp > 1.15 the DPVLR curve, cor- 
responding to the varying Reynolds ‘^numbers, clearly lies below 
the measurements made by Senoo and Ishlda, the DPVLR results for 
R/R 2 = 1.15 increase comparitlvely more sharply up to the 
impeller outlet. This difference is presumably a result only 
of measuring techniques, since the total pressure measurements 
in the diffuser inlet region used to determine the wall 
friction coefficient also involve the momentum loss which arises 
due to the mixing of the disturbed flow in the diffuser inlet 
region . 

Por calculating friction losses, Traupel ([fv.4], p. 125) 
recommended Cf = 0.0045-0.006. Dean and Senoo C7*l] and later 
Senoo and Ishida [7*5] used Cf ■ 0.005, independent of the 
radii ratio R/Ro. The shaded curve for Cf (R) in Pig. 7.3 
is based primarily on the calculations in sections 8.2 and 8.3. 
In one example, Cf is set equal to O.OO65. 

Por the interface friction coefficients 5, Dean and Senoo 
[7.1], in their parameter study, gave the values 0.094 and 
0.037 without explaining how these numerical values were 
derived. 

Shearing stress measurements for "free” turbulent shear 
flows are currently available only for a few simple sample cases. 
In ([3.9], Pig. 6-21) Hinze cites measurements by Corrsln made 
in a round free Jet (2.5 cm in diameter) and in (C3.9] Pig. 6-6) 
he cites measurements by Townsend made in the wake region of a 
circular cylinder (1.59 mm in diameter). In these cases, maximum 
shearing stress coefficients c of 0.03 (free Jet) and 0.12 
(cylinder wake) respectively were observed. These values 
coincide well with the theoretical computer results of Rotta 
(C7.9], Pig. 74) for the maximum values of the Reynolds shear 
stresses on a level free Jet. 

Thus the order of magnitude of the c values originally 
used by Dean and Senoo are well confirmed by these results. 

Of course, it must still be taken into consideration that the 
simple cases cited above for purposes of comparison can cer- 
tainly be applied only to a very limited extent to Jet-wake 
shearing surface friction in the rotating system. Therefore, 

Dean also assumes in a qualitative consideration C7.10] — which 
is essentially based on the varJable effect of Coriolis 
forces in Jet and wake regions ~ that the friction coefficient 
on the rear side of the Jet (in the direction of rotation) is 
strongly Increased (5 < 0.5), while on the front side of the Jet 
the turbulent mixing is drastically reduced (c * 0). Therefore 
he proposes an average Jet-wake friction coefficient of 5 “ 0.25 
in place of the value of c “ 0.094 used in the original study. 


1. Pig. 7.3 shows the mean values over time determined from this. 
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In the friction-affected diffuser flow, boundary layers 
occur on the front wall and rear wall of the diffuser which, 
due to their displacement effect which increase with downstream 
distance, cause a decrease in the flow cross-section with 
respect to the geometrical cross-sectional contour. This dis- 
placement effect is characterized by the boundary layer dis- 
placement thickness 6*. Since no boundary lay^r measurements 
are contained in the test program discussed here, we will refer 
to data given in the literature. 

Fig. 7.^ shows experimental results obtained by V.D. Voorde 
and Bos [6.3] in a parallel-walled centrifugal diffuser with a 
superimposed, axial-radial deflection with purely radial flow. 

A striking feature here is the sharper increase in boundary layer 
thickness on the front wall of the diffuser. At the same’ 
time. Pig. contains results of a theory^ given by Scholz 
[7.11] on the Increase in boundary layer displacement thickness 
in a plate diffuser with level source flow. The good agreement 
between theory (broken line) and experimental results from 
[6.3] shows that the simplifying theoretical assumption of a 
linear Increase in the displacement thickness is Justified. 


mm 

t 

Measurement [6.2 
0 .front wall V 



1 , 

A rear wall R 



5-^ 
0 z 

0 Sum V+R 
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Theory [7.111 >* " 
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0 * ‘ 

A 

1 1 

1 
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Fig. 7.^. Increase in boundary layer displacement 
thickness in an unbladed, centrifugal diffuser, 

with respect to the radii ratio of the DFVLR stage 
(Re/Ro - 1.017). 

Theoretical prediction [7.11]: 


a) — — — for the test conditions in [6.3]* 

b) for DFVLR measurements at 1^1,000 rpm 
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1. This theory is based on the "Quadrature Method for Calculating 
Turbulent Friction Layers" by Truckenbrodt . 
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The solid line marks the theoretical prediction for 
operations conditions such as were present with the DPVLR 
measurements at 1^1,000 rpm. This curve of the boundary layer 
displacement thickness was used as a basis for all of the sample 
calculations with the effect of wall friction discussed in 
Sections 8.2 and 8 . 3 . 

7.2 Expanded Model ftr Estimating the Influence of the 

Variable Relative ti’low Direction in the Jet Region and 

Wetke Resion 

According to the Dean-S" noo Theory it had been assumed up 
to now that Jet and wake regions possessed a uniform relative 
flow direction 8. However, in the actual flow the Integral 
values 8s ®T dlffsr distinctly. For example, evaluation 
c our experiments in the measuring plane directly behind 
the impeller outlet shows variations in direction on the order 
of 10« (cf. Fig. 8.22). 

In order to estimate the Influence of this effect which is /158 
not taken into account in the Dean~Senoo theory, a simple, 
supplementary flow analysis — following a proposal by Van Le 
[ 7 . 7 ] — is carried out for the Jet-wake flow model shown in 
Pig. 7 . 5 . 




Pig. 7 . 5 . Plow model with variable relative flow dir- 
direction in the Jet region and wake region. 


Between a place with the Jet-wake velocity distribution in the 
duffuser inlet region (subscript e, or S (Jet) and T (wake)) and 






■f 


the radius, where the disturbed flow has become couplet ely mixed, 
equations- for the mass flow, the angular momentum and the kinetic 
energy were setup. Effects of friction and compressibility are 
not taken into consideration. 

On the basis of this purely Integral approach to the mixing 
of the Jet and wake — between the disturbed initial state and 
the fully mired final state — the discontinuity at the Jet-wake 
boundary surfaces in the model. Pig. 7.5» is also permissible. 

All of the variables used in this analysis are defined in 
accordance with the graph shown in Pig. 7*5* 7*2 and 5.2. 

Continuity Equation ; 


( 7 . 22 ) 

2 w (R b)^ t(l-e) Mg sinSg ♦ * 2 » R b C^. , 

R b Cj. / (R b Cj.g)^ * (W) |1 ♦ c Y «in0j/[(l-«) sin0g]}, 


Por the parallel-walled diffuser it is valid that b » b^, 
and for Rb (equal in area) ■ (Rb)^. 

Angulam Momentum Equation; /159 

(»-€) (B* b ♦ c (B* b . B* b C„. (7.23) 

B* b C, C„ / t(W) (B* b C,3 C^),I • 1 . . 

/ Kl-C) C^J 


with Eq. (7.22a) it follows that 


(7.23b) 

' R C^ / CR •* ^ ® Y f I(l“C)sin0g 

/ (1 ♦ c Y tin0f / C(l-t) tinRgn . 


Finally, using the velocity couponents Cj, from Eq. (7.22a) and 
Cu from Eq. (7.23b), if we calculate the kinetic energy flow of 
the absolute velocity C: 

t ■ £ NMSt * C*/2 

then we obtain the following as an 


•r 






123 


Equation for the Kinetic Energy ; 


«-e) (R b Cj* / 2), * c (R b C,® / 2)^ . (7-2'') 

R b Cj. C* / 2 , 

R b C C* / t(l-e) (R b c , C *).l • 

^ x*3 s 6 (7,24a) 

1 ♦ e Y 8iri0y / C(l-c) 8in0g Cg-J , 

and setting In Eq. (7*22a): 

• « (C / Cg)^ {1 ♦ e Y sinBy / l(l-c) sinBgJ) / (7.2ilb) 

/ |l ♦ c Y sinB^ / l(l-e) sinBg Cg^l} 


The ratio t/t^ given with Eq. (7.24b) thus defines the 
upper limit value of the mixing losses of the partial flows 
discharged from the Impeller with various relative flow directions 
Results of an numerical parameter study using the above equations 
are presented In Section 8.3* 

8. The Behavior of the Jet-Wake Plow in the Unbladed Centrifugal 
Clf'^use'r 

8.1 The Three-Dimensional Development of the Real Flow In 

the Inlet Region oi^ the Diffuser 

In Section 6 we discussed In detail the disturbed velocity 
distributions measured Immediately behind the Impeller outlet 
and Its origination from a flow separation In the Impeller. 

There we also showed that the observed, fundamental flow pro- 
cesses In the entire range of the performance graph studied 
between 10,000 and 18,000 rpm are largely similar so that the 
following, comprehensive discussion of Jet-wake flow development 
In the diffuser Inlet region can be restricted to measurements 
obtalned-at 14,000 rpm. In general, in so doing the curve for 
the operating point In the region of the optimum line, A/^/6o • 
5.31 kg/s. Is described In detail by results obtained In 
measuring planes I-IV (corresponding to R/R2 “ 1.017, I.039, 

1.089 and I.15I; cf. Section 5.1}* Supplementing this Infor- 
mation, the flow condition In measuring planes I and IV In the 
operating point near the stopper border and pump border Is 
compared. 

As In Section 6, the results are plotted for the outlet flow 
region of an Impeller blade channel. This flow region Is 
graphically simplified as a rectangle. In the real flow the 
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interfaces (cf. Section 5.3) be;*ween two adjacent flow regions 
are twisted in three-dimensions, as clearly Illustrated in 
Fig. 8.1 These interfaces are deter ^ ned from the position of 
the blade wake shifts in measuring planes I-IV. For the three 
operating points analyzed — a - c — the graphs of inter- 
face development are largely congruent so that It appears 
permissible to compare the flow curves with one another for the 
va-r*lous operating points in the simplified rectangular graph. 
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Fig. 8.1 Development of the interface between the flow 
regions of two adjacent blade channels in the diffuser 
input. Operating points at n//So ■ 14,000 rpm, 

Ih /Qq/6q: a 4.53 kg/s, b 5.3l kg/s, c 6.07 kg/s/ 


8.1.1 Static Pressure and Tota?. ?ressure 

Fig. 8.2 shows the static pressure curve Interpolated 
from the nonsteady pressure measurements mad» on the front and 
rear wall of the diffuser. In the circumferential direction, 
the flow region of an impeller blade channel is bounded by the 
interfaces mentioned. 

First of all, a striking feature is the rapid reduction in 
the pressure gradient over the channel width z/b (front wall of 
the diffuser z/b « o, rear wall z/b ■ 1) between measuring planes 
I and II (R/Rp * 1.0x7-1.039). This is explained by the transition' ^ 
from the meridional channel deflection (in the in^)eller) into 
the straight centrifugal diffuser. In tne circumferential 
direction, the variationa in static pressure offset each other 
almost con 4 >letely only up to measuring plane IV. 

1 

i 

1 
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Fig. 8.2. Static pressure curve 
P 2 (y,z) in the diffuser ■'nlet 
region (measuring plane'' I-IV), 
with respect to the total in- 
take pressure p^o* 

Operating point: 
n//5^ = 14,000 rpm, 

* 5,31 kg/s. 

in total pressure is clearly 
observed in comparison with the measured values determined at 
the same spot in plane 1 (while, of course, the total pressure 
mean value clearly decreases between the two planes, cf. 8.l8). 

As an example, the Ap^ curve in Fig. 8.5 — which was 
determined from the local total pressures in planes I and III: 


The corresponding development 
of the total pressure ratio 
Pt8(y ,z)/Pto 1*^ measuring planes 
I-iV of the diffuser inlet region 
is shown by the plotter diagram 
in Fig. 8.3. Immediately behind 
the impeller outlet what is 
first of all striking in the 
transition from measuring plane 
I to II is the rapid mixing 
of the total pressure peaks /162 

in the blade wake region. Further 
downstream, in measuring planes 
III and IV at y/t = 0.7 and 
z/b = 0.5-0. 7, a distinct total 
pressure maximum is then formed 
in the flow channel. These 
observations are confirmed in 
Fig. 8.4 which compares 
measurements made in planes I 
and IV for operating points near 
the pump limit and absorption 
limit . 

In particular, in this con- 
nection we should emphasize the 
at first surprising fact that in 
the region of this pressure 
maximum in plane IV, an Increase 


Ap^(y,z) = tp^^jy(y,z) - Pt^i(y,2)l/Pti(y,z) 


( 8 . 1 ) 


shows a clear increase in pressure for the development of total /I 63 
pressure presented in detail in Fig, 8.3 in the suction half of 
the channel. By contrast, the total pressure losses are con- 
centrated in the corner formed by the diffuser front wall and 
the pressure side and on both sides of the flow channel inter- 
face (y/t = 0 or 1). This pronounced increase in total pressure 
in parts of the flov; channel cannot be Interpreted either as 
an effect of the turbulent mixing process or as an effect of 
a secondary movement within the dlstrubed Jet-wake flow. To be 
sure, a flow line shift would be conceivable, in principle, so 
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Fig. 8 . 3 . Plotter diagrams of the total pressure dis- 
tribution measuring planes I-IV. 

Operating point n/ “ 14,000 rpm. 


* ‘^^^0*5*31 kg/8. 


that the "pressure mountain" which appears in Pig. 8.3 at R/R 2 * 
1.107 in the region of the front wall of the diffuser (y/t = 0.2) 
is shifted into the area of maximum total pressure at R/Rp = 1.151. 
Arguing against this, however, is the almost unchanged outer 
contour of the wake region in both measuring planes (cf. Fig. 8.14) 
and the existing ratio of the circumferential component to the ^ 
radial component of the relative velocity which, under favorable 
assumptions, allows a flow line shift of less than 0.2* t on 
the flow path between planes I and IV (cf. Pigs. 8.15 and 8. 16). 
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Fig. 8 . 5 . Change In total pressure between 
planes I and IV. 

Operating point **^*^®^ = 14,000 rpm. 



Pig. 8.6. Isobar graph of the total pressure 
ratio 03 flow region of a blade 

channel in'^the diffuser Inlet (vertical distance 
with respect to the spacing at R/R 2 =1* Increased 
by a factor of 3*2). 

Operating point; n//ST = 14,000 rpm, 

A ■ 5.31 kg/s 
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(y/t s constant). By contrast, in the middle of the channel 
(z/b = 0.5) the total pressure largely follows the total 
enthalpy development, for which applies the connection described 
in Eq. (6.1) with the change in static pressure over time. 

In actuality, what results here in the region of y/t = O.^l-O.B — 
corresponding to the increase in static pressure per unit time 
shown in Fig. 8.2 — is an increase in total pressure along the /166 
relative flow lines, whereas, for example, from R/R2 = 1.089, 
a distinct drop in pressure follows a negative dp/dt close 
to the pressure side (also see the discussion on the q/pto 
curve in Pig. 8.20). 

The above figures clearly illustrate an exchange of energy 
between areas with a different relative velocity level. 

Undoubtedly this mechanism also Influences the mixing of the 
disturbed flow in the diffuser inlet region. Senoc and Ishida 
discuss comparable findings behind an Impeller with backward 
bent blades in [7*5] . There the attempt was also made to 
quantitatively analyze the reversible exchange of energy based 
on Eq. (6.1). To be sure, this showed considerable dispersions 
in the numerical evaluation. Because of the pronounced wall 
friction effects in the diffuser which was only 14 mm wide, the 
increase in total pressure was at most only Ap^ « 0.01 (cf. [7-5], 

Fig. 11) between R/R2 = 1.014-1.1, as against Ap^ j^ax ^ 

Fig. 8.5. * 


8.1.2 Total Pressure Fluctuation and the Effect of 

Flow Line Bending on the Degree of Mixing 

Further data on the development of the mixing process, 
especially on the degree of mixing in the disturbed jet-wake 
flow, can be derived from the distribution of the total pressure 
fluctuation coefficient f, defined according to Eq. (6.2), in 
the various measuring planes behind the impeller outlet. 

First of all, the fluctuation distributions in planes I and 
IV are compared for two operating points in Fig. 8.7* compari- 
son to the conditions immediately behind the impeller outlet 
(R/R2 ® 1.017), in plane IV (R/R2 = 1.151), we observe the 
following; 

- only small changes in the fluctuation intensity in the 
Jet region, i.e. the pressure side of the channel; 

- a distinct increase in f in the core of the wake 
region; 

- an increase in the region with a high fluctuation /167 
Intensity in the direction of the rear wall of the 
diffuser and the suction side boundary of the flow 
channel. 
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8.7« Distribution of ths total pressurs 
fluctuation coefficient f in measuring planes 
I and IV. — 

Oper^ing point at = 1^,000 rpm, 

* ® ^S/s, © 6.07 kg/s/ 


same phenomena also show up in the detailed illustration /I 67 

of the fluctuation development in planes I-IV, as presented in 

Fig. b.a for the operating point at 14,000 rpm, 5-31 kg/s. In 
addition, in the transition from plane I to II, we observe in 
the entire flow channel a distinct decrease in fluctuation in- 
tenslty before up to plane IV a higher level of total pressure 
1 luc t u3it 1 ons fGS 1 3,b llsh6s Itsslf *•— 3 .S 2 Llr 63 .dy In Fig. 8.7, 


In order to verify that this development does not result 
from the curve of the reference total pressure p^^(y,s), in 
Fig. b.9 the area-averaged value of the absolute fluctuation 
Intensity f» (numerator in definition equation Eq. (6.2)) was 
total intake pressure p^q* The graph over the 
radii ratio R/R 2 likewise shows a distinct minlmum^ln the 
region of plane II (R/R 2 « 1.039). 


At the present time, an exact reason for this fluctuation 
curve based on fluid flow physics cannot be given. This is 
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Fig. 8.8. Distribution of the total pressure 
fluctuation coefficient f in measuring planes 
I-IV (legend, see Fig. 8.7b). 

Operating point: n//oJ^ “ 14,000 rpm, 

A /eJ/«o * 5,31 kg/8. 

primarily due to the possibly very different causes for the 
fluctuation in total pessure — such as nonsteady relative flow 
or axially asymmetric diffuser f lov. . To be sure, certain es- 
timates on the development of turbulence along the channel or 
Jet-wake interfaces in the disturbed flow field can be made. 

The Influence of normal acceleration on the turbulence 
structure of planes of shear within the Impeller flow was already 
discussed in Section t'.4.2. Obviously, the same mechanisms are 
also to be found in a consideration of the relative flow in the 
diffuser inlet region. 

In a plane perpendicular to the rotor axis the following 
equation is regarded as a stability clrterlon for a plane of 
shear under the effect of Coriolis acceleration: 


S^, « 2(r/R2>^ - (r/Rg) (3Wj/3r) , (3 2) /169 
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Pig. 8.9. Radial curve of the 
area-averaged, absolute fluctu- 
ation coefficient 
Operating point * 

li|,000 rpm, ^ 

A /9^/«0 * 5,31 kg/8. 


or the following equation under 
the effect of normal acceleration 
due to relative flow line 
bending: 


(Wj/U2) (r/Ug) (3Wj/3r) ♦ (Wj/Ug) 

and for the overlapping of both 
effects : 


with 

Wg = planes of shear - 
^ relative velocity 

r = bending radius of the 
plane of shear and the 
measuring plane in 
question 


Uo ® circumferential velocity of the impeller in the 
outer radius R2* 


Positive values for Sq, Sj^ and S in each case imply turbu- 
lent damping under the effect of Coriolis or normal acceleration 
or their overlapping. ConverseD.y, negative values indicate 
turbulence excltalon, l.e. Increase mixing in the region of the 
plane of shear. 

S corresponds to the "Brunt-Valsala frequency 
given by Bradshaw in [6.22] for the two-dimensional case 
assumed here, i.e. the mean shearing stress, flow line bending 
and rotation lie in one plane. Equations (o.2)-(8.4) were also 
derived in detail by Hill, and by Senoo and Ishida in the dis- 
cussion section of [7*5]* 


In the following estimate cf the effect of normal acceler 
accelerations on the plane of shear stability behind the impel- 
ler in question here, the following, partially simplified, 
conditions are to be noted: 


1. The flow flexd was evaluated for the test point 
at 14,000 rpm and a mass flow of 5*31 kg/s. 

2. Plane of shear flow lines are considered at z/b ■ 

0.3. 
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3 . The results In a region close to the impeller — 
based on measurements In planes I and II — and 
in a region remote from the impeller (planes III 
and IV) are compared. 

The conditions in the Jet-wake Interfaces on the 
front side of the jet SVS (y/t « 0.5) (in the 
direction of rotation) and on the rear side of 
the Jet SR3 (y/t == 1.0) are considered. 

5 . Both interfaces are based on the same bending 
curve which is taken from Fig. 8.1 b at z/b = 

0 . 3 . 

6. The data on the relative velocity curve come from 
an exact digitalization with 101 values per 
spacing. Fig. 8.15 shows the relative velocity 
distribution only at selected points. 

The results of this e.^tlmate listed in Table 3.1 show, 
first of all for the region close to the Impeller (1.0 < R/R 2 < 

1 . 05 ) that the turbulent mixing of both interfaces is dampened 

under the dominating effect in the case in question of 

Coriolis acceleration (IVS) or normal acceleration (SR3). Whether 

the perceptible drop in total pressure fluctuation intensity at 

R/R 2 = 1.039 (Figs. 8.8 and 8.9) is connected with this cannot, /171 

to be sure, be answered unamblgously with these results. 

Afterwards, in the Impeller- remote region (1.05 < R/R 2 1*2) 
with a simultaneous decrease in the flow line bending in each 
case, Che effect of the Coriolis acceleration predominates. 

This causes severe damping on the front side of the Jet, but on 
the rear side of the Jet it causes an excitation of turbulent 
exchange . 

With respect to the disturbed Jet-wake velocity dlstrlbtion, 
from this it follows that the equalization due to the turbulent 
mixing of the two regions takes place primarily in the area of 
the rear of the Jet starting at R/R 2 * 1.05. By contrast, the 
velocity gradient at the front of the Jet — starting at the 
impeller inlet — becomes increasingly stabilized under the 
effect of Coriolis acceleration. 

In the meridian , lane within the In^eller, the radial 
deflection primarily causes damping of the turbulent exchange 
movement, while behind the impeller outlet a renew increase is 
to be expected with insignificant flow line bending. 

These expectations, derived from a simple estimation, are /172 
surprisingly well verified by the test results. In actuality, 
an Incrase in total pressure fluctuation in the direction towards 
the corner formed by the suction side of the channel and the 
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TABLE 8.1 EVALUATION OF THE PROBLEM OF JET-WAKE INTERFACE 
STABILITY. OPERATING POINT: n/vOZ = 14,000 RPM, ih /^/6r, = 

5.31 kg^s ^ ^ 


z/b “ 0.3 Con^mted data Results 


measur- 

circumfer- 

r 


r dw. 




ing 

ential 


Jl 


Sc 


S 

planes 

position 

Rg 

“2 

U 2 3r 

c 


I-II 

SVS 

0,30 

0,37 

-2,36 

0,89 

-0,74 

♦0,15 

SRS 

0,30 

0,35 

2,29 

0,51 

0,92 

♦0,41 

III-IV 

SVS 

0,75 

0,39 

-4,70 

4,65 

- 1,68 

♦2,97 


SRS 

0,75 

0,42 

4,85 

-2,51 

1 

2,21 

-0,30 


SVS » front of Jet SRS = rear of Jet 


rear wall of the diffuser can be seen in Figs. 8,7 and 8.8, and 
also the equalization of the disturbed velocity distributions 
shown in Figs. 8.12 and 8.13 illustrates the influence of the 
flow mechanisms in question. 

8 .1.3 The Velocity Field in the Absolute and Relative 

System 

In order to simplify the representation of the complex 
mixing processes within the disturbed flow field the development 
of the absolute flow in the diffuser inlet region is shown 
separately for the flow direction a and the radial velocity com- 
ponent Cv. 3 . The direction distribution in planes I-IV can be 
seen in Fig. 8.10 for the operating point at 14,000 rpm and 
5.31 kg/s. On the other hand. Fig. 8.11 compares the conditions 
in planes I and IV for the operating point near the pump limit 
and absorption limit. In an analogous manner, Figs. 8.12 and 
8.13 show the development of the velocity distributions CJ. 3 /U 2 . 

The radial development of the velocity field between planes 
I and IV (Figs. 8.10 and 3.12) clearly Illustrates the Jet-wake 
mixing process. Already we clearly see at the transition from 
plane I to II a decrease in the steep velocity and angle gradients 
over the diffuser width and at the rear of the Jet (y/t » 0 or 1), 
while the gradient at the front of the Jet ,!y/t « 0.4-0. 5) even 
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R/Rg * 1,151 R/Rj * 1,089 



Vorderwand * front wall 

Fig. 8.10. Isocline graph of the absolute flow 
direction o in planes I-IV. 

Operating point: n//eZ » 14,000 rpm, 

A ■5-31 kg/s." 

in plane IV is still conspicuously evident. This velocity curve 
underlines the effect of the mixing mechanisms discussed above in 
section 8.1.1 and 8.1.2 — i.e. reversible energy exchange and 
turbulent mixing. These suggested that the wake region replen- 
ished Itself from the region of the Jet flow close to the rear 
wall of the diffuser and the suction side of the channel. 

These observations are verified in Pigs. 8.11 and 8. 13 
where the experimental results for the operating point close /1?4 

to the pump limit and absorption limit are conqpared. A comparison 
of the velocity distribution in plane IV (R/R 2 • 1.151) shows 
that in the operating point with the greatest mass flow — 

6.07 kg/s — the distrubance in the circumferential direction 
is also still clearly perceivable in the vicinity of the rear 
wall of the diffuser, while for the maos flow of 4.53 kg/s 
a largely evened out velocity curve has been established. 

In addition, for all of the operating points shown we see 
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Vorderwand = front wall 


Fig. 8.11. Isocline graphs of the absolute flow 
direction a in planes I and IV. 

Operating point at = 1H,C00 rpm, 

A ^ ^*53 kg/s, b 6.07 kg/s 

that between planes I and IV there occurs a shift in the mass 

flow towards the rear wall of the diffuser. As a result of 

this, the boundary layer, generally visible at R/R 2 " 1.017 

close to the rear wall of the channel (z/b « 0.9) is rapidly /I 75 

replenished. 

Fig. 8.1^ also shows this in another form. Here for the 
three analyzed operating points at 14,000 rpm, the wake dis- 
tributions in planes I and IV were compared. The shape of the 
wake regions was determined with the aid of definition equation 
( 5 . 60 ) from the corresponding velocity curves CJ. 3 /U 2 . This 
also gave the "wake components" along the rear wall of the 
diffuser shown in plane I. Once again the graph refers to the /176 

flow region of an impeller blade channel. Of course, in the 
simplified rectangular form of the graph, the channel division 
and channel width here were adapted to the existing geometric 
conditions in planes I and IV. By this means, we obtain an 
impression whic‘ approximately represents the area and shape of 
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a/Rg • 1,151 R/R 2 • 1»089 



Vorderwand ■ front wall 


Fig. 8.12. Plotter diagrams of the velocity 
distribution planes I TV. 

Operating point: *»//^ » 14,000 rpm. 

± - 5.31 kg/s. 
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Fig. 8.14. Distribution of wake 
regions in planes I and IV. 
Operating points at n/ZoT s 
it. yS~/A rpm, 0 

m r^Qf^Q . Q 5 QY icg/s, 

© 5-31 kg/s, 0 4.53 kg/s. 


the wake development in the 
diffuser Inlet region. 

Finally, a striking 
feature of the wake regions is 
the fact that the Jet-wake 
separating line in the region 
of the front of the Jet (y/t == 
0.4-0. 6) lines up increasingly 
parallel to the rotor axis as 
the mass flow decreases. Fig. 
8.14 a->-c. This expresses the 
variable influence of the 
secondary flows on the flow 
field in the impeller (cf. 

Figs. 6.5 and 6.6). 

The intensity of the 
secondary flows is increasingly 
reduced by the kinematic con- 
ditions existing behind the 
impeller outlet and under the 
effect of shearing stresses on 
the diffuser walls. This 
can be seen in Pigs. 8.13 and 
8.14. The vector graphs in 
Pig. 8.15 illustrate the 
relative velocity development 
between planes I and IV for the 
operating point at 14,000 rpm, 
5.31 kg/s, while Fig. 8.16 
again shows the corresponding 


curves in planes I and IV for the operating point near the 
pump limit and absorption limit. We can clearly see how, 
on the way through the diffuser inlet region, both the 
differences in velocity magnitude and in dlrectln for the 
relative velocity vectors largely cancel themselves out. 
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8.2 Comparison between Theoretical and Experimental Results 
Using as a Basis the Jet-Wake Flow Model 

Experimental results on the increase in static pressure 
p/p^Q in the vaneless centrifugal diffuser are compared with 
theoretical pressure curves in Pig. 8.17. The theoretical 
curves here were calculated with expanded Dean-Senoo theory 
for compressible flows (Section 7.8). The symbols in the 
figure stand for mean values of nonsteady pressure transducer 
measuremetns which are defined by Eq. (5.38); and they in 
dicate results which were determined with "conventional” wall 
pressure holes. For the most part, these are mean values 


140 





n/n2 - H/K2 * l»O0y 



Vorderwand = front wal.l 


Plg. 8.15. Development of the relative velocity 
curve W 0 /U 2 in the diffuser inlet region, planes 

,JS~ 

Operating point: " 14,000 rpm, 

in /5J/60 " 

based on several measurements on the same radius — on the front 
and rear walls the diffuser, l.e. up to seven circumferential 
positions. 

The figure legend — lower left in Fig. 8.17 — expresses 
in con<^e fo^n the parameter variation in the calculated sample 
cases (D - © . It also applies for the following figures up 
to and including Pig. 8.23. The meanings of these symbols are 
given below: 

S/T + : Disturbed Jet-wake flow; the input data for the 

theory were determined experimentally in 
measuring plane I (Rg/R2 “ 1.017); 
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Vorderwand = front wall 

Fig. 8 . 16 . Comparison of relative velocity 
curves v ^/U 2 in planes I and IV, scale, see 

Operating point at = 14,000 rpm, 

• 0 '*-53 kg/s, © 6.07 kg/s. 

: axially symmetrical substitute flow with the 
same mass flow and inlet torque as the momentum- 
averaged S/T flow. 
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+ : Wall friction coefficient (R) corresponding 

to the broken-line curve in Fig. 7-3, taken 
into account; in additloiij^for the mass 
6.07 kg/s a sample case Clj with c^ = O.OO 65 = 
constant was calculated; 


- : without the effect of wall friction. 

C + : Interface friction coefficient C = 0.25; 

- : wltnout Interface friction. 
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Pig. 8.17. Comparison of theoretical and 
measured values: curve of static pressure 
P/PtO 1ft. the centrifugal diffuser. Remarks: 
curve (S is not shown here, it coincides 
with © ; curve (^) (6.07 kg/s) corresponds 
to @ , only the Cf curve differs. 

Key: A) Measurement C) Conventional 
B) Nonsteady 


p(R) + : Compressible flow; 

: Incon^-ressible flow, p * constant = p . 

^e 

+ : Boundary layer displacement thickness corres- 

ponding to the solid-line curve in Fig. 7.4, 
taken into account; 

: without boundary layer Influence. 

This yielded the f -..i wing conditions for the theoretically 
calculated sample cases: 
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(p - disturbed, friction-affected, compressible (like ; 

^ - undisturbed, friction-affected, compressible (inter- 
face friction not applicable); 

disturbed, friction-affected, incompressible; 

disturbed, frictionless, compressible (boundary layer 
not a nlicable); 

0 - like 0, but without boundary layer displacement; 

0 - like 0 , but without interface friction. 

As Fig. 8.17 confirms, the measured static pressure curve 
is very well approximated by curve (^ in all three operating 
points. It is striking that the measured values, in general, 
lie al^ve the solution for the axially symmetrical substitute 
flow 0 . This results from the larger inlet momentum of the 
disturbed flow. As the mass flow Increases, the irregularity 
of the velocity distribution in the diffuser inlet increases 
(cf., for example. Fig. 6.8 and the £ 1 ^ 2 ( 0 ) values in Table 6.1), 
and, accordingly, the differences between (0 and (0 Increase. 

The difference between curves 0 and 0 reflects the 
effect of compressibility at quiescent mach numbers of the ab- 
solute ^loclty in the range “ 0.?3 -O. 78 . Comparison of 

curve ( 2 ) with curve 0 illustrates the effect of wall 
friction and interface friction^, as well as the effect of the / 18 I 
boundary layer displacement thickness on the build up of static 
pressure in the d^fuser. In addition, for the mass flow 5.31 
5.31 kg/s, curve gives the individual effect of the displace- 
ment thickness the case of friction-affected flow. 

Finally, curve for 6.07 kg/s makes it possible to estimate 

the influence of modified wall friction coefficient — Cf = O.OO 65 . 

In summary, the following can be inferred from Fig. 8.17: 

In order to make an exact determination of the static 
pressure in the vaneless diffuser of a highly loaded centrifugal 
compressor stage the following factors — in order of Importance — 
must be taken into consideration: 

- the compressibility, 

- the disturbed Jet-wake flow and 

1. A graph of curve (0 — effect of c^- was dispensed with 
here; the deviations with respect to lie in the range of 
the accuracy of the graphs. 


0 - 
® - 



- the boundary layer displacement thickness. 

The corresponding total pressure curve Pt/PtO “ f(R/R2) is 
graphed in Fig. 8 . 18 . In the Dean-Senoo Theory the total pressure 



R/R2 



Pig. 8 18. Comparison between theoretical and 
measuiv ’ values: radial curve of the total pres- 
sure raiio p^/PtO* 

Key: A) Measurement C) Conventional 
B) Nonsteady 


is calculated from the static pressure and the kinetic flow 
energy. With the two-dimensional jet-wake flow model — especially 
given the assumption of a uniform relative flow direction 6 — it 
is naturally not possible to fully approximate the real flow. 

This causes the difference between the measured and calculated 
total pressures at Rg/R? ® I.017 (curve ^ ?''ows the results for 
the incompressible flow). In order to m^e it easier tocompare 
the theoretical pressure curves, in addition solution ^ for the 
axially symmetrical substitute flow — which indeed has a smaller 
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inlet momentum than the disturbed flow — was referred to the 
theoretical initial total pressure of the disturbed, compressible 
flow (e.g. shift for (2) in the case of 5* 31 kg/s of Pt/PtO ~ 

+ 0 . 012 ). 

In curve for the disturbed, friction-affected P;pw the 
sharper drop in total pressure dp^/dR with respect to in /l82 
the range of the disturbed flow up to the mixing radii ratio 
R/R 2 = 1.25 -I .3 (of. Fig. 8.20) is clearly evident. This is 
caused by the higher wall friction losses in the wake area 
(absolute velocity > Cs, see Fig. 8.23) and by the Interface 
friction. Theoretically the total efficiency between measuring 
points ” 0 ” (compressor inlet) and "4" (R/R 2 * I. 687 ) is t^reby 
reduced with respect to the undisturbed substitute flow (2) only 
by an amount of 3 = -0.008, -O.OO 8 and 0.01 for the mass flows 
4 . 53 , 5.31 and 6 . 07 * kg/s (cf. Fig. 8.25). The influence of the 
interface friction a^one, as^xpressed in Fig. 8. 18 by the 
difference between (p and (S' (5*31 kg/s), amounts to Ant,s “ 

-0.002, hence an insignificant order of magnitude. These small 
loss components are not surprising in view of the relatively 
small relative kinetic energy which comes to bear in the mixing 
of the disturbed flow in the sample cases presented. In addition, 
however, mention must also be made here of the short comings of /183 

the simple theoretical model which will be discussed in summary 
at the end of this section. 

Basically, curves ® . © and reveal the dominating 
effect of wall friction on the diffuser losses. In the theoretical 
calculations presented here, Cf«(R) takes into account the actual 
wall friction. Usually, in a one-dimensional analysis of dif- 
fuser flow one introduces a dissipation coefficient c^j, which 
then in the form of an apparent wall friction also approximately 
allows for the dissipation of energy due to turbulent mixing and 
Irregular velocity distributions. With the Cf curve used here, 
the Dean-Senoo Theory Ignores the turbulent energy dissipation. 

If as a rexative measure of the turbulence at the diffuser inlet 
one uses the total pressure fluctuation distributions from Fig. 6.4, 
then for the operating point at 14,000 rpm one sees a distinct 
increase in the average fluctuation Intensity with increasing , 
mass flow. Analogously to this, in Fig. 8. 18 one observes in- 
creasing deviations^etween the measured total curves and the 
theoretical curve ^ . Therefore it is obvious to suspect the 
turbulent energy dissipation, which is not taken into account, as 
the essential cause of these deviations. For a theoretical esti- 
mate of this effect, for the mass flow of 6. 07 kg/s, the case of 
^ was calculated with an "apparent" wall friction coefficient 
a Cf » O.QO 65 ““ but with otherwise identical conditions with 
respect to (p . This curve closely approximates the measured 
drop in tot^ pressure in the diffuser with the exception of the 
immediate inlet region. 
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Traupel [6.^] suggests - 1.5-2 Cf. With this, the 
experimental value is well confirmed. It was not attempted, 
however, to further adapt the dissipation coefficient to the 
measured values. The theory used Indeed offers precisely the 
possibility of analyzing a few individual effects more precisely, 
which otherwise are considered as a lump sum in the Cj values. 

So far we do not have a reliable correlation between the dis- 
sipation coefficient and the flow turbulence in the diffuser. 

The results shown in Pigs. 6.^, 8.3* 8.4, 8. 7-8.9 and 8.18 are 
indicative values for this. 

In Pig. 8.19 the experimentally and theoretically determined /184 
angular momentum curves in the disturbed and mixed flow region 
of the centrifugal dlffusgr are compared for the three operating 
points selected. Curve again represents the solution for 

the axially symmetrical substitute flow — with the same mass 
flow and inlet angular momentum as the disturbed flow. The other 
curves refer to Jet-wake velocity distributions at the diffuser 
inlet. In the cases shown, the angular momentum of the wake 
region clearly exceeds the values for the jet region. It is 
striking that in the inlet radii ratio Rg/R2 “ 1.017 the initial 
value of the theoretical curves for the disturbed flow do not 
exactly argee with the measured values of the Jet region and 
wake region. Again, the reason for this is the assumption made 
in the model, for purposes of simplification, of a uniform relative 
flow direction 6. In actuality, differences in direction of up 
to 10® between Jet areas and wake areas occur in the real flow 
(cf. Pig. 8.22). 

First of all, if we consider the frictionless solution 6) /185 

for the compressible, disturbed flow, it turns out that after a 
while the differences in angular momentum, even without the 
effect of friction, cancel themselves out and an undisturbed 
flow becomes established. This results from the rise in static 
pressure in the centrifugal diffuser and, coupled with this, the 
reversible exchange of energy between Jet and wake regions (cf. 

Section 7.1.1) . 

First of all, we observe sn Increase in the angular momentum 
of the wake in the diffuser inlet region, l.e. the Jet region 
transfers energy to the wake and thus also increases the dif- 
ferences in the Jet-wake energy level. Then the direction of the 
energy transfer reverses — the wake performs work on the Jet 
region and thus contributes to the homogenization of the flow. 

Pig. 8.19 confirms that the radial distance of the zone of 
> Increasing angular momentum differences likewise increases with 

increasing majs fiow. 

I 

I Under the effect of wall friction and interface friction, 

i the mixing of the disturbed flow is compj^eted considerable more 

quioKly than for the frictionless case @ . The basic mechanism 
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1.0 1.2 1.4 1^V) V U 1.61.0 i; 1.4 1.6 

R/Rj—^ 


Fig. 8.19. Comparison between theoretical and 
measured values : radial angular momentum development 
for Jet, wake and total flow. 

Key: A) Measuement D) Wake 

B) Nonsteady E) Conventional 

C) Jet 


of the energy exchange, however, can be perceived in all of the 
theoretical solutions for the disturbed flow^For themass flow 
of 5.31 kg/s, the difference between curves and (o) Illustrates 
the considerable effect which interface friction alone has on the 
homogenization of flow. At the same time, it turns out that the 
selection of an Interface friction coefficient on the order of 
C = 0.?5 is necessary in order to approximately predict the drop 
in the angular momentum of the wake which was found experimentally. 

Curve ( 5 } was calculated without boundary layer obstruction. 

It revealed That the decrease in the diffuser area ratio (outlet/ 
inlet) causes an acceleration in the mixing process. The...-.^ 
difference in wall friction coefficients between curves 
(cf » 0.0065) and @ has, on the other hand no decisive perceivable 
effects in the Jet-wake angular momentum curve — as is obvious 
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for the mass flow of 6.07 kg/s. The effect of the higher Cf 

value first makes itself felt in the region of undisturbed flow /186 

in the form of a sharper drop in angular momentum. 

We did not graph the incompressible curve . It nearly 
coincides with curve (^. 

The comparison of experimental values with theoretical curve 
(D shows that the theory indeed to a large degree accurately 
reproduces the tendency of the angular momentum curve for the 
disturbed flow. In places, however, there are sometimes larger 
deviations. In particular, it is striking that the theory — 
because of the restrictive handicap of a uniform relative flow 
direction 3 — obviously over estimates the increase in angular 
momentum of the wake in the diffuser inlet. Only for the mass 
flow of 6.07 kg/s do the measured values in the wake region pro- 
duce a curve which suggests such a shape. It remains to be 
noted, of course, that these data, as mean values, represent 
the entire wake region. In Pigs. 8. 3-8. 6, on the basis of the 
nonsteady total pressure measurements, the reversible exchange of 
energy in partial regions of the flow can also clearly be 
detected experimentally. 

In Fig. 8.20 experimental results on the radial development 
of the most important parameter for chara(^rising the Jet-wake 
flow are compared with theoretical curve Q , The graphs show 
the curve of the surface component of the wake e = A<j/A, the 
relative velocity ratio y * W 1 J/W 3 and the interference parameter 
5 = e(l-Y). In general, the experimental values measured in the 
diffuser inlet region are in good agreement with the theoretical 
prediction. A striking feature is that the mixing of the dis- 
turbed flow is accomplished primarily due to the approximation of 
relative velocities, while the surface component of the wake of 
the total flow up to complete mixing is technically reduced only 
by 1 ^- 16 $. As already mentioned in connection wit-.h the angular 
momentum distributions, a distinct relationship between the 
reduction in disturbance and mass flow can be seen in the <D curve. 

Also plotted in Fig. 8.20, as a result of the theoretical 
calculations was the radial curve of the fluctuation in pressure 
q related to the total pressure of the compressor inlet Pto* 
definition (cf. discussion in text for Fig. 7.1b), q is positive /I 87 
if the static pressure at the front of the Jet (in the direction 
of rotation) is greater than at the rear of the Jet, and thus 
at the same time the Jet transfers energy to the wake. Thus, 
the pressure fluctuation q is a measure of the Intensity and 
direction of the transfer of energy between the two flow regions. 

For example, in the radial q/pto curve for the mass flow of 
5.31 kg/s we see a change in sign to negative values at R/R 2 ■ 

1.12. If now, as a basis of comparison, we refer to the nonsteady 









Pig. 8 . 20 . Comgarlson between theoretical 
I values (curve (p ) and measured values: 

curves of the surface component of the wake 
I e, the relative velocity ratio y, the inter- 

* ference parameter 6 and the pressure fluctuation 

'• Q/PtO diffuser inlet region 

I Key: A) rpm 

i B) Prom nonsteady measurement 

[ C) Jet-wake mixing (theory) 

i 

i 

I measured static pressure curves — Pig. 8.2 ~ than in the two 

' measuring planes immediately behind the impeller outlet (R/R2 ■ 

I 1.017 and 1.039) we see pressure curves which basically resemble 

i the curve of values measured in the impeller — pressure 

i maximum near the pressure side DS of the flow channel and minlmun 

; in the vicinity of the suction side. Prom this it follows that 

I ~ as a result of viscosity forces — the transfer of energy in 

the impeller affects the static pressure curve behind the impeller 
i up to about R/R2 * 1.05. It is conceivable that these greater /l 8 b 

pressure variations here conceal those pressure pulsations which, 

\ theoretically, must occur as a result of interactions between 

Jet and wake regions. 
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Not until the two following measuring planes (R/R 2 * I.O 89 
and 1 . 151 ) does one observe extreme pressure values also In the 
middle of the channel division, as corresponds to the theoretical 
model. For R/R 2 = I.O 89 , a pressure maximum actually occurs at 
y/t - 0.4 In the front region of the Jet, and progressing towards 
R/Rp = l.l4l a pressure minimum Is found at about the same cir- 
cumferential position; l.e. the change In sign between these 
two measuring planes and also the order of magnitude of the 
static pressure fluctuations measured In Fig. 8.2 correspond 
exactly to the theoretical curve shown In Fig. 8.20. Tho 
same phenomenon was detected by Senoo and Ishlda In their 
measurements [7.5]. 


Thus, supported by the theoretical prediction and Independent 
of the other measurements, this gives us further Indication of the 
reversible exchanr^es of energy within the disturbed flow which 
Is Important for the mixing of the flow. 



/139 


R/f^* 


Fig. 8.21. Comparison between theory and experiment: 
radial curve of the absolute flow direction a In the 
Jet region and wake region. 

Key: A; rpm C) Jet 

B) Nonsveady D) Wake 

measurement 
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Pig. 8.22. Comparison between theory and 
experiment: radial curve of the relative flow 
direction 3 for Jet, wake and total flow. 

Key: A) rpm C) Jet 

B) Nonsteady D) Wake 

measurement E) Theory 


Fig, 8.21 gives information about the mean values of the /188 

absolute flow direction measured in the diffuser inlet region 
or the for the Jet and wake, and it also gives the corresponding 
thgoretlcal curve with and without the effects of friction 
(© and ®)^as well the curve without the effect of interface 
friction C ( (6/ ) . As a result of theoretical estimates based on 
the assumption of an insignificant wake mass flow. Dean and Senoo 
[7.1] and later Hess [8.1] reached the conclusion that the dis- 
turbed flow up to R/R2 » 1.1 have largely mlxea. Pig. 8.21 
confirms that this finding cannot be generalized, and at least 
given the existing, relalstic facts in this cas<. (blade angle 
3«2 • 90®, finite wake mass flow) at R/R 2 ■ 1.3 considerable 
fluctuations in direction still occur. Certianly these would 
also strongly affect the operational behavior of an incipient, 
bladed diffuser in this case. 

As a final exang>le of this series of comparisons between 
thf »ry and experiment, in Pig. 8.22 the relative flow directions /190 

g measured in the Jet and wake are compared with the curve of the 
momentum mean value 3 of theoretical curve . A striking 




feature of the measured values for 5*31 and 6.07 kg/s Is that 
3 at the beginning of the waxe at first exceeds by about 10® the 
angle of the Jet region. But then the radial angle curve in the 
wake drops off more sharply and at R/R 2 * 1.151 the wake values 
are clearly less than the Jet values. The theory ignores this 
effect — the posslMe consequences of this limitation have 
Indeed already been referred to several times above — however it 
quite closely approximares the prevailing Jet curve for the mean 
momentum value. Pig. 8.15 shows the result. e of e, parameter 
study (cf. Section 7.3) which allows us to estimate the upper 
limit of possible losses due to variable Jet-wake relative flow 
directions. 



Pig. 8 . 23 . Exauflple for t^ 
theoretical calculation (p of 
the disturbed Jet-wake v^ocity 
curve in the centrifugal dif- 
fuser. 


developed by Dean and Senoo also 
cooipressor stages. 


Plnally, as a product of 
the friction-affected, compres- 
sible theory (p , Pig. 8.23 
shows the radial curve of all 
velocity components and f]ow 
directions in the absolute and 
relative system for the Jet, 
wake and the mixed flow (also 
designated with the subscript 
S « Jet). This figure Illustrates 
in the concise form of a typical 
example the development of the 
flow field during mixing of the 
Jet-wake distribution in the 
vaneless diffuser of a highly- 
loadeci centrifugal cor.pressor 
stage. 

If we take into account 
the measurement uncertalnt ' js 
(see Table 5.1) tnen In con- 
clusion we can say that the 
comparison between experimental 
findings and the expanded theory 
for compressible flows carried 
out in Pigs. 8.17-8.22 produced 
a good agreement. In particular, 
this is true for the prediction 
of the characteristic inter- 
ference variable e, y and 6, 
the rise in static pressure and 
the momentum curve in the 
disturbed diffuser flow. Prom 
this standpoint, the present 
Inve tigations yield a con- 
firmation of the sl.nple two- 
dimensional Jet-waice flow model 
for highly loaded centrifugal 




With certain qualifications, this is true for the prediction 
of diffuser losses. Essentially the theory accurately predicts 
the increased drop jn total pressure in the region of the 
disturbed flow. Moreover, the possibility of an Improved loss 
estimate was illustrated by introducing a turbulence-dependent 
dissipation coefficient. Larger deviations — up to 1 . 5 % with 
respect to the measured total pressure ratio — are seen 
immediately behind the Impeller outlet, 1.0 < R/R 2 < I. 05 . The 
theory ignores here, in particular, the momentum loss due to 
the rapid evening out of the absolute velocity distribution 
as is revealed, for example, in Fig. 8.3 in the development of 
total pressure between R/R 2 = I.OI 7 and 1.039* In addition, 
differences result from the simplified two-dimensional jet- 
wake model. Taking as a basis the actual form of the interface — 
as it appears in Figs. 8.14 and 8.1 — an Increase in the jet- 
wake interface friction to two to three times the value of the 
simple theory can be predicted. An improved loss estimate could 
be obtained by a refinement of the model in which the total flow 
would be divided into several bands of constant, relative 
diffuser width Az/b in order to better approximate the three- 
dimensional character of the jet-wake distribution in the dif- 
fuser inlet region. 

8.3 Theoretical Parameter Study on the Jet-Wake Mixing 

Process 


With the aid of the jet-wake theory, it is attempted in 
Fig. 8.24 to give indications of the variable mixing rates behind 
centrifugal Impellers with backward bent and radially ending 
blades. Assuming that the diffuser inlet radius R 2 , or diffuser 
inlet width b 2 and Impeller blade number z correspond to the 
dimensions of the stage being studied here, the realistic 
assumptions (specified in the upper right hand corner of 
Fig. 8.24) about a' "normal disturbance" are made. Here the 
angular momentum parameter 2 stands for the ratio of the 
circumferential compongnt £0 the radial component of the ab- 
solute velocity: X 2 = c^ 2 /cr 2 • addition, for the impeller 
with a blade angle of Bs 2 = 45°, the mean jet-wake relative 
flow direction 83 was estimated at 42.5°, and correspondingly 
for the 90 ° impeller, B 2 = 75° • 

The legend in the upper left hand corner of tne figure refers 
to the calculated cases. The solid-line curve is for the equal- 
area diffuser (R*b = constant). Moreover, the same assumption: 
are made for it as for the theoretical case 1 in Section 8.2 
above. If we compare these results for both Impellers, then 
it clearly follows from the curve of the disturbance parameter 
6 = e ( 1 -y) and the curve of the relative velocity ratio 
Y * Wiji/Ws that the disturbance behind the Impeller with back- 
ward bent blades quickly subsides in the radial direction. 
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Between the equal-area and parallel-walled diffuser (b = 
constant) there are only insignificant differences in the 
6 curve. To be sure, in the latter case the disturbance is 
completely eliminated only at infinity. The two remaining 
curves for Cf = f(R), ? = 0 and Cf = 0 , ? = 0 (calculation 
without boundary layer displacement) moreover clearly confirm 
the importance of the Jet-wake interface friction 5 for a 
rapid mixing of the flow. 


To illustrate the efficiency losses through the Jet-wake 
diffuser flow. Fig. 8.25 shows by what amount Aht s the total 
efficinecy of the stage decreases in the case of disturbed 
flow in comparison with an axially symmetric substitute flow. 
The total efficiency is defined as ni ^04 Eq. (4.3). This 
theoretical estimate is based on the same conditions as the 
sample cases and at 14,000 rpm in Section 8.2. As a 
function of surface component £2 of th§ wake, the fol- 
lowing equivalencies were setup: X 2 ~ 3 and 5 ; 

and Y 2 = Wip/Ws = 0.2 and 0.6. From this one can assume that 
the greatest portion of the Jet-wake disturbances present, for 
all practical purposes, behind highly-loaded centrifugal im- 



R/R:— • 


Fig. 8.24. Comparison of the 
radial development of a Jet- 
wake disturbance behind a cen- 
trifugal impeller with backward 
bent and radially ending blades: 
6 = e(1-y), y = W^/Vg. 


pellers lie within these limits. 

Thus, it .turns out that the 
results for the 3 operating 
points at 14,000 rpm can well 
be incorporated into the 
picture of this parameter 
study. !'■ appears relalstic to 
assume a value of 3 - 
- 0.02 as the upper limit for the 
Jet-wake mixing losses in the 
optimum operating range of a 
compressor stage. Of course, 
it should be noted that with 
pronounced Jet-wake velocity 
differences, Y2 - 0.02, and 
relatively large mass flows, 

Yp = 1, the losses can Increase 
sharply. Values in this order /194 

of magnitude result when the 
absorption limit of the com- 
pressor performance graph is 
approached. The results shown 
here suggest, then, that 4 to 8 
efficiency points of the sharp 
drop in efficiency frequently 
observed in this operational 
range be attributed to the 
mixing losses. 
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Within the scope of this loss analysis we will also discuss 
the results of the theoretical estimation treated in Section 
7.2. By means of an expanded model it becomes possible to 
localize the hertofore ignored effect of variable Jet-wake 
relative flow directions on the loss development. In Fig. 8.26 
the measured disturbance data BeS» ^e were given for 

the operating point at 1*4,000 rpm, 5.31 kg/s, in two radii 
ratios Re/Rp “ 1.017 and 1.151, and in addition the wake relative 
flow direction Bgip was systematically varied between 0° and 90®. 


Pig. 8.25 Fig. 8.26 




Fig. 8 . 25 . Study on the effect of the jet- 
wake disturbance variables C 2 and y 2 » well 
as the angular momentum parameter A 2 on the 
total efficiency of the stage t \^. . 

U } o 

Fig. 8 . 26 . Effect of va.’iable jet-wake relative 
flow directions 6 on the angular momentum R‘Cu 
and the kinetic energy flow E in the mixing 
radius R, values in () followed by an * are for 
undisturbed flow = Bgs* 


Fig. 8.22 had already shown that at R/R 2 “ 1.017 the wake angle 

BgT “ 80 . 5 ® is greater than the corresponding value in the jet 

region, but at R/R 2 “ 1.151 with BgT * 34. ^® it is clearly 

smaller than the corresponding value in the wake region. Under 

the fonditipns mentioned, it can be deduced from the curve of 

(r2*e)/(r 2*E)» in Fig. 8.26, l.e. the ratio of the kinetic 

energy flow E in the mixing radius R at (B©'? BeS^/^^eT “ ®eS^» /195 

that for the sample case selected these additional losses 

amount to less than 1 % of the kinetic flow energy. Altogether, 
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these results show that the consideration of the effect of 
variable Jet-wake relative flow directions would have changed 
only slightly the loss predictions of the Dean-Senoo theory. 


Finally, Fig. 8.27, for the geometry of the centrifugal com- 
pressor stage being studied here and for the theoretical as- 
sumptions corresponding to case ^ (cf. Pig. 8.17) shows the 
effect which the angular momentum parameter X 2 and the dis- 
turbance variables C 2 and Y2 have on the radial distance of 
the mixing radius R* from the Impeller outlet (related to blade 
division t 2 ). Prom a practical standpoint, the flow for y* > 0*9 
is regarded as "equalized.” 



Fig. 8 . 27 . Study on the effect 
of the angular momentum para- 
meter X2» 0^* the Jet-wake 

disturbance variables c2 Y2 
on the Jet-wake mixing radius 
R». 


A surprising feature of 
these findings is the small 
Influence of the relative 
velocity ratio Y 2 . With an 
approximately constant inlet 
angular momentum, the mixing 
radius ratio R*/R 2 primarily 
depends on the mass flow, the /I 96 
Impeller blade number and the 
surface component £2 of the 
wake . 

Additional findings of 
more extensive parameter 
studies on Jet-wake mixing for 
Impellers with backward bent 
blades arc found in Dean and 
Senoo [ 7 . 13 , and for .90®- 
Impellers in Hess [8.1] and 
Johnston and Dean [7.2]. Of 
course, it should be noted 
that the assumption made in 
[ 7 . 2 ] and [8.1] of £3 = 0 
is not very realistic in view 
of the experimental results 
presented here. 


9 . Summary 


With new measuring and evaluation techniques with high time 
resolution the development of disturbed, compressible flow be- 
hind a highly loaded centrifugal compressor Impeller was exactly 
determined and theoretically analyzed in comparison with cal- 
culations based on a two-dimensional Jet-wake flow model. The 
results enable basic conclusions to be drawn on the complex flow 
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pattern and the loss distribution in high-performance centrigual 
compressors. Thus they yield Important pre-conditions for the 
development of an Improved stage design. 

New measuring techniques had to be developed for the experi- 
mental analysis of the complex flow discharged from the Impeller. 

In so doing, we succeeded in developing a measuring system which 
allows the disturbed flow pattern to be determined with quanti- 
tative exactness under very difficult measurement conditions.^ 

The following points should be emphasized; 

- the development of a measuring technique for determ^’ning 
the nonsteady direction field behind the impeller outlet 
which, for the flow analysis, makes possible the important 
transition into the relative system; 

- the use of the synchronized sampling technique as an im- 
portant prerequisite for the exact, digital processing 
of a large number of nonsteady test signal curves from a 
strongly fluctuating flow; 

- the study on the operational behavior of miniature semi- 
conductor transducers, as a precondition for the development 

of exact methods for nonsteady wall pressure and total /I 98 

pressure measurements. 

With this measuring technique the disturbed flow behind a 
centrifugal compressor Impeller with radially ending blades was 
studied in detail in several operating points at rotations speeds 
between 10,000 and 18,000 rpm and mass flows between 3*8 and 7*^ 
kg/s. Results of measurements made in four measurement planes 
staggered radially in the inlet region of the vaneless diffuser 
were presented. These include: 

- the nonsteady wall pressure curves, 

- the distributions of total pressure and total pressure 
fluctuation intensity, 


1. As a by-product of these measuring technique investigations, 
a few findings resulted which did not find any direct exj'”esslon 
in this study. They may, however, be of Interest, for example 
for industrial development projects. These include, for example, 
estimates on the error reading of conventional pressure and 
temperature measuring methods in the nonsteady flow region or 
simple measuring systems which make possible a quick, qualitative 
analysis of the disturbed flow behind the impeller outlet. These 
results were Included in Report C9«l]« 



- the distributions of the absolute flow direction, and 
derived from this 

- the distributions of the radial component of the absolute 
velocity vector, and 

- the distributions of the relative velocity vector. 

Immediately behind the impeller output the relative flow 
possesses a clearly separate Jet-wake distribution. This 
becomes more pronounced with increasing rotation jpeed and mass 
flow. The development of the disturbed flow in the Impeller 
Made channel — starting from a separation of flow in the 
region of the meridional bend in the channel — is illustrated 
and discussed with the aid of sample measurements made with a 
laser technique for measuring velocity. The wake forms in the 
region of the front wall of the channel near the blade suction 
side. For the operating points studied and with respect to 
the total flow, it possesses a surface component of 35-4135 
at the impeller outlet and a mass flow component of 19-265#. The 
ratio of the Jet-wake relative velocities ranges between 0.46- 
0.55. This small variation in disturbance parameters underlines 
t‘he fact that the relative importance of individual flow 
mechanisms, which contribute to the creation and development 
of the region of separation in the Impeller, is largely 
maintained in a large operational range Characteristically, /199 
the distribution of the relative total ,.ressure behind the 
Impeller outlet show a high loss concentration in the wake 
region, while the Jet region possesses an Isentroplc flow 
core. 

The radial development of the measured velocity distributions 
clearly illustrates the equalization and mixing process within 
the disturbed flow. With the condition studied here, however, 
the disturbance — above all in the region of the front wall 
of the diffuser — is not yet completely eliminated up to the 
radii ratio R/R 2 =1.15. By comparing the total pressure 
distribution in the various measurement radii behind the im- 
peller outlet, the reversible excha .. :e of energy between areas 
with high and low kinetic energy in the relative system was 
detected. Results of a theoretical study on the effect of 
Coriolis accelerations and normal accelerations on the turbulent 
mixing intensity at the Jet-wake Interfaces are in good agreement 
with the experimentally determined curves. 

The development of the disturbed flow in the diffuser inlet 
region was exactly analyzed in a comparison of experimental 
results with calculations based on the two-dimensional Jet-wake 
flow model of Dean and Senoo. In contrast to the original 
theory, the considerable effects with high-performance cen- 
trifugal compressors — flow compressibility and boundary layer 
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displacement — on the conversion of energy in the vaneless 
diffuser were also taken into consideration. V/itn this expanded 
theory, calculated radial curves of the characteristic distur- 
bance parameters and of the static pressure were In good agreement 
with measured values. By contrast, by the assumption of an 
axially symmetric flow the prediction for the rise In static 
pressure would be too low. The angular momentum distribution of 
the disturbed flow illustrates the effect of the reversible 
exchange of energy and of the jet-wake interface friction on 
the homogenization of the flow. 

The radial total pressure development is clearly influenced 
by wall friction effects. Immediately behind the impeller out- 
let the disturbed flow leads to increased total pressure losses. 

In the operating points studied here, the total efficiency of /200 
the stage decreases by less than 2 points under the effect of 
the disturbed flow. A parameter study shows that these losses 
— above all in the vicinity of the performance graph absorption 
limit — may climb to ^-8 efficiency points. 

Finally, it must be said that the results of the detailed 
measurements behind a highly-loaded centrifugal compressor 
impeller are closely approximated by the expanded Jet-wake 
theory. At the present time, the application of this theory 
is still affected by the disadvantage that the starting con- 
ditions of the disturbed flow can be exactly determined only 
by constly measurements (in case one is not satisfied with 
estimates or results from simple measuring techniques [9*1]). 

In the meantime, as a result of the introduction of the 
dual-focus laser method for measuring velocity, it has been 
possible also to exactly measure the flow pattern within 
cnetrlfugal impellers. Theoretical flow models derived from 
such measurements should be capable of predetermining the 
condition of flow at the impeller outlet. Thus appears the 
possibility of simulating the complex flow in high-performance 
centrifugal compressors in a simplified model system based on 
the physics of fluid flow, and from this derive fairly optimal 
stage configurations. 

Dietrich Eckardt, Certified Engineer 
5,000 Cologne 90, Linder Ht5he 
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I. Investigation of Temperature Effects on Miniature 
Semiconductor Pressure Transducers 


Experimental Setup and Measurement Program 

For these investigations and for calibrating the 
total temperature probe a thermostat was designed and built 
which meets particularly nigh demands for measurement accuracy 
and control accuracy of the experimental temperature. Fig. I.l 
shows the main construction of the apparatus. B; separating the 
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Fig. I.l. Experimental setup for investigating temperature effects 
on miniature semiconductor pressure transducers. A - Thermostat 
adapter for pressure calubration at variable ambient temperature; 

B - adapter for Investigating the temperature gradient effect; 

1- propeller; 2 - straightener ; 3 - hot/cold control; 4 - circu- 
lating pump; J - insulation; 6 - condensor; DA - pressure trans- 
ducer; NTC - thermistor, Pt - resistance thermometer; p - test 
Px-’essure; - atmospheric pressure. 







measuring bath and the heating bath and building in a stiring 
propeller (1) and straightener (2) in the measuring bath, and 
buidling in a circulating pump (4) in the heating bath, very 
uniform temperature fields are produced within the test region 
. A special bath fluid^makes it possible to operate the 
^ermostat at calubration temperatures between -20° c and +250°C. 
The sensitivity of response of the hot/cold control (3) is 
around 0.5K and the measurement accuracy of the calibrated Pt 
resistance thermometer with a 5-place digital ohmmeter is 
±0.1K. To study the effect of changing ambient temperatures 
on the zero point drift and the sensitivity of the semiconductor 
pressure transducers the adapter (A) shown in Fig. I.l (upper- 
right) was inserted into the thermostat (modified for pt probes). 
Thus the transducer can be acted upon with the pressure (p) 
at various ambient temperatures, while at the same time the 
temperature in the immediate vicinity of the pressure transducer 
membrane is measured by means of a NTC thermistor. 

Prior to the actual pressure calibration, the uncooled total 
pressure transducers were carefully exposed to alternating 
temperatures several times (Including approximately 2 hours at 
+150°C) and thus the effect of artificial aging was achieved. 

After complete cooling of the bath, pressure calibration 
curves were then plotted at each of 4 ambient temperatures (20°C, 
60°C, 100°C, 140°C) for all the total pressure probes, Kulite 
XBL-9-093-25 G (maximum operating temperature, 17b°C). In 
addition, four Kulite XQL-093-25 wall pressure transducers (max- 
imum operating temperature 120°C), which were installed in the 
centrifugal compressor in water-cooled adapters ^ were calibrated 
for a general check of the operating behavior of semiconductor 
pressure transducer — 20°C, 60°C and 100°C. 

Deviating from the U.S. test conditions (cf. [3*6, 3*7]), 
which when making temperature changes suggest gradients of 
approximately lOK per hour, a difference of fourty K between 
the test temperatures Was bridged in about thirty minutes.^ 

But then the test temperature, as specified, is neld constant 

for 45 minutes before the pressure calibration is made. During 
the calibration, the readings of the Pt resistance thermometer 
and the NTC thermistor differ by less than 0.2 K. 

For investigating temperature gradient effects on the 
zero point reading of semiconductor pressure transducers adaptor 
(B) shown in Fig. I.l (lower right) was used. It could alter- 


1. Marlotherm L, Chemishe Werke HUls. 

2. Based on empirical values associated wit.; the change in 
operating point of the centrifugal comnressor where the steady 
state is reached in less than 30 minutes. 
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nately be connected to the water circuit of two thermostats. After 
the pressure transducer (DA) inserted in the Cu cylinder had 
be». n held constant at 295 K for 45 minutes, the switch was made 
to the second thermostat with a higher water temperature. During 
this process, the curve of the pressure transducer output voltage 
and the temperature on the measuring membrane were continuously 
recorded (cf. Fig. 3*7)* 

II. The Synchronized Sampling Technique /212 


Principles and Basic Equipment Design 

In the signal processing of rapidly changing, periodically 
recurrent and very noisy test signals by means of the synchro- 
nized sampling technique (SST), two aspects of the technique 
are of prime importance: 

1. The test signal is exp<anded over time (expansion 
factor up to 10^^) and the Fourier spectrum of the 
signal is transformed downward to frequencies which 
lie in the operating range of a normal x-y recorder. 

2. The average, periodically recurrent signal component 
is separated from stochastically distributed signal 
fluctuations. This considerably improves the 
signal-to-nolse ratio. 



time ► 

Fig. II. 1. Principle of the expansion over time of a 
periodic test signal with the aid of the simplified 
sampling technique, expansion factor a = 6. 
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Expansion Over Time of the Periodic Test Signal 

This process Is explained by referring to Fig, II. 1. If one 
scans the Instantaneous voltage of the test signal In the time 
Intervals t = tj. + At, whereby At Is steadily Increased, this 
produces the output signal represented by the circles. If the 
scanning points are sufficiently dense, then the SST output 
signal gives an Image of the test signal expanded over time. 

The expansion factor Is a = tj,/At and can be made any size with 
At -*• 0 . 

Thus the synchronli: d sampling technique Is an electronic 
"Ar.alogon" to the familiar optical stroboscope technique In 
which a periodically recurrent pattern Is "apparently" thereby 
slowed down so that it Is not observed continuously but only 
momentarily In certain time Intervals, whose reference is fixed 
to the repitition rate of the original process. 



Fig. II. 2. Block diagram of the "sample/hold/ 
integration" unit of the SST apparatus. 


To accomplish the pointwise scanning requires a selection 
logic and a rapid switch (cf. Fig. II. 2). In order to achieve 
high resolution the opening time of the switch tg must be 
small with respect to the period duration cf the test signal T. 

So that the scanned voltage between the measurement periods holds 
good. It is stored on a condensor Cg. Tge following low~pass 
filter evens the Jumps between the individual scanning points an 
and Is used for averaging. 

Taking the Mean 

Fig. II. 3 will help explain the SST averaging principle. 

At position a It shows a n vjiiot ccidy signal voltage curve with 
a period of T = t^.. 
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i- i. JL i- i. 


tlTrigotr 


(§) 



Fig. II. 3. Diagram of the 
sampling process: 

^ - periodic input signal 
(T: phase relationship of the 
samples); (£) - sampler output 
(pulse height corresponds to 
the instantaneous signal value 
at t); 0 - output signal of 

the voltage stabilizer. 

Key: A) Test signal synchro- 
nization trigger. 

corresponding to the mean value 


The beginning of the period 
is marked by trigger impulses 
which are picked up by a trigger 
disc with the aid of a feromag- 
netic magnetoresistor trans- 
ducer. This turns around with 
the compressor shaft. The 
scanning takes place in each case 
at point T after the trigger 
pulse, i.e. in this example the 
measurement is made without any 
scanning shift (At = 0). This 
restriction simplifies the 
description and produces prac- 
tically the same result as a 
treatment with 0 ^ At << t^. 

At each signal scanning the 
condensor Cg is charged with the 
instantaneous signal voltage, 
position . It holds this 
voltage un€ll the next scanning 
and what arises is the so-called 
box car train signal, position 

This signal train consists /215 
of the d.c. comr ■'nent U, which 
gives the mean value, and an 
a.c. component u, which originates 
from stochastically distributed 
disturbances. The Interferrlng 
a.c. voltage is suppressed in 
the low-pas'j filter after C.g, 
so that the d.c. component 
stands out more sharply. 


The effective value of the fluctuation intensity of the 
test signal is formed at the a.c. voltage output of the 
voltage \lu^ stabilizer. 


In the language of statistics, the sampling scanning process 
produces a sequenc"' of scattering individual measurements Uj_. 


In a first approximation the RC term averages over the 
measured values Ui which were recorded within the time constants 
RC (Integration time tj_). The number of measurements is n = 
t ^/t . 


For stochastically distributed signal deviations the slgnal- 
to-nolse ratio is improved by a factor of /n. Therefore, the 
signal-to-nolse ratio can be improved at will if the time con- 


172 


stant of the RC term is made sufficiently large a large ti, 
however, the total measuring time tm must also be correspondingly 
Increased in order to prevent a reduction of the upper limit 
frequencies of the measuring system. Thus the upper limit for 
signal-to-nolse Improvements is set by practicable measuring 
times and the long term stability of the compress-^r operating 
condition. 

Further information on the SST measuring principle can be 
obtained from the quantitative appraisal in Table II. 1. 

Fig. II. 4 shows the principle arrangement of the sampling 
apparatus used in the form of a block diagram. In addition, a 



Fig. II. 4. Block diagram of the Sampling 
Apparatus for averaging periodic test signals. 


cathode ray oscillogi-aph is required, on whose screen the 
signal train to be t:\-nsferred is adjusted by means of the 
synchronous trigger. The SST scanning measuring time for the 
entire screen width (10 cm) is generally tj^ = 2 minutes. The 
position in time of the sample is set by a voltage Uj,gf« while 
tpj steadily Increased. A comparator compares this voltage with 
tne internal sawtooth voltage of the oscillograph which is 
triggered anew each time the beam passes, and when there is 
coincidence on the rising side of the saw tooth it gives off 
an Impulse which controls a momentarily closing switch (gate 
time tg) through which the Instantaneous signal value Is 
connected to the averaging condenser. Since a quasi d.c. voltage 
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TABLE II. 1. QUANTATIVE APPRAISAL ON THE USE OF THE 
SYNCHRONIZED SAMPLING TECHNIQUE 

Problem: 

The nonsteady wall pressure curve for the region of a blade 
spacing is to be transfered to a x-y recorder and digitalized 
in parallel with this. 

Exapple : 

Rotation speed n = 18,000 rpm. Blade number z = 20 

Signal frequency 6 kHz (167 ys/t) 

(Blade spacing t) 


Option: 


Horizontal time sweep on the cathode ray 
oscillograph 20 ps/cm 
Screen width 10 cm 


Option: 


1,2 signal divisions per screen width 
Beam transit time tg = 200 ys 
Sweep frequency 3 kHz (t^, = 0.333 ms) 


SST - Integration time 
SST - Gate time 


t^ = 0.5 s 
= 2 ys 


Number of samples per mean value n=tj^/t^ = 1,500 


Improvement of the signal-to- 
noise ratio by a factor of 


/n = 38.7 (31,8 dB) 


Instrument constant: SST measuring time for screen width 

t = 120 s 
m 

^ Digital signal resolution (tjj^/t^) < 240 values 


Time expansion 


t /t, = 6-10^ 
m s 


174 


n/ ’IV 


... . fS 


I ■ 


is measured at the condenser, the measurement uncertainty, 
under normal conditions even for signals up into the MHz range, 
is less than ±0.555. The accuracy with respect to time — with 
appropriate constancy of rotation speed — is limited only by 
the sawtooth of the oscillograph. By calibration, an accuracy 
of better than 15? can easily be obtained here. 


III. Total Temperature Measurement in the Disturbed Flow behind 
the Impeller Outlet 


Total temperature measurements with slow stagnation temper 



ature probes immediately behind 
the centrifugal compressor im- 
peller showed larger temperature 
increases — up to about 10^ — 
than '•"'re measured downstream 
in the vicinity of the diffuser 
outlet. As an example. Fig. 

Fig. III.l shows total temper- 
ature profiles (with correction 
of an average recovery factor 
Rf) measured in two radii ratios 
over the diffuser width. The 
difference in these two curves 
is especially strong in the 
region of the front wall of the 
channel . 

The reason for the higher 
temperature reading given by 
the probe installed immediately 
behind the Impeller outlet is 
basically the particularly 
pronounced irregularity of the 
relative flow over the channel 
spacing in this reason. 


z/b 


Fig. III.l. Curve of the total 
temperature increase, plotted 
over diffuser width, measured 
in two radii ratios behind 
the Imoeller outlet. Operating 
point : n//5T = 16,000 rpm, 
m/5"/6o *^6.88 kg/s. 


Pfau [III.l] studied the 
effect of various velocity 
distribution on the temperature 
readings of impellers with 
backward bent blades. 

For Impellers with radially 
ending blades, Moffat and Dean 
[III. 2] carried out detailed 
theoretical analyses on the heat 
transfei condition at the sensor 


of a total temperature probe with nonsteady oncoming flow. The 
following, simplified estimate, taking as a basis the "block” 
velocity distribution shown in Fig. III. 2 above, will indicate 
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the order of magnitude of ^,he error in the temperature reading 
as a function of the disturbance parameters genrally used in 
this study, i.-e. the wake component e and the relative velocity 
ratio of the wake and w,j/ws.^ 

For the mass averaged total temperature Increase, i.e. the 
real temperature Increase after the termination of the mixing 
process, the following equations are valid: 






(III. la) 


(Assumption: density p = constant) and 

AT!j. s w u^/ Cp . (III. lb) 

The latter equation ignores components of the temperature 
increase due to Internal disc friction and turbulent exchange 
movements at the impeller outlet. In the following equations 
it is also only valid for the Jet region and wake region: 

*®t.S • “ ' t* / Op . 

*®t,T * “ (« - c, ,, / tg g) / Cp . 

whereby in keeping with Pig. III. 2 it is assumed: 

® = ®S " ®T 

The symbols are defined as follows: 

Aft total temperature Increase T^^ “ ^cO* 

Cp specific heat at constant pressure, 

Cjjj mean meridian velocity, 
e relative wake component of a spacing. 


(III. 2a) 7220 

(III. 2b) 


1. The derived system of equations largely correspone to an 
unpublished system worked out by Conrad- Daimler-Benz A.G. 
1971. 
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y reduced output factor, 
u circumferential velocity, 

6 direction of the relative flow with respect to the 
negative circumferential tangent. 

On the sensor of the total temperature probe the temperature 
increase indicated is Under the simplifying assumption 

of identical heat transfef* coefficients in the Jet region and 
wake region, the following equation applies: 


"t.« • ‘ * "t.T ♦ “ - «t.S 


(III. 3) 


If If we substitute Eqs. (III. 2a, b) in Eq. (III. 3), then 
the following equation follovis with the flow coefficinet 
f = Cjn/u; 


"t,« * / t* S) / Cp 


( 111 . 4 ) 


By using Eq. (III. lb) we obtain the relationship between 
the indicated and mass averaged total temperature increase: 


^^t,a / * (1 • 5 / tg 5) / w , (III. 5) 


whereby the following relationship is used: 




In order to produce a reference for the disturbed velocity 
distribution Eqs. (III. la, b) are set equal and Eqs. (III. 2a, b) 
are used. If we solve for tg 0 we obtain the following: 




If we substitute for Cm s using Eq. (III. 6), then the 
following equation follows: * 
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A 


or with the following notation 

8L * C_ m't 2 


^ “r^-e I WZ 


nil. 7) 


we finally obtain 


, _ - 1 ♦ X 
tg f - ^ 1 - u 


(III. 8 ) 


This substituted into Eq. (III. 5) then gives; 


AT 


t.a 


AT. 


1 ♦ X / W 
1 ♦ X 


(III. 9) 


To introduce the disturbance parameter y - W 1 J 1 /W 3 , 

Eq. (III. 7) is rearranged with Eq. (III. 6 ) and the meridian 



error in the total temperature reading which, for a specific 
reduced output factor y, depends only on the disturbance 
parameters e and y: 

X » ** c)(l * 

[1 - e (1 - y))^ (III. 10) 


As an example of such a parameter study. Fig. III. 2 shows 
the curve of AT^ as a function of the disturbance variable 

l»wip/ws for thre4 wake components e. An evaluation was made in 
each case for a value of the r'^duced output factor y which is 
typical for imepllers with a ^iven blade angle of Bg 2 . It turns 
out that the error in the temperature reading Increase with 
Icreaslng e and l-wip/ws: or with decreasing y. 

The order of magnitude of the error in the reading I'eveals 
that exact ’measurements of total temperatures in the flow 
discharge^’ *rom the Impeller are only possible in areas in 
which tl disturbances have largely subsided. 

Conversely, it is possible, from the comparison of tempera- 
ture measurements at dlffereni" distances behind the Impeller, 
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Fig. III. 2. Parameter study 
r estimating the error in 
reading of total te..ipera- 
‘ are p’'.'»bes with irregular 
oncor'. 2 flow. 


to infer the degree of non- 
uniformity, even if this can 
only be qualitative because of 
the large number of simplifying 
assumptions . 

For example, analysis of 
the findings shown in Fig. III.l 
permits us to state that the 
flow behind the impeller outlet 
in the backward section of the 
channel (z/b > 0.6) is largely 
equalized in the circumferential 
direction, l.e. the circum- 
ferential disturbance of the flow 
field (wake region) lies in the 
region of the front wall of the 
channel and — assuming wip/wg “ 

0 — there it assumes a wake 
component of at most about 50? 
of the spacing. 

In view of the simple 
measuring method, these state- 
metns coincide rem.arkably well 
with the results of the detailed 
flow studies using measuring 
techniques with high time 
resolution. 

IV. Definition of the Relative 
I*otal Pressure Pt,rel 


Starting with the energy 
equation in the relative system, 
o jllowing equation holds for the relative total enthalpy 

,rel 


t,r«l 


h ♦ w^/2 - u^/2 


with h ■ static enthalpy; w ■ relative velocity; 
ferentlal velocity. 

For a relative flow without energy i.ransfer 
of "internal” dlsi rictlon on the housing wall) 
equation is valid a^.ong a flow line; 

N.r.1 * «p • ■'t.ptl • «"»»• 


(IV.l) 

u ■ clrcum- 

(as a result 
the following 

(IV. 2) 
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with ** relative total temperature. 

In particular, it follows from Eqs. (IV.l) and (IV. 2), 
with nontwisting oncoming flow to the impeller (w2/2-u2/2 »• 
0 ^/ 2 ) that: 


^t,rel * ^tO,N 


(IV. 3) 


with T^q jg ■ (standard) intake total temperature. 


t 



Fig. IV.l. T,s diagran f any pi^sure change 
between the total states @ and (^ . 


In the T,s diagram. Pig. I V,1 , the relative total pressure /225 
Pt,rel defined by the state : 

- it lies on tue isotherm valid in tne relative system 
according to Eq. (IV. 3) and 

- it has the same entropy as the actual final compression 
pressure Pt2*^ 

1. The same definition was used by Dean and co-worlcei's in ^[6.5], 
appendix VII). The derivation of the relative total pressure 
Pt.rel from the "relative total enthalpy relationship h^ ^ « 

(cont. on p. I8l) 
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With the isentroplc relationship, the following, equation 
is valid: 




«c/(ic-l) 


(IV.4) 


Along the isotherms pgi ■ constant , the Gibbs equatlon 
reads as follows: * 


T ds * Cp ^ - ap/p, and 
dt/R * - dp/(pRT) « - dp/p 


(IV.5) 


(TV.6) 


with s entropy; R » gas constant; p ■-density; ;<p ■ specific 
heat at constant pressure; < ■ isentroplc exponent. 

In the* inlet standard condition — p^q M» "^tO N — 
entropy Sq is set equal to zero analogous t6 Eq. T$*l8). If 
we integrate £q. (IV, 6) from state 0 to state 0*, then we 
obtain the following: 


(ttft* - 8 q)/R « 8/R * - ln<Pt,rel'PtO.N^» 


(TV. 7) 


or the following equation for tne relative total pressure 
ratio in the impeller: 


Pt,rel'PtO,N * • 


-s/R 


(IV. 8) 


constant is sensible fro;a the energy standpoint, since for an 
observer rotating with the system the impeller shaft is at 
rest, hence no work of any kind is transmitted. In the case of 
a conceivable definition of the "total pres jv re of the relative 
flow** p't rel — based on the "total enthalpy of the relative 
flow" n't*r»ei “ h + w2/2 ■ ht^rel ^ u^/ii — this condition is 
net met. *The U.S. llterature*dlcc.ingulshes — more clear'y in 
words — between p^ (relative total pi*essure) and p'^ ^gj^ 
(relative stagnatloA pressure). * 
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